Hi!
I tried to compute an ideal barbell-shaped object's dynamics, but my results were wrong.
My Langrangian is:
## L = \frac{m}{2} ( \dot{x_1}^2 + \dot{x_2}^2 + \dot{y_1}^2 + \dot{y_2}^2 ) - U( x_1 , y_1 ) - U ( x_2 , y_2 ) ##
And the constraint is:
## f = ( x_1 - x_2 )^2 + ( y_1 -...