Riemann's Zeta Function - another approach (RPL)

[Dieter]

(07-02-2017 11:53 PM)Gerson W. Barbosa Wrote:  Great! It needs to be accurate only for x > 0.5. Then the functional equation can handle the rest if Gamma is available.

I have tried also this on the '41, using the FACT function so that at least integer arguments can be handled, and it worked fine. Now, if there only was a Gamma function...

BTW, does anyone know if there was a ROM that included Gamma in Mcode (!) ? Except the ones by Ángel, that is. ;-) The Advantage ROM was said to update the '41 to features that else were only available on the 34C and 15C (Solve, Integrate, Matrix operations), but AFAIK this did not include Gamma.

(07-02-2017 11:53 PM)Gerson W. Barbosa Wrote:  Your code can be nicely pasted into Free42, except line 91 (because the space between the last digit and E).

The listing was generated by hp41uc. But remember that all the math is designed for 10-digit precision. On the 42s you may be better off with a constant 5 E–9. Try it and see what you get. Maybe you want to round the results to 10 digits by SCI 9 RND ALL. ;-) On the 42s you can also apply RCL-arithmetics which will save a few steps.

(07-02-2017 11:53 PM)Gerson W. Barbosa Wrote:  LBL 2 is not necessary, unless there's a reason to keep it:

The reason is a common exit point at the end of the program. This way it can be restarted with a simple R/S. ;-)

Dieter