Gaussian quadrature for matrix valued functions on the real line

Parallel versions of the block Chebyshev algorithm to generate the recursion coefficients of orthonormal matrix polynomials and the Gaussian quadrature algorithm to approximate matrix integrals on the real line are implemented on an SP1.