The authors propose a numerical solution to the eigenvalue problem for the biharmonic operator and the bending plate problem using a spectral Legendre-Galerkin method. Using quadruple-precision arithmetic, accurate computations of the nodal curves are obtained for the principal eigenfunction near the corners of the plate. Comparisons between known asymptotic expansions of the solution near a corner provide evidence of an accurate computation of the eigenfunction over the entire domain.