Write
r1=(x1,y1),
r2=(x2,y2),
ρ1=∣r1∣, and
ρ2=∣r2−r1∣. Then
L=2m1∣r˙1∣2+2m2∣r˙2∣2−2k1(ρ1−ℓ)2−2k2(ρ2−ℓ)2−m1gy1−m2gy2.
With
ey=(0,1), the Euler--Lagrange equations are
m1r¨1=−k1(1−ρ1ℓ)r1+k2(1−ρ2ℓ)(r2−r1)−m1gey,
m2r¨2=−k2(1−ρ2ℓ)(r2−r1)−m2gey.
Solved by gpt-5.6-sol high.