By using the generalized variational principles of thin plate bending problem, a numerical method is proposed for a mixed boundary value problem of circular plate.