Google
 

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