Advanced Search

Question: Interpolation for large number of points

Question:Interpolation for large number of points

manjees 360 Maple
convert plot
June 26 2009
0
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?
Please Wait...