This is pretty fun. How many digits of accuracy do you want? I think that needs to be specified or else there is no
point to it.
Rich
--------------------------------------------------
From: "Robert Dodier" <robert.dodier at gmail.com>
Sent: Tuesday, March 16, 2010 12:38 AM
To: "Luigi Marino" <luigi_marino2 at alice.it>
Cc: <maxima at math.utexas.edu>
Subject: Re: [Maxima] numerical integration
> On 3/15/10, Luigi Marino <luigi_marino2 at alice.it> wrote:
>
>> Can Maxima integrate numerical
>> the function x^x^x from 1 to 6 or from 5 to 6 ?
>> It is a problem of the ISSAC'97 Systems Challenge.
>
> You can try the bromberg (bigfloat romberg) function.
> load(brmbrg) loads it. It is undocumented, sorry, but it
> works like romberg.
>
> Just playing around a little, I get a result for
> bromberg(x^x^x,x,5,5.01) => 4.834235552810633b2240
> (no idea how accurate that is) but it fails with
> "too many iterations" or something like that for
> bromberg(x^x^x,x,5,5.1). Maybe you can increase
> brombergit (maxima number of iterations).
>
> Or maybe there is another approach -- perhaps by
> combining symbolic and numerical results somehow.
> Just a guess.
>
> Have fun,
>
> Robert Dodier
> _______________________________________________
> Maxima mailing list
> Maxima at math.utexas.edu
> http://www.math.utexas.edu/mailman/listinfo/maxima
>