x,y,z,t=var('x y z t')
a=4
b=4
c=2
d=0.28
e=2
init_x=4
init_y=2
init_z=3
p=desolve_system_rk4([x*(a-b*y),-y*(c-d*x),c*sin(t)],[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
