benchmarks/timebench2c.gp

   1 n=100
   2 lm=9
   3 m=10^lm
   4 \precision=lm+5
   5 root3=sqrt(3)
   6 root5=sqrt(5)
   7 p1 = sum(0,j=0,2*n,divres(floor(root5*m*j),m)[2]*(-1)^j*x^j)
   8 p2 = sum(0,j=0,n,divres(floor(root3*m*j),m)[2]*x^j)
   9 #
  10 p1*p2
  11 p1*p2
  12 divres(p1,p2)
  13 divres(p1,p2)
  14 gcd(p1,p2)
  15 gcd(p1,p2)