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