
PROG

(PARI) \Continued fraction rational approximation of numeric constants f. m=steps. cfracnumprime(m, f) = { default(realprecision, 3000); cf = vector(m+10); x=f; for(n=0, m, i=floor(x); x=1/(xi); 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(ispseudoprime(numer), print1(numer, ", ")); ) }
