We propose a spectral collocation method to approximate the exact boundary control of the wave equation in a square domain. The idea is to introduce a suitable approximate control problem that we solve in the finite-dimensional space of polynomials of degree \(N\in {\mathbb {N}}\) in space. We prove that we can choose a sequence of controls \(f^N\) associated with the approximate control problem in such a way that they converge, as \(N\rightarrow \infty \) , to a control of the continuous wave equation. Unlike other numerical approximations tried in the literature, this one does not require regularization techniques and can be easily adapted to other equations and systems where the controllability of the continuous model is known. The method is illustrated with several examples in 1-d and 2-d in a square domain. We also give numerical evidence of the highly accurate approximation inherent to spectral methods.
