home *** CD-ROM | disk | FTP | other *** search
/ Power-Programmierung / CD1.mdf / fortran / library / ssp / gauquad / qg4.for < prev    next >
Encoding:
Text File  |  1985-11-29  |  1.4 KB  |  47 lines
  1. C
  2. C     ..................................................................
  3. C
  4. C        SUBROUTINE QG4
  5. C
  6. C        PURPOSE
  7. C           TO COMPUTE INTEGRAL(FCT(X), SUMMED OVER X FROM XL TO XU)
  8. C
  9. C        USAGE
  10. C           CALL QG4 (XL,XU,FCT,Y)
  11. C           PARAMETER FCT REQUIRES AN EXTERNAL STATEMENT
  12. C
  13. C        DESCRIPTION OF PARAMETERS
  14. C           XL     - THE LOWER BOUND OF THE INTERVAL.
  15. C           XU     - THE UPPER BOUND OF THE INTERVAL.
  16. C           FCT    - THE NAME OF AN EXTERNAL FUNCTION SUBPROGRAM USED.
  17. C           Y      - THE RESULTING INTEGRAL VALUE.
  18. C
  19. C        REMARKS
  20. C           NONE
  21. C
  22. C        SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED
  23. C           THE EXTERNAL FUNCTION SUBPROGRAM FCT(X) MUST BE FURNISHED
  24. C           BY THE USER.
  25. C
  26. C        METHOD
  27. C           EVALUATION IS DONE BY MEANS OF 4-POINT GAUSS QUADRATURE
  28. C           FORMULA, WHICH INTEGRATES POLYNOMIALS UP TO DEGREE 7
  29. C           EXACTLY.
  30. C           FOR REFERENCE, SEE
  31. C           V.I.KRYLOV, APPROXIMATE CALCULATION OF INTEGRALS,
  32. C           MACMILLAN, NEW YORK/LONDON, 1962, PP.100-111 AND 337-338.
  33. C
  34. C     ..................................................................
  35. C
  36.       SUBROUTINE QG4(XL,XU,FCT,Y)
  37. C
  38. C
  39.       A=.5*(XU+XL)
  40.       B=XU-XL
  41.       C=.4305682*B
  42.       Y=.1739274*(FCT(A+C)+FCT(A-C))
  43.       C=.1699905*B
  44.       Y=B*(Y+.3260726*(FCT(A+C)+FCT(A-C)))
  45.       RETURN
  46.       END
  47.