$\psi_{pq}(x,y)=X_p(x)\,Y_q(y)$
$ \omega_{pq}=\sqrt{D/\rho h}\,(\beta_p^2+\gamma_q^2) $
$$ Y(\omega)=\frac{i\omega}{\rho h L_x L_y}\sum_p\sum_q \frac{X_p^2(x_0)\,Y_q^2(y_0)\,L_xL_y}{N_xN_y\,[\omega_{pq}^2(1+i\eta)-\omega^2]} $$