Google
 

Trailing-Edge - PDP-10 Archives - decuslib10-02 - 43,50145/qh5.ssp
There are 2 other files named qh5.ssp in the archive. Click here to see a list.
C                                                                       QH5   10
C     ..................................................................QH5   20
C                                                                       QH5   30
C        SUBROUTINE QH5                                                 QH5   40
C                                                                       QH5   50
C        PURPOSE                                                        QH5   60
C           TO COMPUTE INTEGRAL(EXP(-X*X)*FCT(X), SUMMED OVER X FROM    QH5   70
C                               -INFINITY TO +INFINITY).                QH5   80
C                                                                       QH5   90
C        USAGE                                                          QH5  100
C           CALL QH5 (FCT,Y)                                            QH5  110
C           PARAMETER FCT REQUIRES AN EXTERNAL STATEMENT                QH5  120
C                                                                       QH5  130
C        DESCRIPTION OF PARAMETERS                                      QH5  140
C           FCT    - THE NAME OF AN EXTERNAL FUNCTION SUBPROGRAM USED.  QH5  150
C           Y      - THE RESULTING INTEGRAL VALUE.                      QH5  160
C                                                                       QH5  170
C        REMARKS                                                        QH5  180
C           NONE                                                        QH5  190
C                                                                       QH5  200
C        SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED                  QH5  210
C           THE EXTERNAL FUNCTION SUBPROGRAM FCT(X) MUST BE FURNISHED   QH5  220
C           BY THE USER.                                                QH5  230
C                                                                       QH5  240
C        METHOD                                                         QH5  250
C           EVALUATION IS DONE BY MEANS OF 5-POINT GAUSSIAN-HERMITE     QH5  260
C           QUADRATURE FORMULA, WHICH INTEGRATES EXACTLY WHENEVER       QH5  270
C           FCT(X) IS A POLYNOMIAL UP TO DEGREE 9.                      QH5  280
C           FOR REFERENCE, SEE                                          QH5  290
C           V.I.KRYLOV, APPROXIMATE CALCULATION OF INTEGRALS,           QH5  300
C           MACMILLAN, NEW YORK/LONDON, 1962, PP.129-130 AND 343-346.   QH5  310
C                                                                       QH5  320
C     ..................................................................QH5  330
C                                                                       QH5  340
      SUBROUTINE QH5(FCT,Y)                                             QH5  350
C                                                                       QH5  360
C                                                                       QH5  370
      X=2.020183                                                        QH5  380
      Z=-X                                                              QH5  390
      Y=.01995324*(FCT(X)+FCT(Z))                                       QH5  400
      X=.9585725                                                        QH5  410
      Z=-X                                                              QH5  420
      Y=Y+.3936193*(FCT(X)+FCT(Z))                                      QH5  430
      X=0.                                                              QH5  440
      Y=Y+.9453087*FCT(X)                                               QH5  450
      RETURN                                                            QH5  460
      END                                                               QH5  470