The Museum of HP Calculators

HP Forum Archive 21

[ Return to Index | Top of Index ]

HP Prime (emulator): Zeta function query
Message #1 Posted by Paul Townsend (UK) on 29 Sept 2013, 1:39 p.m.

The HP Prime emulator seems to have difficulty with Zeta() of an odd positive integer, retaining the function call unevaluated even after Eval(). In a plot view such values show as "Undefined". What's going on? Zeta() is well-defined at these values and my WP34s (the real thing) has no difficulty evaluating Zeta(3) for example to be 1.20205...

Zeta() of an even integer is correctly handled returning the appropriate expression involving pi, and Zeta() of non-integral values is also correctly handled.

      
Re: HP Prime (emulator): Zeta function query
Message #2 Posted by Helge Gabert on 29 Sept 2013, 2:49 p.m.,
in response to message #1 by Paul Townsend (UK)

In CAS mode, hit approx key, after the Prime returns Zeta(3), then you also get the desired numerical result (1.202...).

I agree, in plot view you get "undefined".

            
Re: HP Prime (emulator): Zeta function query
Message #3 Posted by Paul Townsend (UK) on 29 Sept 2013, 3:41 p.m.,
in response to message #2 by Helge Gabert

I wonder if it is evaluating Zeta(n) recursively using Zeta(n-2)? Starting with Zetas(3) would result in the algorithm boggling on Zeta(1) and bombing out with an error.

Asking for Zeta(1) correctly returns infinity.

                  
Re: HP Prime (emulator): Zeta function query
Message #4 Posted by Helge Gabert on 29 Sept 2013, 3:51 p.m.,
in response to message #3 by Paul Townsend (UK)

Yes, I bet that's it! The algorithm should change depending on the user input - - but since this "only" affects odd integers > 1 it was probably either overlooked or deemed to be too much trouble to implement.

      
Re: HP Prime (emulator): Zeta function query
Message #5 Posted by Paul Dale on 29 Sept 2013, 6:34 p.m.,
in response to message #1 by Paul Townsend (UK)

Sounds like this is related to zeta(2n) being expressible relatively simply in terms of PI and the Bernoulli numbers whereas zeta(2n+1) are not.

Zeta(3) is Apéry's constant which has been proven to be transcendental, there are infinite series expressions for it but nothing really simple and certainly nothing simpler than just zeta(3).

- Pauli

            
Re: HP Prime (emulator): Zeta function query
Message #6 Posted by Gerson W. Barbosa on 29 Sept 2013, 6:50 p.m.,
in response to message #5 by Paul Dale

That's what I thought, but that works for non-integer values and this is defined for integers only.

http://people.math.sfu.ca/~cbm/aands/page_807.htm

(Abramowitz and Stegun: Handbook of Mathematical Functions, page 807, formula 23.2.16)

Gerson.

P.S.: This HP 50g program has the same behaviour regarding integer arguments.

%%HP: T(3)A(D)F(.);
\<< -22 CF \-> n
  \<< '(2.*\pi)^n/(2.*n!)' \->NUM n IBERNOULLI ABS INV /
  \>>
\>>

Edited: 29 Sept 2013, 6:54 p.m.

            
Re: HP Prime (emulator): Zeta function query
Message #7 Posted by mjcohen on 29 Sept 2013, 8:05 p.m.,
in response to message #5 by Paul Dale

I don't think it has been proven transcendental - only irrational.

                  
Re: HP Prime (emulator): Zeta function query
Message #8 Posted by Paul Dale on 29 Sept 2013, 8:17 p.m.,
in response to message #7 by mjcohen

You are correct, the proof was for irrationality only.

- Pauli


[ Return to Index | Top of Index ]

Go back to the main exhibit hall