A Legendre spectral collocation method is presented for the solutionof the biharmonic Dirichlet problem on a square. The solution andits Laplacian are approximated using the set of basis functions suggestedby Shen, which are linear combinations of Legendre polynomials. A Schurcomplement approach is used to reduce the resulting linear system to oneinvolving the approximation of the Laplacian of the solution on the twovertical sides of the square. The Schur complement system is solved bya preconditioned conjugate gradient method. The total cost of the algorithmis O(N 3). Numerical results demonstrate the spectral convergence of the method.
