d01alf

d01alf © Numerical Algorithms Group, 2002.

Purpose

D01ALF 1-D quadrature, adaptive, finite interval, allowing for singularities at user-specified break-points

Synopsis

[result,abserr,w,iw,ifail] = d01alf(f,a,b,points,epsabs,epsrel<,lw,liw,ifail>)

Description

 
 D01ALF is very similar to D01AJF, but allows the user to supply
 'break-points', points at which the function is known to be 
 difficult. It is an adaptive routine, using the Gauss 10-point 
 and Kronrod 21-point rules. The user-supplied 'break-points'
 always occur as the end-points of some sub-interval during the 
 adaptive process.

Parameters

d01alf

Required Input Arguments:

f                                     function (User-Supplied)
a                                     real
b                                     real
points (:)                            real
epsabs                                real
epsrel                                real

Optional Input Arguments:                       <Default>

lw                                    integer  d01alf11
liw                                   integer  d01alf13(lw)
ifail                                 integer  -1

Output Arguments:

result                                real
abserr                                real
w (lw)                                real
iw (liw)                              integer
ifail                                 integer