Solution (source code)

= Solution

Since
$$
\omega=\sqrt{c_s^2\kappa^2+\omega_n^2},
$$
the phase and <group> <velocities> are
$$
\boxed{
c=\frac{\omega}{\kappa}
=\sqrt{c_s^2+\frac{\omega_n^2}{\kappa^2}},
\qquad
c_g=\frac{d\omega}{d\kappa}
=\frac{c_s^2\kappa}{\sqrt{c_s^2\kappa^2+\omega_n^2}}
=\frac{c_s^2}{c}.}
$$