login
Primes in the numerator of the continued fraction rational approximation of zeta(3).
0

%I #5 Oct 01 2013 17:57:42

%S 5,113,1987,552493,628313002458512784191921,

%T 40755082849497410605337341,6681921617166540622940410282864619819

%N Primes in the numerator of the continued fraction rational approximation of zeta(3).

%H Cino Hilliard, <a href="http://groups.msn.com/BC2LCC/continuedfractions.msnw">Continued fractions rational approximation of numeric constants</a>.

%o (PARI) \Continued fractions rational approximation of numeric functions cfrac(m,f) = x=f; for(n=0,m,i=floor(x); x=1/(x-i); print1(i,",")) cfraczeta(m,f) = { cf = vector(100000); x=f; for(n=0,m, i=floor(x); x=1/(x-i); cf[n+1] = i; ); for(m1=0,m, r=cf[m1+1]; forstep(n=m1,1,-1, r = 1/r; r+=cf[n]; ); numer=numerator(r); denom=denominator(r); if(isprime(numer),print1(numer,",")); ) }

%K easy,nonn

%O 0,1

%A _Cino Hilliard_, Aug 05 2003

%E The next term is too large to include.