Hello, 

I am trying to analyse a nonlinear ode. To see the behaviour,

I am using implicit plot (suggested by @Preben Alsholm).

The problem is that every time I need to search for the value of N

(=infinity) to have some sort of graph (needed for analysis).

The other thing is, if we have an output of the implicitplot, How we can be sure

that it is the correct one?

restart:with(plots):with(plottools):with(stats):
unprotect(D):
a[1]:=0:a[2]:=0.2:a[3]:=-0.15:a[4]:=0.85:
b[1]:=0:b[2]:=-0.26:b[3]:=0.28:b[4]:=0.04:
c[1]:=0:c[2]:=0.23:c[3]:=0.26:c[4]:=-0.04:
d[1]:=0.16:d[2]:=0.22:d[3]:=0.24:d[4]:=0.85:
e[1]:=0:e[2]:=0:e[3]:=0:e[4]:=0:
f[1]:=0:f[2]:=0.2:f[3]:=0.2:f[4]:=0.2:
p[1]:=0.01:p[2]:=0.08:p[3]:=0.15:
N:=10000:x[0]:=0:y[0]:=0:
randomize():
for n from 0 to N do
r[n]:=rand()/(10)^(12);
if r[n]then i:=2 elif r[n...

Does any one has any idea about "Ndynamics" package?

Whenever, I try to call it, maple gives this error

restart:

with(Ndynamics);

"Error, invalid input: with expects its 1st argument, pname, to be of type {`module`, package},

but received Ndynamics"

Is this an old package and replaced by a new one?

Thanks

Hello every one,

I got an error while trying to solve numerically a nonlinear system of odes using 

the dsolve command. This error is very common,

"Error, (in dsolve/numeric/bvp) initial Newton iteration is not converging".

I tried the different technique explained in 

http://www.maplesoft.com/support/help/Maple/view...

This question was asked before but because of the curiosity, I 

bring it in the light again. We have a system of odes

restart:with(plots):

a := 1; b := .5; d := 1; omega := .4; h1 := 1+a*cos(x); h2 := -d-b*cos(x+omega);

F := Q-1-d;

de:={alpha*(diff(f(y), y, y, y, y))+G*(diff(theta(y), y, y))+B*(diff(phi(y), y, y))

+6*beta*(diff(f(y), y, y))*(diff(f(y), y, y, y))^2+3*beta*(diff(f(y), y, y, y, y))*

(diff(f(y), y, y))^2 = 0,