This is currently only for the isotropic cone model. A basic and rough Monte Carlo integration
method has been hand coded to replace the advanced scipy.integrate.tplquad() quadratic integration
function. Not only are the number of integration steps massively decreased, but all PCS values can
be integrated simultaneously. I.e. there is one iteration of the integration algorithm for all PCS
values of all atoms and alignments, rather than one iteration per PCS.
