Return to bench.gp CVS log | Up to [local] / OpenXM_contrib / pari / examples |
1.1 ! maekawa 1: { ! 2: u=v=p=q=1; ! 3: for (k=1, 1000, ! 4: w = u+v; u = v; v = w; ! 5: p *= w; q = lcm(q,w); ! 6: if (k%50 == 0, ! 7: print(k " " log(p)/log(q)) ! 8: ) ! 9: ) ! 10: }