# dqmomo.f

          SUBROUTINE DQMOMO (ALFA, BETA, RI, RJ, RG, RH, INTEGR)
    C***BEGIN PROLOGUE  DQMOMO
    C***PURPOSE  This routine computes modified Chebyshev moments.  The K-th
    C            modified Chebyshev moment is defined as the integral over
    C            (-1,1) of W(X)*T(K,X), where T(K,X) is the Chebyshev
    C            polynomial of degree K.
    C***LIBRARY   SLATEC (QUADPACK)
    C***CATEGORY  H2A2A1, C3A2
    C***TYPE      DOUBLE PRECISION (QMOMO-S, DQMOMO-D)
    C***KEYWORDS  MODIFIED CHEBYSHEV MOMENTS, QUADPACK, QUADRATURE
    C***AUTHOR  Piessens, Robert
    C             Applied Mathematics and Programming Division
    C             K. U. Leuven
    C           de Doncker, Elise
    C             Applied Mathematics and Programming Division
    C             K. U. Leuven
    C***DESCRIPTION
    C
    C        MODIFIED CHEBYSHEV MOMENTS
    C        STANDARD FORTRAN SUBROUTINE
    C        DOUBLE PRECISION VERSION
    C
    C        PARAMETERS
    C           ALFA   - Double precision
    C                    Parameter in the weight function W(X), ALFA.GT.(-1)
    C
    C           BETA   - Double precision
    C                    Parameter in the weight function W(X), BETA.GT.(-1)
    C
    C           RI     - Double precision
    C                    Vector of dimension 25
    C                    RI(K) is the integral over (-1,1) of
    C                    (1+X)**ALFA*T(K-1,X), K = 1, ..., 25.
    C
    C           RJ     - Double precision
    C                    Vector of dimension 25
    C                    RJ(K) is the integral over (-1,1) of
    C                    (1-X)**BETA*T(K-1,X), K = 1, ..., 25.
    C
    C           RG     - Double precision
    C                    Vector of dimension 25
    C                    RG(K) is the integral over (-1,1) of
    C                    (1+X)**ALFA*LOG((1+X)/2)*T(K-1,X), K = 1, ..., 25.
    C
    C           RH     - Double precision
    C                    Vector of dimension 25
    C                    RH(K) is the integral over (-1,1) of
    C                    (1-X)**BETA*LOG((1-X)/2)*T(K-1,X), K = 1, ..., 25.
    C
    C           INTEGR - Integer
    C                    Input parameter indicating the modified
    C                    Moments to be computed
    C                    INTEGR = 1 compute RI, RJ
    C                           = 2 compute RI, RJ, RG
    C                           = 3 compute RI, RJ, RH
    C                           = 4 compute RI, RJ, RG, RH
    C
    C***REFERENCES  (NONE)
    C***ROUTINES CALLED  (NONE)
    C***REVISION HISTORY  (YYMMDD)
    C   820101  DATE WRITTEN
    C   891009  Removed unreferenced statement label.  (WRB)
    C   891009  REVISION DATE from Version 3.2
    C   891214  Prologue converted to Version 4.0 format.  (BAB)
    C***END PROLOGUE  DQMOMO