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