Question: Error using dsolve to solve orbital equations

I'm trying to use dsolve to get a numeric solution to a series of differential equations I have. I've got:

G := 1;

M := 1;

pT := 2;

OrbitEnergy :=  (r0) -> -1/r0+(1/2)*L^2/r0^2-L^2/r0^3;

OrbitDEs := diff(r(tau), [tau$2]) = -G*M/r(tau)^2+L^2/r(tau)^3-3*G*M*L^2/r(tau)^4, diff(phi(tau), tau) = L/r(tau)^2, diff(t(tau), tau) = sqrt(r(tau)/(r(tau)-2*G*M)+(diff(r(tau), tau))^2*r(tau)^2/(r(tau)-2*G*M)^2+r(tau)^3*(diff(phi(tau), tau))^2/(r(tau)-2*G*M));

OrbitICs := r(0) = r0, phi(0) = 0, t(0) = 0;

OrbitSolution := dsolve({OrbitDEs, OrbitICs}, numeric, {r(tau), phi(tau), t(tau)}, range = 0 .. pT)

However, I get an error message when I try and run it. "Warning, The use of global variables in numerical ODE problems is deprecated, and will be removed in a future release. Use the 'parameters' argument instead (see ?dsolve,numeric,parameters )"

Note: Right now, I've just got pT as an arbitrary number. I plan on making it variable later.

Please Wait...