HP 50g Romberg Integration

[peacecalc]

Hello Gershon,

Quote:but I've read the legacy algorithm is not the best one around (it's only there because it's been tried and tested since the HP-34C, the Saturn code since the HP-71B, I think).

That maybe is true, but it is not fair, because when you use the Romberg method with the iteration depth of 20, the program have to calculate 2^20 - 1 function values, this are at least over a million, no wonder about the time you have to wait...
I don't know wether a side effect occures: so many additions can cause a summation error???

Did you try the lower iteration depths and how the results behavier is?

And your test function is a hard piece of cake because of its limits at 0 and 1, both are divergent.

Greetings
peacecalc