Advanced Search
All Posts : MaplePrimes Commons General Technical Discussions

numsolve ode

franky 0
May 01 2006
0
I would like to solve the following problem as stated below:

> ode3:=J=diff(E(x),x)*E(x)-diff((diff(E(x),x)),x);


> dsolve(ode3);



> odeplot(dsolve({ode3,D(E)(0)=2,D(E)(1)=2000},numeric,range=0..1),refine=1);

Error, (in dsolve/numeric/BVPSolve) initial Newton iteration is not converging

I have taken "J=2" doesn't matter too much for now. As you will see from the second command Maple can find an analytical solution to the problem. From my trials it seems that the initial condition D(E)(1)=2000 makes the problem, smaller values are ok. I have tried "Digits:=100" but this didn't help me. In the end it would be nice to go to larger (about 7 orders of magnitude) values for D(E)(1) in comparison to D(E)(0). Can you help me find a way to solve this problem?

I have maple 9.0

Best regards,

Frank
plot

Please Wait...