Return to contfrac.gp CVS log

Up to [local] / OpenXM_contrib / pari-2.2 / examples

Annotation of OpenXM_contrib/pari-2.2/examples/contfrac.gp, Revision 1.1.1.1

1.1       noro        1: period(D) =
                      2: {
                      3:   local(u,v,j,r,s);
                      4:
                      5:   if (type(D) != "t_INT" || D < 2, return(-1));
                      6:   u = sqrtint(D); v = D-u^2;
                      7:   if (!v, return(0));
                      8:   s = v;
                      9:   r = u; j = 0;
                     10:   until (u==r && v==s,
                     11:     u = (r+u)\v * v - u;
                     12:     v = (D-u^2)\v; j++;
                     13:   ); j
                     14: }