This paper presents boundary element methods for the solution of elliptic partial differential equation with Neumann boundary problem in R~2. Total details of error analysis are given for general condition, and optimal estimates are obtained. For simplicity Laplace's equation is discussed in illustration.