Originally Posted by

**darii** I have other problem...(i do it how you say me, thanks!).

I have other differential equations:

dx = y(4)+y(2);

dy = y(5)-y(1);

dz = y(6);

dpx = y(5) - ((1-mu)/(r1.^3))*(y(1)-mu) - (mu/(r2.^3))*(y(1)-mu+1);

dpy = -y(4) -((1-mu)/(r1.^3) + mu/(r2.^3))*y(2);

dpz = -((1-mu)/(r1.^3) + mu/(r2.^3))*y(2);

and i want to resolve the system, i do:

[t,y]= ode45('test2',tf, cond_initials);

but, i want plot this..."y" have 6 columns...i want to do something like that:

plot(y(:,1),y(:,2),y(:,3),y(:,4),y(:,5),y(:,6))

but it is wrong...what can i do?

thank you very much!