[R] Integration with "adapt"

[Martin Maechler]

>>>>> "Rodrigo" == Rodrigo Drummond <rduarte at unicamp.br>
>>>>>     on Tue, 27 Jul 2004 15:31:07 -0300 (BRT) writes:

    Rodrigo> Hi all, I need to calculate a multidimensional
    Rodrigo> integration on R. I am using the command
    Rodrigo> "adapt" (from library adapt), although

it's a "package", not a library.

    Rodrigo> sometimes I get the following error message:

    Rodrigo> Ifail=2, lenwrk was too small. -- fix adapt() !
    Rodrigo> Check the returned relerr! in: adapt(3, linf, lsup,
    Rodrigo> functn = Integrando1)

If you could give as a *reproducible* example,

we (the adapt authors) would have chance to do what the above
message says, namely "fix adapt()" ..


    Rodrigo> I guess it happens because the domain of
    Rodrigo> integration is too small,

maybe, maybe not. We need an example we can reproduce, see above.

    Rodrigo> although I tried a change of variables to avoid
    Rodrigo> this problem and it didn’t help. The command adapt
    Rodrigo> calls a fortran routine, but I don’t know fortran
    Rodrigo> enough to fix the problem.

Martin Maechler, ETH Zurich