A new fast and accurate protection zone determination methodology for lightning rod is presented. The methodology is based on Quasi Monte Carlo Integration technique applied to Method of Moments (MoM) solution of Integral Equations. As an example, solution of the integral equation for unknown charge distribution on lightning rod is obtained. The electric field in the region surrounding the rod is then computed and the protection zone plotted accordingly. It is observed that the proposed method is capable of dealing with the problem of singularity efficiently.