Ask Your Question

Revision history [back]

click to hide/show revision 1
initial version

how to plot solution of differential system

Hi, everyone! I got a question. I have a system of three Newton's equations x''=-x/(xx+yy+zz)-F(t) y''=-x/(xx+yy+zz) z''=-x/(xx+yy+z*z)

I use odesolve() to solve numerically the system by using y[0] for x, y[1] for y, y[2] for z like in example in the manual without using jacobian: j_1= lambda t,y: [ y[3], y[4], y[5], -y[0]/(y[0]y[0]+y[1]y[1]+y[2]y[2])^(1.5)-Fcos(wt), -y[1]/(y[0]y[0]+y[1]y[1]+y[2]y[2])^(1.5), -y[2]/(y[0]y[0]+y[1]y[1]+y[2]*y[2])^(1.5)]

How may I plot the curve (x(t), y(t)) or (x(t), z(t))? T.plot_solution(i=0, filename=outfile) allows to plot only x(t), y(t) or z(t) Thank you

how to plot solution of differential system

Hi, everyone! I got a question. I have a system of three Newton's equations x''=-x/(xx+yy+zz)-F(t) z)^3/2-F(t) y''=-x/(xx+yy+zz) z)^3/2 z''=-x/(xx+yy+z*z)y+z*z)^3/2

I use odesolve() to solve numerically the system by using y[0] for x, y[1] for y, y[2] for z like in example in the manual without using jacobian: j_1= lambda t,y: [ y[3], y[4], y[5], -y[0]/(y[0]y[0]+y[1]y[1]+y[2]y[2])^(1.5)-Fcos(wt), -y[1]/(y[0]y[0]+y[1]y[1]+y[2]y[2])^(1.5), -y[2]/(y[0]y[0]+y[1]y[1]+y[2]*y[2])^(1.5)]

How may I plot the curve (x(t), y(t)) or (x(t), z(t))? T.plot_solution(i=0, filename=outfile) allows to plot only x(t), y(t) or z(t) Thank you

how to plot solution of differential system

Hi, everyone! I got a question. I have a system of three Newton's equations x''=-x/(xx+yy+zz)^3/2-F(t) y''=-x/(xx+yy+zz)^3/2 z''=-x/(xx+yy+z*z)^3/2

I use odesolve() to solve numerically the system by using y[0] for x, y[1] for y, y[2] for z like in example in the manual without using jacobian: j_1= lambda t,y: [ y[3], y[4], y[5], -y[0]/(y[0]y[0]+y[1]y[1]+y[2]y[2])^(1.5)-Fcos(wt), -y[1]/(y[0]y[0]+y[1]y[1]+y[2]y[2])^(1.5), -y[2]/(y[0]y[0]+y[1]y[1]+y[2]*y[2])^(1.5)]

How may I plot the curve (x(t), y(t)) or (x(t), z(t))? T.plot_solution(i=0, filename=outfile) allows to plot only x(t), y(t) (x(t), t), (y(t), t) or z(t) (z(t), t) Thank you

click to hide/show revision 4
No.4 Revision

updated 2012-09-21 16:44:12 -0600

kcrisman gravatar image

how to plot solution of differential system

Hi, everyone! I got a question. I have a system of three Newton's equations x''=-x/(xx+yy+zz)^3/2-F(t) y''=-x/(xx+yy+zz)^3/2 z''=-x/(xx+yy+z*z)^3/2equations

x''=-x/(x*x+y*y+z*z)^3/2-F(t)
y''=-x/(x*x+y*y+z*z)^3/2
z''=-x/(x*x+y*y+z*z)^3/2

I use odesolve() odesolve() to solve numerically the system by using y[0] y[0] for x, y[1] y[1] for y, y[2] y[2] for z like in example in the manual without using jacobian: jacobian:

j_1= lambda t,y: [ y[3], y[4], y[5],  -y[0]/(y[0]y[0]+y[1]y[1]+y[2]y[2])^(1.5)-Fcos(wt), -y[1]/(y[0]y[0]+y[1]y[1]+y[2]y[2])^(1.5), -y[2]/(y[0]y[0]+y[1]y[1]+y[2]*y[2])^(1.5)]

-y[0]/(y[0]*y[0]+y[1]*y[1]+y[2]*y[2])^(1.5)-F*cos(w*t), -y[1]/(y[0]*y[0]+y[1]*y[1]+y[2]*y[2])^(1.5), -y[2]/(y[0]*y[0]+y[1]*y[1]+y[2]*y[2])^(1.5)]

How may I plot the curve (x(t), y(t)) or (x(t), z(t))? T.plot_solution(i=0, filename=outfile) filename=outfile) allows to plot only (x(t), t), (y(t), t) or (z(t), t) Thank you

click to hide/show revision 5
No.5 Revision

how to plot solution of differential system

Hi, everyone! I got a question. I have a system of three Newton's equations

x''=-x/(x*x+y*y+z*z)^3/2-F(t)
y''=-x/(x*x+y*y+z*z)^3/2
z''=-x/(x*x+y*y+z*z)^3/2

I use odesolve()ode_solver() to solve numerically the system by using y[0] for x, y[1] for y, y[2] for z like in example in the manual without using jacobian:

j_1= lambda t,y: [ y[3], y[4], y[5],  -y[0]/(y[0]*y[0]+y[1]*y[1]+y[2]*y[2])^(1.5)-F*cos(w*t), -y[1]/(y[0]*y[0]+y[1]*y[1]+y[2]*y[2])^(1.5), -y[2]/(y[0]*y[0]+y[1]*y[1]+y[2]*y[2])^(1.5)]

How may I plot the curve (x(t), y(t)) or (x(t), z(t))? T.plot_solution(i=0, filename=outfile) allows to plot only (x(t), t), (y(t), t) or (z(t), t) Thank you