Occasionally Maple does not make me happy, if I want numerical integrals with more than the 14 digits which are supplied through the NAG libraries - that may be rather slow. So I wrote my own solution using LCC-WIN32 (a (free) compiler system allowing 104 decimal points of precision), where I coded the double exponential integration method. That gives me what I want to have in reasonable time. Details are sketched in the uploaded zip-file (it contains all what is needed to run the stuff). Some draw-back: this is for Windows only and because of OpenMaple at least Maple 9 is need. Feed back and further test results are appreciated, mostly I tested that approach for cumulative normal distributions in dim LT 3. The main differences against the posted Pari-method: using callbacks for the integrand allows to use all Maple functions (the external lib need not to have them, but that of course slows down the thing a bit) and precision is fixed (well, LCC uses its 104 decimal places for the involved data type). And I am (almost) sure one can use Watcom instead of MSVC (and I think the LCC solution is technical more stable as my somewhat hacked DLL for Pari). Just try, what's better for your purpose (if ever needed). Download 102_external-numerical-integration-with-many-digits-using-openmaple.zip
View file details

Please Wait...