      SUBROUTINE GHCN48
C     ------------------------------------------
C     SET CONSTANTS FOR GAUSS-HERMITE QUADRATURE
C     LAST REVISION: 02/04/91 BY JCH
C     ------------------------------------------
      IMPLICIT DOUBLE PRECISION (A-H,O-Z)
C6    IMPLICIT REAL*8 (A-H,O-Z)
      DIMENSION PH(48),WH(48),WHLN(48),ZH(48)
      COMMON/LIST6/PH,WH,ZH,WHLN
      COMMON/LIST9/SR2,PI,SRPI,PIB2
C
      PH( 1)= 0.897531508193168706215358223476D1
      PH( 2)= 0.831075219070478413034784782670D1
      PH( 3)= 0.775929551976577463923770264126D1
      PH( 4)= 0.726604655416435040282062777183D1
      PH( 5)= 0.681006457807414138661792830624D1
      PH( 6)= 0.638056409618641062386710795477D1
      PH( 7)= 0.597107222501354540770712331894D1
      PH( 8)= 0.557731698122372862626895168709D1
      PH( 9)= 0.519628771879236454097986825550D1
      PH(10)= 0.482575722813320948395976669772D1
      PH(11)= 0.446401454693445890384147501246D1
      PH(12)= 0.410970460356059023650054924395D1
      PH(13)= 0.376172649022835778750118455640D1
      PH(14)= 0.341916596936388461135547839399D1
      PH(15)= 0.308124898864510584970840484976D1
      PH(16)= 0.274730862482238321948145392074D1
      PH(17)= 0.241676090487321645904643205900D1
      PH(18)= 0.208908666094427643599567785406D1
      PH(19)= 0.176381757989530000622984196672D1
      PH(20)= 0.144052522013756518650995755277D1
      PH(21)= 0.111881215240215656333400026648D1
      PH(22)= 0.798304627778562231219097373434D0
      PH(23)= 0.478646337594496098233148980980D0
      PH(24)= 0.159492935848862470507111125726D0
C
      WH( 1)= 0.793555146077399683862306435342D-35
      WH( 2)= 0.598461269331387842583125604487D-30
      WH( 3)= 0.368503608015066987944041141543D-26
      WH( 4)= 0.556457746890228483635284540092D-23
      WH( 5)= 0.318838732350513844269483296100D-20
      WH( 6)= 0.873015960118667654492535811908D-18
      WH( 7)= 0.131515962265840851189411917719D-15
      WH( 8)= 0.119758986547917935294761992879D-13
      WH( 9)= 0.704693258154588908419691347901D-12
      WH(10)= 0.281529653783816910883414782888D-10
      WH(11)= 0.793046749516538231102728561821D-9
      WH(12)= 0.162251413589576983699842464296D-7
      WH(13)= 0.246865899366975047772528943937D-6
      WH(14)= 0.284725869173484808307111697909D-5
      WH(15)= 0.252859902774848890255132137567D-4
      WH(16)= 0.175150431801172828786962777660D-3
      WH(17)= 0.956392319819415276734092656065D-3
      WH(18)= 0.415300491197755245156800170082D-2
      WH(19)= 0.144449615749810994172466181382D-1
      WH(20)= 0.404796769846038489872299218363D-1
      WH(21)= 0.918222970792851793180659667553D-1
      WH(22)= 0.169204471945641106711260999076D0
      WH(23)= 0.253961542664759097659625200207D0
      WH(24)= 0.311001030377963078884809767240D0
C
      DO 10 I=1,24
        PH(I+24)= -PH(I)
        WH(I+24)=  WH(I)
C       --------------------------------------------------------------
C       TRANSFORM TO WEIGHT FUNCTION EXP(-Z^2/2) INSTEAD OF EXP(-Z^2);
C       MUST DIVIDE INTEGRAL BY SQRT(PI) TO EFFECT INTEGRATION W.R.T.
C       STANDARD NORMAL DENSITY
C       --------------------------------------------------------------
        ZH(I)   =  SR2*PH(I)
        ZH(I+24)= -ZH(I)
        WHLN(I)=DLOG(WH(I))
        WHLN(I+24)= WHLN(I)
   10 CONTINUE
      RETURN
      END
C
