5taras's profile - activity

The last procedure doesn't plot 3D image if I take 10000 points. I get only black square.

xyzsoln=[[soln[1][0],soln[1][1],soln[1][2]] for soln in a] list_plot(xyzsoln)

thank you for your answer

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

ok, thank you. I did compile sage once but it took time.

I downloaded precompiled sage5.0 and can't make it work in natty. It requires libc6 newer than I have in the kernel. What can I do to solve the situation?