[t,y] = ode45(@threelqrrhs,[0 10],[0 0 0 0 0 1]);

plot(t,y(:,1),'b');
hold on;
plot(t,y(:,3),'r');
plot(t,y(:,5),'k');
