The most common command for numerical integration is evalf(Int(f, x=a..b)) where the integration command is expressed in inert form to avoid first invoking the symbolic integration routines. It is also possible to invoke evalf on an unevaluated integral returned by the symbolic int command, as in evalf(int(f, x=a..b)), if it happens that symbolic int fails (returns an unevaluated integral).

All numerical integration calling sequences can also be accessed directly from the int command by using the numeric option.

You can enter the command evalf/Int using either the 1-D or 2-D calling sequence. For example, evalf(Int(1/(x^2+1), x=0..infinity)) is equivalent to $\mathrm{evalf}\left({\int }_0^{\mathrm{\infty }}\frac{1}{x^2+1}\ⅆx\right)$.

The integrand f may be another unevaluated integral, that is, multiple integrals are supported. A special list syntax (see below) can be used to specify multiple integrals, rather than using nested integrals. Integrals expressed in the standard non-list notation are referred to as 1-D (one-dimensional) integrals including the case of nested 1-D integrals.

If the integrand f is specified as a procedure or a Maple operator, then the second argument must be a range a..b and not an equation, that is, a variable of integration must not be specified.

Various levels of user information are displayed during the computation if infolevel[`evalf/int`] is assigned values between 1 and 4.

Additional options may be specified as equations. (For backward compatibility some options are accepted as values rather than equations, as specified below.) An option is one of the following forms:

The specification method = <name> (or simply <name>) indicates a particular numerical integration method to be applied. The methods that can be specified are described below. By default, a hybrid symbolic-numeric strategy is applied.

The specification digits = <posint> (or simply <posint>) indicates the number of digits of precision for the computation. Some additional guard digits are carried during the computation to attempt to achieve a result with <posint> correct digits (although a larger tolerance can be specified by using the 'epsilon' option). By default, the Maple environment variable Digits specifies the precision for the computation.

The specification epsilon = <numeric> specifies the relative error tolerance for the computed result. The routines attempt to achieve a final result with a relative error less than this value. By default, the relative error tolerance which the routines attempt to achieve for the final result is

where digits is the precision specified for the computation. In attempting to achieve this accuracy, the working value of Digits is increased as deemed necessary. It is an error to specify 'epsilon' smaller than the default value above, and for any value larger than 1e-3 the value 1e-3 is used instead if the method in use is deterministic (i.e. not the MonteCarlo or Cuba methods).

Note: For some integrands, the numerical accuracy attained when computing values of the integrand may be insufficient to allow the value of the integral to be computed to the default tolerance $\mathrm{eps}$ (even though the computation is using some number of guard digits). In such cases, specifying a larger tolerance (relative to the setting of digits) via the 'epsilon' option may be helpful. Alternatively, increasing Digits and fixing 'epsilon' may provide the desired answer (see the end of the examples section).

The specification methodoptions = <list> specifies a list of zero or more options that are specific to a method selected with the method option. In particular, if the _d01ajc or _d01akc method is selected, one can supply an option of the form methodoptions=[maxintervals = <posint>] to specify a maximal number of subintervals that can be used internally by those methods.

For backward compatibility, the option maxintervals = <posint> for the _d01ajc and _d01akc methods can also be specified as a separate option, as an argument to Int directly, rather than in the methodoptions option.

In the default case (no particular method specified), the problem is first passed to NAG integration routines if Digits is not too large (that is, if Digits <= evalhf(Digits)). The NAG routines are in a compiled C library and hence operate at hardware floating-point speed. If the NAG routines cannot perform the integration, then some singularity handling may be performed and control may pass back to the NAG routines with a modified problem. Native Maple routines are invoked if the NAG routines cannot solve the problem (for example, if Digits is too large or if the integrand involves functions for which hardware floating-point evaluation is not supported).

The native Maple hybrid symbolic-numeric solution strategy is as follows. The default numerical method applied is Clenshaw-Curtis quadrature (_CCquad). If slow convergence is detected then there must be singularities in or near the interval of integration (perhaps in the complex plane). Some techniques of symbolic analysis are used to deal with the singularities. For problems with non-removable endpoint singularities, an adaptive double-exponential quadrature method (_Dexp) is applied.

If singularities interior to the interval are suspected, then an attempt is made to locate the singularities in order to break up the interval of integration. Finally, if still unsuccessful, then the interval is subdivided and the _Dexp method is applied, or if the method was already _Dexp or _Sinc then an adaptive Gaussian quadrature method (_Gquad) is applied.

For the limits of integration, the values infinity and/or -infinity are valid, and a symbolic-numeric strategy attempts to deal with singularities. Techniques employed include variable transformations, subtracting out the singularity, and integration of a truncated generalized series near the singularity.

No singularity handling is attempted in the case where the integrand f is specified as a procedure or a Maple operator.

A numerical multiple integration problem may be specified in a natural way using nested one-dimensional integrals, for example:

where the integrand f depends on x1, x2, ..., xn. Such a problem may also be specified using the following special multiple integration notation with a list as the second argument:

Additional optional arguments may be stated just as in the case of 1-D integration. Also as in 1-D integration, the integrand f may be specified as a procedure in which case the second argument must be a list of ranges: [a1..b1, a2..b2, ..., an..bn].

Whether a multiple integration problem is stated using nested integrals or using the list notation, the arguments will be extracted so as to invoke the same numerical multiple integration routines.

The optional argument method = <name> (or simply <name>) accepts the following method names.

$\mathrm{evalf}\left(\mathrm{Int}\left(\frac{\exp \left(-x^3\right)}{x^2+1}\&comma;x=0..1\right)\right)$

$0.6649369431$

$\mathrm{evalf}\left(\mathrm{Int}\left(\frac{1}{x^2+1}\&comma;x=0..\mathrm{\infty }\right)\right)$

$1.570796327$

$\mathrm{evalf}\left(\mathrm{Int}\left(\sin \left(x\right)\ln \left(x\right)\exp \left(-x^3\right)\&comma;x=0..\mathrm{\infty }\right)\right)$

$\mathrm{-0.1957885158}$

$\mathrm{e1}≔\frac{1}{\mathrm{\Gamma }\left(x\right)}\&colon;$

$\mathrm{Int}\left(\mathrm{e1}\&comma;x=0..2\right)=\mathrm{evalf}\left(\mathrm{Int}\left(\mathrm{e1}\&comma;x=0..2\&comma;\mathrm{digits}=20\&comma;\mathrm{method}=\mathrm{\_Dexp}\right)\right)$

${\int }_0^2\frac{1}{\mathrm{\Gamma }\left(x\right)}\&DifferentialD;x=1.6263783986861406145$

$\mathrm{e2}≔\frac{\exp \left(v-\frac{v^2}{2}\right)}{1+\frac{1}{2}\exp \left(v\right)}\&colon;$

$\mathrm{Int}\left(\mathrm{e2}\&comma;v=0..\mathrm{\infty }\right)=\mathrm{evalf}\left[20\right]\left(\mathrm{Int}\left(\mathrm{e2}\&comma;v=0..\mathrm{\infty }\right)\right)$

${\int }_0^{\mathrm{\infty }}\frac{{\&ExponentialE;}^{v-\frac{1}{2}{v}^2}}{1+\frac{{\&ExponentialE;}^v}{2}}\&DifferentialD;v=1.3055168991185060654$

$\mathrm{e3}≔\frac{1}{1+\ln \left(1+x\right)}\&colon;$

$\mathrm{Int}\left(\mathrm{e3}\&comma;x=0..1\right)=\mathrm{evalf}\left[32\right]\left(\mathrm{Int}\left(\mathrm{e3}\&comma;x=0..1\&comma;\mathrm{method}=\mathrm{\_Gquad}\right)\right)$

${\int }_0^1\frac{1}{1+\ln \left(1+x\right)}\&DifferentialD;x=0.73716070962368003213791626905536$

$r≔\mathrm{int}\left(\mathrm{sech}\left(x\right)\exp \left(-x^2\right)\&comma;x=-\mathrm{\infty }..\mathrm{\infty }\right)$

$r≔{\int }_{-\mathrm{\infty }}^{\mathrm{\infty }}\mathrm{sech}\left(x\right){\&ExponentialE;}^{-x^2}\&DifferentialD;x$

$\mathrm{evalf}\left(r\right)$

$1.479061171$

$\mathrm{evalf}\left[25\right]\left(r\right)$

$1.479061171449575890854454$

f := proc(x) if x < 2 then 2*x else x^2 end if; end proc;

$f≔\mathbf{proc}\left(x\right)\mathbf{if}x<2\mathbf{then}2\&ast;x\mathbf{else}x\&Hat;2\mathbf{end\; if}\mathbf{end\; proc}$

$\mathrm{evalf}\left(\mathrm{Int}\left(f\left(x\right)\&comma;x=0..3\right)\right)$

Error, (in f) cannot determine if this expression is true or false: x < 2

$\mathrm{evalf}\left(\mathrm{Int}\left(f\&comma;0..3\right)\right)$

$10.33333333$

$\mathrm{evalf}\left(\mathrm{Int}\left('f'\left(x\right)\&comma;x=0..3\right)\right)$

$10.33333333$

$\mathrm{Int}\left(\mathrm{Int}\left(\mathrm{Int}\left(\mathrm{`/`}\left(\exp \left(\mathrm{`+`}\left(\mathrm{`+`}\left(x\&comma;y\right)\&comma;z\right)\right)\&comma;\mathrm{`*`}\left(\mathrm{`*`}\left(\mathrm{`+`}\left(\mathrm{`*`}\left(5\&comma;x\right)\&comma;1\right)\&comma;\mathrm{`+`}\left(\mathrm{`*`}\left(10\&comma;y\right)\&comma;2\right)\right)\&comma;\mathrm{`+`}\left(\mathrm{`*`}\left(15\&comma;z\right)\&comma;3\right)\right)\right)\&comma;\mathrm{`=`}\left(x\&comma;0..4\right)\right)\&comma;\mathrm{`=`}\left(y\&comma;0..3\right)\right)\&comma;\mathrm{`=`}\left(z\&comma;0..\mathrm{sqrt}\left(2\right)\right)\right)$

${\int }_0^{\sqrt{2}}\left({\int }_0^3\left({\int }_0^4\frac{{\&ExponentialE;}^{x+y+z}}{\left(5x+1\right)\left(10y+2\right)\left(15z+3\right)}\&DifferentialD;x\right)\&DifferentialD;y\right)\&DifferentialD;z$

$\mathrm{evalf}\left(\right)$

$0.9331611325$

$d≔1-w^2{x}^2{y}^2{z}^2\&colon;$

$g≔d\cos \left(wxyz\right)-dwxyz\sin \left(wxyz\right)$

$g≔\left(-w^2{x}^2{y}^2{z}^2+1\right)\cos \left(wxyz\right)-\left(-w^2{x}^2{y}^2{z}^2+1\right)wxyz\sin \left(wxyz\right)$

$\mathrm{evalf}\left(\mathrm{Int}\left(g\&comma;\left[w=0..1\&comma;x=0..1\&comma;y=0..1\&comma;z=0..1\right]\right)\right)$

$0.9717798177$

$h≔\frac{1}{2+\sin \left(\mathrm{\pi }\mathrm{sqrt}\left(87\right)\left(\mathrm{x1}+\mathrm{x2}+\mathrm{x3}+\mathrm{x4}+\mathrm{x5}+\mathrm{x6}\right)\right)}$

$h≔\frac{1}{2+\sin \left(\mathrm{\pi }\sqrt{87}\left(\mathrm{x1}+\mathrm{x2}+\mathrm{x3}+\mathrm{x4}+\mathrm{x5}+\mathrm{x6}\right)\right)}$

$\mathrm{evalf}\left(\mathrm{Int}\left(h\&comma;\left[\mathrm{x1}=-1..1\&comma;\mathrm{x2}=-1..1\&comma;\mathrm{x3}=-1..1\&comma;\mathrm{x4}=-1..1\&comma;\mathrm{x5}=-1..1\&comma;\mathrm{x6}=-1..1\right]\&comma;\mathrm{method}=\mathrm{\_MonteCarlo}\&comma;\mathrm{\epsilon }=0.005\right)\right)$

$36.91495206$

$\mathrm{evalf}\left[3\right]\left(\right)$

$36.9$

$\mathrm{igrand}≔\frac{1}{2-\sin \left(\mathrm{\pi }x\right)-\sin \left(\frac{355}{113}x\right)}$

$\mathrm{igrand}≔\frac{1}{2-\sin \left(\mathrm{\pi }x\right)-\sin \left(\frac{355x}{113}\right)}$

$\mathrm{evalf}\left(\mathrm{eval}\left(\mathrm{igrand}\&comma;x=0.5\right)\right)$

$Float\left(\mathrm{\infty }\right)$

$\mathrm{evalf}\left(\mathrm{Int}\left(\mathrm{igrand}\&comma;x=0..1\right)\right)$

${\int }_{0.}^{1.}\frac{1}{2.-1.\sin \left(3.141592654x\right)-1.\sin \left(3.141592920x\right)}\&DifferentialD;x$

$\mathrm{evalf}\left[15\right]\left(\mathrm{eval}\left(\mathrm{igrand}\&comma;x=0.5\right)\right)$

$1.11111111111111\times {10}^{14}$

$\mathrm{evalf}\left[20\right]\left(\mathrm{eval}\left(\mathrm{igrand}\&comma;x=0.5\right)\right)$

$1.1241778044582643369\times {10}^{14}$

$\mathrm{evalf}\left[25\right]\left(\mathrm{eval}\left(\mathrm{igrand}\&comma;x=0.5\right)\right)$

$1.124177605944406775029070\times {10}^{14}$

$\mathrm{evalf}\left[25\right]\left(\mathrm{Int}\left(\mathrm{igrand}\&comma;x=0..1\&comma;\mathrm{\epsilon }=1.\times {10}^{\mathrm{-10}}\right)\right)$

$1.499451605234141071490295\times {10}^7$

$\mathrm{igrand}≔\frac{\sin \left(\exp \left(\mathrm{abs}\left(x\right)\right)\right)}{1+x^2}$

$\mathrm{igrand}≔\frac{\sin \left({\&ExponentialE;}^{\left|x\right|}\right)}{x^2+1}$

$\mathrm{evalf}\left(\mathrm{Int}\left(\mathrm{igrand}\&comma;x=-10..10\&comma;\mathrm{method}=\mathrm{\_d01akc}\right)\right)$

Error, (in `evalf/int`) NE_QUAD_MAX_SUBDIV:
  The maximum number of subdivisions has
  been reached: max_num_subint = 500

$\mathrm{evalf}\left(\mathrm{Int}\left(\mathrm{igrand}\&comma;x=-10..10\&comma;\mathrm{method}=\mathrm{\_d01akc}\&comma;\mathrm{methodoptions}=\left[\mathrm{maxintervals}=1000\right]\right)\right)$

$1.235076653$