Routine ID: D101 | |
---|---|
Author(s): K.S. Kölbig | Library: MATHLIB |
Submitter: | Submitted: 15.01.1988 |
Language: Fortran | Revised: 15.03.1993 |
Function subprograms SIMPS and DSIMPS use Simpson's rule to compute an approximate value of the integral
On CDC or Cray computers, the double-precision version DSIMPS is not available.
Structure:
FUNCTION subprograms
User Entry Names: SIMPS, DSIMPS
Files Referenced: Unit 6
External References: MTLMTR (N002), ABEND (Z035)
Usage:
In any arithmetic expression, SIMPS(F,A,B,N) or DSIMPS(F,A,B,N)
has the approximate value of the integral I, where SIMPS is of type REAL and DSIMPS is of type DOUBLE PRECISION, and F, A, B have the same type as the function name. N is of type INTEGER.
equally-spaced points , with and .
Error handling:
Error D101.1:
or
odd.
The function value is set equal to zero, and a message is written on
Unit 6, unless subroutine MTLSET (N002) has been called.