Interpolation for large number of points
I have this code:
ade:= diff(y(r),r,r)+2*diff(y(r),r)/r+486*polylog(3/2,-exp(1-1/43*y(r))) = 0;
ys:= dsolve({de, y(0) = 0, D(y)(0) = 0}, y(r), series, order=10);
yn := convert(evalf(ys), polynom);yn1 := eval(yn, r = .1); dyn1 := eval(diff(yn, r));
dyn2 := eval(dyn1, r = .1);S := dsolve({de, y(.1) = 1.266393114, (D(y))(.1) = 25.18200436}, y(r), numeric);with(plots); odeplot(S, [r, y(r)], r = 0.1e-1 .. 200);
# I get a plot from this step, but I need an expression for y(r) also. A crude method is to get S(1);S(2) and so on and copy paste those values for interpolation.
Question: Would someone mind helping me with a generic code for that step?
Loading Comments & Answers
