clear
all
timerange=
[
0
5
]
;
%seconds
vinitial=
[
0
0
0
0
]
;
%
[
t,v
]
=ode45
(
@dof3,timerange,vinitial
)
v
(
:,
2
)
=v
(
:,
2
)
.*
180
/
pi
;
plot
(
t,v
(
:,
1
)
,t,v
(
:,
2
)
,t,v
(
:,
3
)
,t,v
(
:,
4
)
)
ylabel
(
'Rate'
)
xlabel
(
'time(s)'
)