(HP71B) Bernoulli Polynomials

[John Keith]

Quite an impressive program. Congratulations for representing Bernoulli numbers as rational numbers on a system that does not support rationals. However, ...

(07-17-2023 07:16 PM)C.Ret Wrote:  If someone wants to use this program on an HP-71B without this ROM, they will have to adapt by defining a user function:
5 DEF FNC0(N,K)=FACT(N)/FACT(K)/FACT(N-K)
and replacing every use of COMB(n,k) with FNC0(n,k) in the code.

Computing binomial coefficients this way is numerically unstable and will likely undo the advantages of maintaining Bernoulli numbers as integer ratios. This section of the Wikipedia page on binomial coefficients gives several useful methods.

Finally, not to nit-pick, but it's Ada, not Eva.