Numerical solution of two dimensional harmonic boundary problems containing singularities by conformal transformation methods