The combination method, which makes use of both the finite element method and the boundary element method is applied for analyses of nonlinear magnetic field. Using magnetic vector potential, the basic formulation for combining the both methods in the closed space problem is given. For treating nonlinearities arising due to material properties, the Newton-Raphson and over-relaxation schemes are used and the convergence of the solutions are examind. The results are compaired with those obtained using finite element analysis. It is shown that owing to the sparceness of the whole matrix, the convergency of the combination method is not so good as the finite element method, however, the accuracy of the converged solutions are sufficiently high especially in space region, then, this method is useful for the computation of two-dimensional nonlinear magnetic field.