Solution (source code)

= Solution

Fix $u=u_0$. Initial data $p(u_0),p'(u_0)$ and $q(u_0),q'(u_0)$ may be prescribed independently for the two oscillator equations. In particular, values of $p(u_0)$ and $q(u_0)$ generate arbitrary translations in $x$ and $y$, while their derivatives merely add controllable $w$ components. The independent Killing field $\partial_w$ supplies any remaining translation in $w$. These Killing fields span the tangent space of each constant-$u$ wavefront, so their flows act transitively. Therefore an isometry maps any point $r$ on such a surface to any other point $s$, realizing a <homogeneous wavefront of a plane gravitational wave>.