|
|
# from sage.calculus.desolvers import desolve_odeint
x,y=var('x,y')
mu= -1.5 # de trage richting in knooppunt (mu,0) met mu>0 verandert rond mu=1
send=3.
f=[mu*x-x^2,-y]
p1=plot_vector_field(f,(x,-2.5,2),(y,-2.5,2.5),color='black',headaxislength=5)
sol=desolve_odeint(f,[mu,2],srange(0,100,0.1),[x,y])
p2=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[mu,-2],srange(0,100,0.1),[x,y])
p2n=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[0,2],srange(0,100,0.1),[x,y])
p3=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[0,-2],srange(0,100,0.1),[x,y])
p3n=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[mu/2,2],srange(0,100,0.1),[x,y])
p4=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[mu/2,-2],srange(0,send,0.1),[x,y])
p4n=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[mu-0.1,2],srange(0,1.3,0.1),[x,y])
p5=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[mu-0.1,-2],srange(0,1.3,0.1),[x,y])
p5n=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[-mu,2],srange(0,send,0.1),[x,y])
p6=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[-mu,-2],srange(0,send,0.1),[x,y])
p6n=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[-mu,2],srange(0,send,0.1),[x,y])
p7=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[mu+0.1,2],srange(0,send,0.1),[x,y])
p8=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[mu+0.1,-2],srange(0,send,0.1),[x,y])
p8n=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[mu+0.1,0],srange(0,send,0.1),[x,y])
p9=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[-mu+0.1,0],srange(0,send,0.1),[x,y])
p10=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[mu-0.1,0],srange(0,1.3,0.1),[x,y])
p11=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
#show(p1+p2+p2n+p3+p3n+p4+p4n+p5+p5n+p6+p6n+p8+p8n+p9+p10+p11)
show(p1+p2+p2n+p3+p3n+p4+p4n+p5+p5n+p6+p6n+p8+p8n+p9+p10+p11, ticks=[[],[]],axes=false)
|
# from sage.calculus.desolvers import desolve_odeint
x,y=var('x,y')
mu= -0.5 # de trage richting in knooppunt (mu,0) met mu>0 verandert rond mu=1
send=3.
f=[mu*x-x^2,-y]
p1=plot_vector_field(f,(x,-0.75,0.5),(y,-2,2),color='black',headaxislength=5)
sol=desolve_odeint(f,[mu,2],srange(0,100,0.1),[x,y])
p2=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[mu,-2],srange(0,100,0.1),[x,y])
p2n=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[0,2],srange(0,100,0.1),[x,y])
p3=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[0,-2],srange(0,100,0.1),[x,y])
p3n=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[mu/2,2],srange(0,100,0.1),[x,y])
p4=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[mu/2,-2],srange(0,send,0.1),[x,y])
p4n=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[mu-0.1,2],srange(0,1,0.1),[x,y])
p5=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[mu-0.1,-2],srange(0,1,0.1),[x,y])
p5n=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[-mu,2],srange(0,send,0.1),[x,y])
p6=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[-mu,-2],srange(0,send,0.1),[x,y])
p6n=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[-mu,2],srange(0,send,0.1),[x,y])
p7=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[mu+0.1,2],srange(0,send,0.1),[x,y])
p8=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[mu+0.1,-2],srange(0,send,0.1),[x,y])
p8n=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[mu+0.1,0],srange(0,send,0.1),[x,y])
p9=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[-mu,0],srange(0,send,0.1),[x,y])
p10=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[mu-0.1,0],srange(0,1,0.1),[x,y])
p11=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
show(p1+p2+p2n+p3+p3n+p4+p4n+p5+p5n+p6+p6n+p8+p8n+p9+p10+p11, ticks=[[],[]],axes=false)
|
# from sage.calculus.desolvers import desolve_odeint
x,y=var('x,y')
mu= 0.5 # de trage richting in knooppunt (mu,0) met mu>0 verandert rond mu=1
send=10
f=[mu*x-x^2,-y]
p1=plot_vector_field(f,(x,-0.5,1),(y,-2,2),color='black',headaxislength=5)
sol=desolve_odeint(f,[mu,2],srange(0,100,0.1),[x,y])
p2=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[mu,-2],srange(0,100,0.1),[x,y])
p2n=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[0,2],srange(0,100,0.1),[x,y])
p3=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[0,-2],srange(0,100,0.1),[x,y])
p3n=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[mu/2,2],srange(0,100,0.1),[x,y])
p4=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[mu/2,-2],srange(0,send,0.1),[x,y])
p4n=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[mu-0.1,2],srange(0,send,0.1),[x,y])
p5=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[mu-0.1,-2],srange(0,send,0.1),[x,y])
p5n=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[-mu/2,2],srange(0,1,0.1),[x,y])
p6=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[-mu/2,-2],srange(0,1,0.1),[x,y])
p6n=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[-mu,2],srange(0,send,0.1),[x,y])
p7=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[mu+0.4,2],srange(0,send,0.1),[x,y])
p8=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[mu+0.4,-2],srange(0,send,0.1),[x,y])
p8n=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[mu+0.5,0],srange(0,send,0.1),[x,y])
p9=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[-mu/2,0],srange(0,1,0.1),[x,y])
p10=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[0.1,0],srange(0,10,0.1),[x,y])
p11=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
show(p1+p2+p2n+p3+p3n+p4+p4n+p5+p5n+p6+p6n+p8+p8n+p9+p10+p11, ticks=[[],[]],axes=false)
lsoda-- at current t (=r1), mxstep (=i1) steps
taken on this call before reaching tout
in above message, i1 = 500
in above message, r1 = 0.1386294039478D+01
lsoda-- at current t (=r1), mxstep (=i1) steps
taken on this call before reaching tout
in above message, i1 = 500
in above message, r1 = 0.1386294039478D+01
|
# from sage.calculus.desolvers import desolve_odeint
x,y=var('x,y')
mu=1.5 # de trage richting in knooppunt (mu,0) met mu>0 verandert rond mu=1
send=10
ss=10
f=[mu*x-x^2,-y]
p1=plot_vector_field(f,(x,-1,2.5),(y,-2,2),color='black',headaxislength=5)
sol=desolve_odeint(f,[mu,2],srange(0,ss,0.1),[x,y])
p2=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[mu,-2],srange(0,ss,0.1),[x,y])
p2n=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[0,2],srange(0,ss,0.1),[x,y])
p3=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[0,-2],srange(0,ss,0.1),[x,y])
p3n=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[mu/2,2],srange(0,ss,0.1),[x,y])
p4=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[mu/2,-2],srange(0,send,0.1),[x,y])
p4n=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[mu-0.1,2],srange(0,send,0.1),[x,y])
p5=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[mu-0.1,-2],srange(0,send,0.1),[x,y])
p5n=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[-mu/4,2],srange(0,0.6,0.1),[x,y])
p6=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[-mu/4,-2],srange(0,0.6,0.1),[x,y])
p6n=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[-mu,2],srange(0,send,0.1),[x,y])
p7=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[mu+0.8,2],srange(0,send,0.1),[x,y])
p8=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[mu+0.8,-2],srange(0,send,0.1),[x,y])
p8n=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[mu+0.8,0],srange(0,send,0.1),[x,y])
p9=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[-mu/2,0],srange(0,0.25,0.1),[x,y])
p10=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
sol=desolve_odeint(f,[0.1,0],srange(0,ss,0.1),[x,y])
p11=line(zip(sol[:,0],sol[:,1]),color='blue',thickness=2)
show(p1+p2+p2n+p3+p3n+p4+p4n+p5+p5n+p6+p6n+p8+p8n+p9+p10+p11, ticks=[[],[]],axes=false)
lsoda-- at current t (=r1), mxstep (=i1) steps
taken on this call before reaching tout
in above message, i1 = 500
in above message, r1 = 0.4620980364892D+00
lsoda-- at current t (=r1), mxstep (=i1) steps
taken on this call before reaching tout
in above message, i1 = 500
in above message, r1 = 0.4620980364892D+00
|
|
|
|
|