D106: Gaussian Quadrature with Five- and Six-Point Rules

Author(s): F. James Library: MATHLIB
Submitter: Submitted: 01.12.1994
Language: Fortran Revised:

Subroutine subprograms RGS56P and DGS56P calculate an approximation and uncertainty for the integral

displaymath70

equal respectively to the mean value and the difference of the results tex2html_wrap_inline72 and tex2html_wrap_inline74 obtained by the five- and six-point Gaussian integration rules.

On CDC and Cray computers, the double-precision version DGS56P is not available.

Structure:

SUBROUTINE subprograms
User Entry Names: RGS56P, DGS56P
External References: User-supplied FUNCTION subprogram.

Usage:

For tex2html_wrap_inline76 (type REAL), tex2html_wrap_inline78 (type DOUBLE PRECISION), 

    CALL tGS56P(F,A,B,RES,ERR)
F
(type according to t) Name of a user-supplied FUNCTION subprogram, declared EXTERNAL in the calling program. This subprogram must set tex2html_wrap_inline80 .
A,B
(type according to t) End-points of integration interval. Note that B may be less than A.
RES
(type according to t) The calculated approximation for I, i.e. tex2html_wrap_inline84 ,
ERR
(type according to t) An estimated uncertainty on this approximation, i.e. tex2html_wrap_inline86 .
tex2html_wrap_inline88

Michel Goossens Tue Jun 4 23:42:36 METDST 1996