x,y,z,t=var('x y z t')
a=10
b=28
c=8/3
init_x=2
init_y=3
init_z=4
p=desolve_system_rk4([a*(y-x),x*(b-z)-y,x*y-c*z],[x,y,z],ics=[0,init_x,init_y,init_z],ivar=t,end_points=20)
q1=[[i,j] for i,j,k,1 in p]
lp1=list_plot(q1,plotjoined=true)
q2=[[i,k] for i,j,k,1 in p]
lp2=list_plot(q2,plotjoined=true, color='red')
q3=[[i,1] for i,j,k,1 in p]
lp3=list_plot(q3,plotjoined=true, color='green')
lp1+lp2+lp3