HP Prime (emulator): Zeta function query



#9

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.


#10

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".


#11

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.


#12

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.

#13

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


#14

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.

#15

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


#16

You are correct, the proof was for irrationality only.

- Pauli


Possibly Related Threads...
Thread Author Replies Views Last Post
  Prime Emulator Connection Problem John Colvin 3 525 12-14-2013, 11:00 PM
Last Post: Han
  HP Prime Emulator download John Colvin 2 406 12-14-2013, 05:54 PM
Last Post: John Colvin
  HP50g: Writing a function that returns a function Chris de Castro 2 448 12-10-2013, 06:49 PM
Last Post: Han
  IFERR function on HP Prime Mic 2 382 12-02-2013, 01:33 AM
Last Post: cyrille de Brébisson
  HP Prime: Dirichlet's eta function recognized but not numerically evaluated Helge Gabert 0 220 11-16-2013, 03:41 PM
Last Post: Helge Gabert
  Possible bug with sqrt function in the HP prime Michael de Estrada 6 489 11-15-2013, 12:49 PM
Last Post: Michael de Estrada
  HP-41 MCODE: The Last Function - at last! Ángel Martin 0 226 11-08-2013, 05:11 AM
Last Post: Ángel Martin
  [Prime] any ideas for a undo function? Stefan Dröge (Germany) 8 605 11-04-2013, 04:37 PM
Last Post: Damien
  HP Prime 'where' function bluesun08 11 793 10-29-2013, 06:56 PM
Last Post: Joe Horn
  HP Prime - Defining a function bluesun08 5 465 10-23-2013, 02:43 PM
Last Post: Han

Forum Jump: