<< intprb_crudemc Integration Algorithms intprb_getcrudemc >>

Intprbs >> Intprbs > Integration Algorithms > intprb_getcrudeld

intprb_getcrudeld

Uses a crude Quasi Monte-Carlo to estimate the integral.

Calling Sequence

integr = intprb_getcrudeld ( nprob , n , p , callf )
integr = intprb_getcrudeld ( nprob , n , p , callf , ldseq )

Parameters

nprob:

the problem number

n:

a 1-by-1 matrix of floating point integers, the spatial dimension.

p:

a floating point matrix, the parameters of the problem

callf :

a 1-by-1 matrix of floating point integers, the number of calls to the function (default = 1.e4)

ldseq :

a 1-by-1 matrix of strings, the name of the sequence (Default = "sobolf"). The list of available sequences can be computed by lowdisc_methods().

integr :

a 1-by-1 matrix of doubles, the approximate value of the integral, the mean of the function.

accur :

a 1-by-1 matrix of doubles, the estimated error on the integral.

varf :

a 1-by-1 matrix of doubles, the approximate value of the variance of f.

Description

Uses a crude Quasi Monte-Carlo to approximate the integral corresponding to the nprob problem. The computation is based on a Low Discrepancy sequence. This function is a wrapper on intprb_crudeld.

Examples

// Consider problem #1
nprob = 1;
[n,p] = intprb_getsetup(nprob);
callf = 1.e4;
// Get the exact result
e = intprb_getexpect(n,p,nprob)
// Use a Low Discrepancy sequence
integr = intprb_getcrudeld ( nprob , n , p , callf )

Authors


Report an issue
<< intprb_crudemc Integration Algorithms intprb_getcrudemc >>