(HP15C)(HP67)(HP41C) Bernoulli Polynomials

[John Keith]

Now that we have (mostly) conquered Bernoulli numbers, we can get back to the original purpose of this thread which is Bernoulli polynomials. Albert's second program in post #20 is simple and clearly written, and should be easy to implement on any HP calculator that can handle Bernoulli numbers.

To start off, here is an implementation that will run on any RPL calculator. In this example, BERN can refer to the first program in this post in the other thread, or the program here. In the latter case, BERN should be followed by / (divide).

Code:


\<< \-> m x
  \<< 1. 1. 1. m
    FOR k m k - 1. + * k /
      SWAP x * OVER k BERN * + SWAP
    NEXT DROP
  \>>
\>>