! Gauss-Legendre quadrature roots and weights

  integer,parameter      :: nQuad = 21
  double precision, save :: rQuad(1:nQuad) = &
          (/ 0.00312391468981d0  , 0.0163865807168d0  , 0.0399503329248d0  , 0.0733183177083d0  , 0.115780018262d0   , &
             0.166430597901d0    , 0.224190582056d0   , 0.287828939896d0   , 0.355989341599d0   , 0.42721907292d0    , &
             0.5d0               , 0.57278092708d0    , 0.644010658401d0   , 0.712171060104d0   , 0.775809417944d0   , &
             0.833569402099d0    , 0.884219981738d0   , 0.926681682292d0   , 0.960049667075d0   , 0.983613419283d0   , &
             0.99687608531d0       /)
  double precision, save :: wQuad(1:nQuad) = &
          (/ 0.0080086141288872d0, 0.018476894885426d0, 0.028567212713429d0, 0.03805005681419d0 , 0.046722211728017d0, &
             0.054398649583574d0 , 0.060915708026864d0, 0.066134469316669d0, 0.069943697395537d0, 0.072262201994985d0, &
             0.07304056682485d0  , 0.072262201994985d0, 0.069943697395537d0, 0.066134469316669d0, 0.060915708026864d0, &
             0.054398649583574d0 , 0.046722211728017d0, 0.03805005681419d0 , 0.028567212713429d0, 0.018476894885426d0, &
             0.0080086141288872d0  /)

  integer,parameter      :: nAC = 21
  double precision, save :: rAC(1:nAC) = &
          (/ 0.00312391468981d0  , 0.0163865807168d0  , 0.0399503329248d0  , 0.0733183177083d0  , 0.115780018262d0   , &
             0.166430597901d0    , 0.224190582056d0   , 0.287828939896d0   , 0.355989341599d0   , 0.42721907292d0    , &
             0.5d0               , 0.57278092708d0    , 0.644010658401d0   , 0.712171060104d0   , 0.775809417944d0   , &
             0.833569402099d0    , 0.884219981738d0   , 0.926681682292d0   , 0.960049667075d0   , 0.983613419283d0   , &
             0.99687608531d0       /)
  double precision, save :: wAc(1:nAC) = &
          (/ 0.0080086141288872d0, 0.018476894885426d0, 0.028567212713429d0, 0.03805005681419d0 , 0.046722211728017d0, &
             0.054398649583574d0 , 0.060915708026864d0, 0.066134469316669d0, 0.069943697395537d0, 0.072262201994985d0, &
             0.07304056682485d0  , 0.072262201994985d0, 0.069943697395537d0, 0.066134469316669d0, 0.060915708026864d0, &
             0.054398649583574d0 , 0.046722211728017d0, 0.03805005681419d0 , 0.028567212713429d0, 0.018476894885426d0, &
             0.0080086141288872d0  /)