 Sine collocation of square modified Helmholtz global relations

On the square $[-1,1{]}^2$, write normal traces in ${\varphi }_n(s)=\sin [n\pi (s+1)/2]$. Set $p_n=n\pi /2$, ${\omega }_n=\sqrt{k^2+p_n^2}$, $r_n=({\omega }_n+p_n)/k$. Paired spectral samples $\{r_n,r_n^{-1}\}$, $\{-r_n,-r_n^{-1}\}$, $\{i{r}_n,i/r_n\}$, $\{-i{r}_n,-i/r_n\}$ give adjoint tests $e^{\mp {\omega }_{n}y}{\varphi }_n(x)$ and $e^{\mp {\omega }_{n}x}{\varphi }_n(y)$. Each test vanishes on the adjacent sides, removing their unknown normal derivatives.