In this paper, the numerical method of nonpolynomial spline approximation is used to solve 2nd order Neumann type boundary value problems (bvp's) in electrostatics. This new approach provides more accurate results than the ...