Treatment of the integral divergence in computations of vacuum polarization for photons with the method of approximate values of Dirac δ function