qmomo.f
SUBROUTINE QMOMO (ALFA, BETA, RI, RJ, RG, RH, INTEGR)
C***BEGIN PROLOGUE QMOMO
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 SINGLE 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 REAL VERSION
C
C PARAMETERS
C ALFA - Real
C Parameter in the weight function W(X), ALFA.GT.(-1)
C
C BETA - Real
C Parameter in the weight function W(X), BETA.GT.(-1)
C
C RI - Real
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 - Real
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 - Real
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 - Real
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 810101 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 QMOMO