|
PROGRAM
|
(PARI) s(m)= { if (m==1, return(a1)); if (m==2, return(a2)); r = (3*(2*m - 3)*a2 - (m - 3)*a1)/m; a1=a2; a2=r; return(r); } { a1=1; a2=1; n=0; for (m=1, 300, a=s(m); if (isprime(a), write("b092839.txt", n++, " ", m))); } [From Harry J. Smith (hjsmithh(AT)sbcglobal.net), Jun 21 2009]
|