Question: ODE plot not working

Don't get a DEplot 


 

with(DEtools,DEplot);with(plots,display);

[DEplot]

 

[display]

(1)

 

 

 

init3:=seq([0,k*0.2,0],k=-4..10);

[0, -.8, 0], [0, -.6, 0], [0, -.4, 0], [0, -.2, 0], [0, 0., 0], [0, .2, 0], [0, .4, 0], [0, .6, 0], [0, .8, 0], [0, 1.0, 0], [0, 1.2, 0], [0, 1.4, 0], [0, 1.6, 0], [0, 1.8, 0], [0, 2.0, 0]

(2)

 

The forward plot is p3:

 

p3:=DEplot([sys],[x,y],t=0..2,[init3],x=-3..2,y=-3..2,stepsize=.05,arrows=none,linecolor=black,thickness=1,axes=box,scaling=constrained):

Error, (in f) numeric exception: division by zero

 

he backward plot is np3:

 

np3:=DEplot([sys],[x,y],t=-2..0,[init3],x=-3..2,y=-3..2,stepsize=.05,arrows=none,linecolor=black,thickness=1,axes=box,scaling=constrained):

Warning, plot may be incomplete, the following errors(s) were issued:
   cannot evaluate the solution further left of -.29861232, probably a singularity
Warning, plot may be incomplete, the following errors(s) were issued:
   cannot evaluate the solution further left of -.93147185e-1, probably a singularity
Warning, plot may be incomplete, the following errors(s) were issued:
   cannot evaluate the solution further left of -.23648931e-1, probably a singularity
Warning, plot may be incomplete, the following errors(s) were issued:
   cannot evaluate the solution further left of -.27325542e-2, probably a singularity

 

Error, (in f) numeric exception: division by zero

 

 

display({p1,p2,p3,np3});

Error, (in plots:-display) expecting plot structures but received: {np3}

 

 


 

Download uitzoeken_ode_lesmateriaal_betounes_of_het_werkt.mw

Please Wait...