The semianalytical bolometric method outlined in Ref. [27], where we assume equal energy in neutrinos as in gamma rays and equate the power fluxes, allows us to determine the neutrino flux using: $$\int E_{\gamma }\frac{d{N}_{\gamma }}{d{E}_{\gamma }}d{E}_{\gamma }=\int E_{\nu }\frac{d{N}_{\nu }}{d{E}_{\nu }}d{E}_{\nu }$$While this method will determine the normalization of the neutrino spectrum, we must posit a priori a slope $s_{\nu }$ (assuming we take the neutrino spectrum as a simple power law) as well as an upper threshold. In this case we take $s_{\nu }=-2.0$ and $E_{\nu }^{\max }=50\mathrm{TeV}$, recognizing that in principle we cannot assume a pionic origin for the gamma rays or that the sequence of events that we have chosen to have evolve the spectra is the correct one. The main consequence of using this method is that the variation in energy in the $\gamma$-ray flux at the EGRET energy range causes a large variation in the harder neutrino flux at GeV–TeV energies, potentially greatly increasing the number of detected neutrinos per year. By this method, we find event rates between 2 and 12 neutrinos/year, giving, respectively, 4.9 and 29.4 $\sigma$ after 15 years of operation.