Abstract: Numerical method for solving the free boundary plasma equilibrium equation (a nonlinear partial differential equation of elliptic type) during the evolution of the external maintaining, field is presented. The direct method is employed to solve the corresponding set of difference equation and compared with the over relaxation method.It is shown that the direct method can give better computation precision and can save considerable computer time.