GQLOG: Matlab routines for computing Gauss Quadrature rules with logarithmic weight functions

<p>By means of the Matlab symbolic/variable-precision facilities, routines are developed that generate an arbitrary number of recurrence coefficients to any given precision for polynomials orthogonal with respect to weight functions of Laguerre and Jacobi type containing logarithmic factors. The results are applied to Gaussian quadrature of integrals involving weight functions of the type mentioned.</p>