Return to contfrac.gp CVS log | Up to [local] / OpenXM_contrib / pari-2.2 / examples |
File: [local] / OpenXM_contrib / pari-2.2 / examples / Attic / contfrac.gp (download)
Revision 1.1.1.1 (vendor branch), Tue Oct 2 11:16:59 2001 UTC (22 years, 10 months ago) by noro
Imported pari-2.2.1(alpha). |
period(D) = { local(u,v,j,r,s); if (type(D) != "t_INT" || D < 2, return(-1)); u = sqrtint(D); v = D-u^2; if (!v, return(0)); s = v; r = u; j = 0; until (u==r && v==s, u = (r+u)\v * v - u; v = (D-u^2)\v; j++; ); j }