  In this paper, the iterative method developed by Daftardar-Gejji and Jafari
(DJ method) is employed for analytic treatment of Laplace equation with
Dirichlet and Neumann boundary conditions. The method is demonstrated by
several physical models of Laplace equation. The obtained results show that the
present approach is highly accurate and requires reduced amount of calculations
compared with the existing iterative methods.
