The objective of this work is to extend the use of the Analytical Method for assessing the interference involving non-geostationary satellite networks to situations where the interference entries from the various satellites in a given Constellations are modeled as statistically independent random variables. The analysis of this type of situation requires a large number of convolution computations. After showing that the direct approach to numerically compute convolutions is not feasible and that the use of Monte Carlo simulation techniques require a prohibitive computer time to guarantee the necessary accuracy, the study identifies two alternatives: the indirect use of the Central Limit theorem, and numerical integration using Gauss Quadrature Rules. Finally, the present work presents two numerical example applications of the Analytical Method. The obtained results show an accuracy improvement when compared to those obtained using other methods.