Convert this to a system of first order

ordinary differential equations and

propose an explicit numerical method

to solve for the vertical displacements

of the masses versus time.