HP50G Bisection methode

[acser]

Thanks, Albert. I have not yet looked at the Ridders' method.

Here's the updated HP50g RPL code:

Copy and paste the below code under the "---" to the Emu48 stack then use the HP 50g program ASCIBIN.49 found at https://www.hpcalc.org/details/3648 to convert it to a proper HP 50g program on the stack.

‰ is the <= sign ([LS] [X])
and
Š is the HP >= sign ([LS] [1/x])
---

%%HP: T(1)A(R)F(.);
«
IF OVER F OVER F > THEN SWAP END
0. 0. -> R A B C FC
«
DO A B + 2. / 'C'
STO 'F(C)' EVAL
'FC' STO
IF FC 0. ‰
THEN C 'A' STO
END
IF FC 0. Š
THEN C 'B' STO
END
UNTIL 'ABS(A-B)<
.0000001'
END C
»
»