Due to orbiting [1]}, [2]}, [3]}, [4]},
the integration must be divided into two parts when the energy is less
than \(E^{\star }_{\mathrm {C}}\) , corresponding to \(r^{\star }_{\mathrm {C}}\) ,
which is the minimum of the function
\(y(r^{\star })=V^{\star }(r^{\star })+\tfrac{1}{2}r^{\star }{V^{\star }(r^{\star })}{r^{\star }}.\) 
