Implemented the gaussianQuadrature function that can take n integrals over a function that returns a matrix. Adjusted compute_coeffs_axi for this function.
Implemented the gaussianQuadrature function that can take n integrals over a function that returns a matrix. Adjusted compute_coeffs_axi for this function.

jobovy/galpy,jobovy/galpy,jobovy/galpy,jobovy/galpy