\vec{p_i}=\frac{\partial L}{\partial \dot{\vec{r_i}}}