Abstract
When standard spectral methods are applied to the biharmonic equation the condition number of the resulting algebraic system of equations is very large and grows as 0(N
8
) where N is the degree of the polynomial expansion in each coordinate direction. When direct methods are used to solve these systems we fail to achieve machine precision for the solution due to rounding errors. An improved method which reduces the condition number to 0(N
4
) is described. Numerical results are presented demonstrating the exponential convergence of the approximation.