Gaussian quadrature for matrix valued functions on the unit circle.
We present algorithms for the determination of polynomials orthogonal with respect to a positive weight function multiplied by a polynomial with simple roots inside the interval of integration. We apply these algorithms to search for and calculate all possible sequences of imbedded quadratures of maximal polynomials order of precision for the generalized Laguerre and Hermite weight functions.