x,y,t=var('x y t')
x=function('x',t)
P=desolve_system_rk4([x(1-y),-y(1-x)],[x,y],ics=[0,0.5,2],ivar=t,end_points=20) P=P+desolve_system_rk4([x(1-y),-y(1-x)],[x,y],ics=[0,0.3,2],ivar=t,end_points=20) Q=[ [i,j] for i,j,k in P] LP=list_plot(Q)
Q=[ [j,k] for i,j,k in P]
Q=[ [j,k] for i,j,k in P] LP=list_plot(Q) show(LP) plot(x);plot(y) #Doesn't work