This paper applies the indirect method - potential theory for converting the partial differential equation to the boundary integral equation for harmonic transformation model. Potential theory (PT) is simple and easy to implement comparing with Green's theorem. PT is designed specifically for integral equation. It has been verified that it is efficiency in detecting the singularity point to the harmonic transformation. Experimental results show that the potential theory has good performance comparing with the traditional one.