Question: Subroutine to solve nonlinear DE

Hi all, I have a non-linear differential equation x := diff(y(r), r, r)+2*(diff(y(r), r))/r+ C*polylog(3/2, -exp(C2*(C3-y(r)))) = 0 with boundary values D(y)(0)=0 and y(infinity)=0. The hint says to use a subroutine to solve polylog and then integrate the differential equation with help of a Runge-Kutta like scheme. I was wondering if anyone could help me with this. I have no clue how to use subroutines and specially when y(r) is within polylog. Thanks in advance, MS
Please Wait...