Galerkin Boundary Element Method(GBEM) is a powerful technique to solve partial differential equations. In this work the method is presented by solving the axisymmetric Laplace's equation. Later the Grad-Shafranov Equation(GSE) is solved by means of this powerful technique. The GSE is a simplification of ideal MHD in two dimensions. The resulting equation is a nonlinear, elliptic partial differential equation. The GBEM is well suited for solving this type of equations in 'irregular domains' such as unbounded and irregular shape domains.