\phi=\int\frac{M\dd{r}/r^2}{\sqrt{\frac{2}{m}(E-U(r))-\frac{M^2}{m^2r^2}}}+\text{constant}