L=\frac{1}{2}m\dot{\mathbf{r}}^2-U(r)
