OpenXM_contrib/pari-2.2/examples/contfrac.gp - view

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, 8 months ago) by noro
Branch: NORO
CVS Tags: RELEASE_1_2_3, RELEASE_1_2_2_KNOPPIX_b, RELEASE_1_2_2_KNOPPIX, RELEASE_1_2_2, RELEASE_1_2_1, PARI_2_2
Changes since 1.1: +0 -0 lines

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
}