HP49-50G Willan's formulae to find p(n), the n-th prime number

[Gil]

I tried to put, from Wikipedia, Willan's formulae to find p(n), the n-th prime number, on the HP50G.

Code:


\<< \-> n '1+\GS(ii=1,2^n,FLOOR((n/\GS(j=1,ii,FLOOR(COS(((j-1)!+1)*\pi/j)^2)))^(1/n)))' \->NUM
\>>

For n=1, I get 2, ok
For n=2, I get 3, ok
For n=3, I get 5, ok
For n=4, I get 7, ok
For n=5, I get 33!
For n=6, I get 65!

Could somebody tell me what is wrong in my code above?

NB
However, the following user RPN code works nicely with the loop instructions FOR... NEXT.

Code:


\<< \-> n
  \<< 1 1 2 n ^
    FOR ii 0 1 ii
      FOR j j 1 - ! 1 + j / \pi * COS SQ FLOOR +
      NEXT n SWAP / n INV ^ FLOOR +
    NEXT
  \>>
\>>

Thanks for your help.

Gil