A numerical method for the biharmonic problem