Google
 

Trailing-Edge - PDP-10 Archives - decuslib20-02 - decus/20-0026/ql9.ssp
There are 2 other files named ql9.ssp in the archive. Click here to see a list.
C                                                                       QL9   10
C     ..................................................................QL9   20
C                                                                       QL9   30
C        SUBROUTINE QL9                                                 QL9   40
C                                                                       QL9   50
C        PURPOSE                                                        QL9   60
C           TO COMPUTE INTEGRAL(EXP(-X)*FCT(X), SUMMED OVER X FROM 0    QL9   70
C                               TO INFINITY).                           QL9   80
C                                                                       QL9   90
C        USAGE                                                          QL9  100
C           CALL QL9 (FCT,Y)                                            QL9  110
C           PARAMETER FCT REQUIRES AN EXTERNAL STATEMENT                QL9  120
C                                                                       QL9  130
C        DESCRIPTION OF PARAMETERS                                      QL9  140
C           FCT    - THE NAME OF AN EXTERNAL FUNCTION SUBPROGRAM USED.  QL9  150
C           Y      - THE RESULTING INTEGRAL VALUE.                      QL9  160
C                                                                       QL9  170
C        REMARKS                                                        QL9  180
C           NONE                                                        QL9  190
C                                                                       QL9  200
C        SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED                  QL9  210
C           THE EXTERNAL FUNCTION SUBPROGRAM FCT(X) MUST BE FURNISHED   QL9  220
C           BY THE USER.                                                QL9  230
C                                                                       QL9  240
C        METHOD                                                         QL9  250
C           EVALUATION IS DONE BY MEANS OF 9-POINT GAUSSIAN-LAGUERRE    QL9  260
C           QUADRATURE FORMULA, WHICH INTEGRATES EXACTLY WHENEVER       QL9  270
C           FCT(X) IS A POLYNOMIAL UP TO DEGREE 17.                     QL9  280
C           FOR REFERENCE, SEE                                          QL9  290
C           V.I.KRYLOV, APPROXIMATE CALCULATION OF INTEGRALS,           QL9  300
C           MACMILLAN, NEW YORK/LONDON, 1962, PP.130-132 AND 347-352.   QL9  310
C                                                                       QL9  320
C     ..................................................................QL9  330
C                                                                       QL9  340
      SUBROUTINE QL9(FCT,Y)                                             QL9  350
C                                                                       QL9  360
C                                                                       QL9  370
      X=26.37407                                                        QL9  380
      Y=.3290874E-10*FCT(X)                                             QL9  390
      X=18.83360                                                        QL9  400
      Y=Y+.4110769E-7*FCT(X)                                            QL9  410
      X=13.46624                                                        QL9  420
      Y=Y+.6592123E-5*FCT(X)                                            QL9  430
      X=9.372985                                                        QL9  440
      Y=Y+.3052498E-3*FCT(X)                                            QL9  450
      X=6.204957                                                        QL9  460
      Y=Y+.005599627*FCT(X)                                             QL9  470
      X=3.783474                                                        QL9  480
      Y=Y+.04746056*FCT(X)                                              QL9  490
      X=2.005135                                                        QL9  500
      Y=Y+.1992875*FCT(X)                                               QL9  510
      X=.8072200                                                        QL9  520
      Y=Y+.4112140*FCT(X)                                               QL9  530
      X=.1523222                                                        QL9  540
      Y=Y+.3361264*FCT(X)                                               QL9  550
      RETURN                                                            QL9  560
      END                                                               QL9  570