d01fcf

d01fcf © Numerical Algorithms Group, 2002.

Purpose

D01FCF Multi-dimensional adaptive quadrature over hyper-rectangle

Synopsis

[acc,finval,minpts,ifail] = d01fcf(a,b,functn,epslon<,minpts,maxpts,ifail>)

Description

 
 The routine returns an estimate of a multi-dimensional integral 
 over a hyper-rectangle (i.e., with constant limits), and also an 
 estimate of the relative error. The user sets the relative 
 accuracy required, supplies the integrand as a function 
 subprogram (FUNCTN), and also sets the minimum and maximum 
 acceptable number of calls to FUNCTN (in MINPTS and MAXPTS).
 
 The routine will fail if the requested relative error level has 
 not been attained by the time MAXPTS calls to FUNCTN have been 
 made; or, if the amount LENWRK of working storage is 
 insufficient. If a smaller value is used, and is exhausted in 
 the course of execution, the routine switches to a less 
 efficient mode of operation; only if this mode also breaks down 
 is insufficient storage reported.
 

Parameters

d01fcf

Required Input Arguments:

a (:)                                 real
b (:)                                 real
functn                                function (User-Supplied)
epslon                                real

Optional Input Arguments:                       <Default>

minpts                                integer  0
maxpts                                integer  30000
ifail                                 integer  -1

Output Arguments:

acc                                   real
finval                                real
minpts                                integer
ifail                                 integer