Solution (source code)

= Solution

The initial vector $(d,0)$ has coefficient $d/\sqrt2$ in each normalized <normal mode>. Zero initial velocities remove their sine terms. Therefore
$$
\boxed{x_1(t)=\frac d2(\cos\omega_-t+\cos\omega_+t),\qquad
x_2(t)=\frac d2(\cos\omega_-t-\cos\omega_+t).}
$$
At zero these give $(d,0)$ and both derivatives vanish; substitution into the coupled equations confirms the opposite signs of the two higher-frequency components.