Ignore

Preben Alsholm — Tue, 06 Nov 2012 20:02:20 Z

I tried the following in Maple 16:

restart;
Digits:=40:
de := proc (omega, l, M) options operator, arrow; diff(y(r), `$`(r, 2))+(2*r-2*M)*(diff(y(r), r))/(r*(r-2*M))+(omega^2*r^2/(r-2*M)^2-l*(l+1)/(r*(r-2*M)))*y(r) = 0 end proc;
init:=y(10000)=-.00008944+.00004472*I;
dinit:=D(y)(10000)=-4.4641e-6-8.9e-6*I;
sol:=dsolve([de(1/10, 1, 1), init, dinit], numeric, range = 10 .. 10000);
plots:-odeplot(sol,[Re(y(r)),Im(y(r))],10..10000,refine=1);

The plot looks fine. Since init and dinit in your posted question are shown in an image, I didn't use all the digits.

I think if you did use all the digits you

fpghost — Wed, 07 Nov 2012 14:25:18 Z

I think if you did use all the digits you would recover the same error. I have no idea why the init and dinit keep copying as images. Maybe just make up the remaininder of the 40dps.

MaplePrimes - answers and comments on Question, Solving numerically ode errors

Ignore

I think if you did use all the digits you