File: timebench2c.gp

package info (click to toggle)
cln 0.98-7.1
  • links: PTS
  • area: main
  • in suites: slink
  • size: 12,188 kB
  • ctags: 15,282
  • sloc: cpp: 71,545; ansic: 12,015; asm: 8,431; sh: 3,159; makefile: 886; lisp: 64
file content (15 lines) | stat: -rw-r--r-- 240 bytes parent folder | download | duplicates (6) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
 
n=100
lm=9
m=10^lm
\precision=lm+5
root3=sqrt(3)
root5=sqrt(5)
p1 = sum(0,j=0,2*n,divres(floor(root5*m*j),m)[2]*(-1)^j*x^j)
p2 = sum(0,j=0,n,divres(floor(root3*m*j),m)[2]*x^j)
#
p1*p2
p1*p2
divres(p1,p2)
divres(p1,p2)
gcd(p1,p2)
gcd(p1,p2)