so if [itex]\alpha = 2E/\mu[/itex] and [itex]\beta = L^2\alpha^2/\mu^2[/itex] (not the same alpha, sorry) and bounds [r0,r] I should get:
[tex]\frac{\mu}{2E} \left[\left(\frac{2E}{\mu}r^2 + \frac{L^2\alpha^2 }{\mu^2}\right)^{1/2} - \left(\frac{2E}{\mu}r_0^2 + \frac{L^2\alpha^2 }{\mu^2}\right)^{1/2}\right][/tex]
and this is equal to time so solving r(t) gives a quadratic. Does this make sense for central force motion where [itex]r = ke^{-\alpha\theta}[/itex]?