OFFSET
0,1
LINKS
PROG
(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, ", ")); ) }
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Cino Hilliard, Aug 05 2003
EXTENSIONS
The next term is too large to include.
STATUS
approved