It’s much more powerful than classic Newtonian method to solve more complex system problem, such as multiple particles and exotic coordinate systems

The result is (0.5*k*(2*L + 2*x1(t) – 2*x2(t)) + 1.0*m1*Derivative(x1(t), (t, 2)), 0.5*k*(-2*L – 2*x1(t) + 2*x2(t)) + 1.0*m2*Derivative(x2(t), (t, 2))) or

to further solve these two differential equations: