Google
 

Trailing-Edge - PDP-10 Archives - decus_20tap2_198111 - decus/20-0026/dqh8.ssp
There are 2 other files named dqh8.ssp in the archive. Click here to see a list.
C                                                                       DQH8  10
C     ..................................................................DQH8  20
C                                                                       DQH8  30
C        SUBROUTINE DQH8                                                DQH8  40
C                                                                       DQH8  50
C        PURPOSE                                                        DQH8  60
C           TO COMPUTE INTEGRAL(EXP(-X*X)*FCT(X), SUMMED OVER X FROM    DQH8  70
C                               -INFINITY TO +INFINITY).                DQH8  80
C                                                                       DQH8  90
C        USAGE                                                          DQH8 100
C           CALL DQH8 (FCT,Y)                                           DQH8 110
C           PARAMETER FCT REQUIRES AN EXTERNAL STATEMENT                DQH8 120
C                                                                       DQH8 130
C        DESCRIPTION OF PARAMETERS                                      DQH8 140
C           FCT    - THE NAME OF AN EXTERNAL DOUBLE PRECISION FUNCTION  DQH8 150
C                    SUBPROGRAM USED.                                   DQH8 160
C           Y      - THE RESULTING DOUBLE PRECISION INTEGRAL VALUE.     DQH8 170
C                                                                       DQH8 180
C        REMARKS                                                        DQH8 190
C           NONE                                                        DQH8 200
C                                                                       DQH8 210
C        SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED                  DQH8 220
C           THE EXTERNAL DOUBLE PRECISION FUNCTION SUBPROGRAM FCT(X)    DQH8 230
C           MUST BE FURNISHED BY THE USER.                              DQH8 240
C                                                                       DQH8 250
C        METHOD                                                         DQH8 260
C           EVALUATION IS DONE BY MEANS OF 8-POINT GAUSSIAN-HERMITE     DQH8 270
C           QUADRATURE FORMULA, WHICH INTEGRATES EXACTLY WHENEVER       DQH8 280
C           FCT(X) IS A POLYNOMIAL UP TO DEGREE 15.                     DQH8 290
C           FOR REFERENCE, SEE                                          DQH8 300
C           SHAO/CHEN/FRANK, TABLES OF ZEROS AND GAUSSIAN WEIGHTS OF    DQH8 310
C           CERTAIN ASSOCIATED LAGUERRE POLYNOMIALS AND THE RELATED     DQH8 320
C           GENERALIZED HERMITE POLYNOMIALS, IBM TECHNICAL REPORT       DQH8 330
C           TR00.1100 (MARCH 1964), PP.213-214.                         DQH8 340
C                                                                       DQH8 350
C     ..................................................................DQH8 360
C                                                                       DQH8 370
      SUBROUTINE DQH8(FCT,Y)                                            DQH8 380
C                                                                       DQH8 390
C                                                                       DQH8 400
      DOUBLE PRECISION X,Y,Z,FCT                                        DQH8 410
C                                                                       DQH8 420
      X=.29306374202572440D1                                            DQH8 430
      Z=-X                                                              DQH8 440
      Y=.19960407221136762D-3*(FCT(X)+FCT(Z))                           DQH8 450
      X=.19816567566958429D1                                            DQH8 460
      Z=-X                                                              DQH8 470
      Y=Y+.17077983007413475D-1*(FCT(X)+FCT(Z))                         DQH8 480
      X=.11571937124467802D1                                            DQH8 490
      Z=-X                                                              DQH8 500
      Y=Y+.20780232581489188D0*(FCT(X)+FCT(Z))                          DQH8 510
      X=.38118699020732212D0                                            DQH8 520
      Z=-X                                                              DQH8 530
      Y=Y+.66114701255824129D0*(FCT(X)+FCT(Z))                          DQH8 540
      RETURN                                                            DQH8 550
      END                                                               DQH8 560