Numerical implementation of the boundary element method with point-source approximation of the potential