Solution (source code)

= Solution

Two polarizations and periodic wavevectors give $g(\omega)=A\omega/(\pi c_s^2)$.