A finite difference variational method for bending of plates