HP 50g Romberg Integration

[peacecalc]

Hello Gerson,

first of all, I tried to include with the tags my program, but it doesn't works, the format was spoiled.

Second the example:

The function to integrate is defined in 'Y' (I forgot this to mention in my first post, sorry) (here: << INV >> for 1/x)

Input:
4.
2.
1.

4 is iteration depth
2 upper limit
1 lower limit

output:

6,9412E-1
{ 6,9444E-1 6,9325E-1 6,9315E-1 }
{ 6,9317E-1 6,9315E-1 }
{ 6,9315E-1 }
6,9315E-1

First value: result from pure trapezoid rule in total
Last value: with Richardson extrapolation
the three lists in between are the "table" for the Richardson extrapolation
(this outputs results every DOSUBS call).