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)')