This article presents a method which models printed circuit antennas with very high efficiency. The model is based on a decomposition of the pertinent Green's function, using weighted image theory, into two smooth-valued perturbation terms and a singularity containing term which can be treated analytically. This, in addition to the adoptation of a numerical interpolation technique, leads to minimization of computation time when the method of moments is enforced. An analytic-numeric real axis integration approach is also applied to evaluate the near field quantity of the perturbation terms.