A Brief Introduction to Modular Forms Computation

Size: px
Start display at page:

Download "A Brief Introduction to Modular Forms Computation"

Transcription

1 A Brief Introduction to Modular Forms Computation Magma Supported by GCOE Program Math-For-Industry Education & Research Hub

2 What s this? Definitions and Properties Demonstration H := H P 1 (Q) some conditions k Z: f (g(z)) = (cz + d) k f (z) for all g Γ(N) Γ(N): Congruence subgroup N: with Nebentypus ε M k (Γ(N), ε): (N, k, ε) dim C M k (Γ(N)) < + f = n 0 a nq n q = e 2πiz/N : q- S k (Γ(N), ε): (N, k, ε) f M k (Γ(N), ε) s.t. a 0 = 0

3 Background Definitions and Properties Demonstration eigenform Fermat - Serre mod p Galois etc. Algebraic topology String theory Algebraic combinatorics : e.g. Kissing Number Problem

4 on Number Theory Definitions and Properties Demonstration X 0 (N) e.g. N = 39 Hecke T n q- T n f = a n f a n C (Atkin-Lehner-Li, Miyake) S k (Γ 1 (N)) = M N d (N/M) α d (Sk new (Γ 1 (M)))

5 Construct Space and Hecke Action T 2 on M 2 (Γ 0 (41)). Definitions and Properties Demonstration > M41 := ModularForms(Gamma0(41),2); > T2 := HeckeOperator(M41,2); T2; [ ] [ ] [ ] [ ] > Parent(T2); Full Matrix Algebra of degree 4 over Integer Ring > Ch2 := CharacteristicPolynomial(T2); Ch2; x^4-2*x^3-8*x^2 + 14*x +3 > Factorization(Ch2); [ <x - 3, 1>, <x^3 + x^2-5*x - 1, 1> ]

6 Compute Newforms Definitions and Properties Demonstration S new 2 (Γ 0 (11)). > S := CuspForms(Gamma0(11),2); > N := Newforms(S); N; [* [* q - 2*q^2 - q^3 + 2*q^4 + q^5 + 2*q^6-2*q^7-2*q^9-2*q^10 + q^11 + O(q^12) *], [* 5/12 + q + 3*q^2 + 4*q^3 + 7*q^4 + 6*q^5 + 12*q^6 + 8*q^7 + 15*q^8 + 13*q^9 + 18*q^10 + q^11 + O(q^12) *] *] > Newforms("G0N11k2A"); // LABELS [* [* q - 2*q^2 - q^3 + 2*q^4 + q^5 + 2*q^6-2*q^7-2*q^9-2*q^10 + q^11 + O(q^12) *] *]

7 SMC: Algebraic vs. Analytic Theorem (Khare-Wintenberger, 2007) Q detρ(c) = 1 2 mod p Galois ρ : Gal(Q/Q) GL 2 (F p ) (N(ρ), k(ρ), ε(ρ)). f = a n q n ( q = e 2πiz ) S k(ρ) (Γ, ε(ρ)) n 1

8 Verification Galois Tr(ρ(Frob l )) a l (mod p) for all l pn(ρ): prime Galois (N, k, ε) (N, k, ε) conductor, Serre weight, character level, weight, character

9 Matching example Q, Galois. E : y 2 + xy + y = x 3 + 1, E E. E Galois ρ E,l N(ρ E,l ) = p l,l ord p( E ) { 2 (l ordl ( p, k(ρ E,l ) = E )) l + 1 (otherwise), ε(ρ E,l ) = 1. Serre > E := EllipticCurve([1,0,1,0,1]); Elliptic Curve defined by y^2 + x*y + y = x^3 + 1 over Rational Field > D := Discriminant(E); D; -639 > Factorization(D); [ <3, 2>, <71, 1> ]

10 Matching example E = Galois ρ 3 = ρ E,3 (N(ρ 3 ), k(ρ 3 ), ε(ρ 3 )) = (71, 4, 1) Tr(Frob p (ρ 3 )) E [ 1, 1, 2, 2, 0, -2, 0, 0, 0, -2, -10, -6,... ], (71, 4, 1). S new 4 (Γ 0 (71)). > S71 := CuspForms(Gamma0(71),4); > f := Newforms(S71,1); > [Coefficient(f,p) : p in [1..50] IsPrime(p)]; [ 1, 1, -16, -1, 24, 7, 72, -153, -213, 232, 149, -204, -432, 71, 273 ]

11 Matching example p ρ 1 * f 1 (1) mod 3 p ρ 1 * f 1 (1) ,.

12 for Generalization Magma Hilbert / totally real case - Quaternion algebra Bianchi / imaginary quadratic case - Sharbly complex, Voronoi polyhedron. Bianchi Ver Bianchi.

13 Final Remark System for Algebra and Geometry Experimentation. W. Stein. Magma interpreter sage: magma.setdefaultrealfieldprecision(50) # magma >= v2.12; optional - magma sage: magma.eval( 1.1 ) # optional - magma (omitted)

Mazur [Ma1] Schlessinger [Sch] [SL] [Ma1] [Ma1] [Ma2] Galois [] 17 R m R R R M End R M) M R ut R M) M R R G R[G] R G Sets 1 Λ Noether Λ k Λ m Λ k C Λ

Mazur [Ma1] Schlessinger [Sch] [SL] [Ma1] [Ma1] [Ma2] Galois [] 17 R m R R R M End R M) M R ut R M) M R R G R[G] R G Sets 1 Λ Noether Λ k Λ m Λ k C Λ Galois ) 0 1 1 2 2 4 3 10 4 12 5 14 16 0 Galois Galois Galois TaylorWiles Fermat [W][TW] Galois Galois Galois 1 Noether 2 1 Mazur [Ma1] Schlessinger [Sch] [SL] [Ma1] [Ma1] [Ma2] Galois [] 17 R m R R R

More information

Tabulation of the clasp number of prime knots with up to 10 crossings

Tabulation of the clasp number of prime knots  with up to 10 crossings . Tabulation of the clasp number of prime knots with up to 10 crossings... Kengo Kawamura (Osaka City University) joint work with Teruhisa Kadokami (East China Normal University).. VI December 20, 2013

More information

●70974_100_AC009160_KAPヘ<3099>ーシス自動車約款(11.10).indb

●70974_100_AC009160_KAPヘ<3099>ーシス自動車約款(11.10).indb " # $ % & ' ( ) * +, -. / 0 1 2 3 4 5 6 7 8 9 : ; < = >? @ A B C D E F G H I J K L M N O P Q R S T U V W X Y " # $ % & ' ( ) * + , -. / 0 1 2 3 4 5 6 7 8 9 : ; < = > ? @ A B

More information

- 1 - - 2 - 320 421 928 1115 12 8 116 124 2 7 4 5 428 515 530 624 921 1115 1-3 - 100 250-4 - - 5 - - 6 - - 7 - - 8 - - 9 - & & - 11 - - 12 - GT GT - 13 - GT - 14 - - 15 - - 16 - - 17 - - 18 - - 19 - -

More information

I. (CREMONA ) : Cremona [C],., modular form f E f. 1., modular X H 1 (X, Q). modular symbol M-symbol, ( ) modular symbol., notation. H = { z = x

I. (CREMONA ) : Cremona [C],., modular form f E f. 1., modular X H 1 (X, Q). modular symbol M-symbol, ( ) modular symbol., notation. H = { z = x I. (CREMONA ) : Cremona [C],., modular form f E f. 1., modular X H 1 (X, Q). modular symbol M-symbol, ( ). 1.1. modular symbol., notation. H = z = x iy C y > 0, cusp H = H Q., Γ = PSL 2 (Z), G Γ [Γ : G]

More information

1 = = = (set) (element) a A a A a A a A a A {2, 5, (0, 1)}, [ 1, 1] = {x; 1 x 1}. (proposition) A = {x; P (x)} P (x) x x a A a A Remark. (i) {2, 0, 0,

1 = = = (set) (element) a A a A a A a A a A {2, 5, (0, 1)}, [ 1, 1] = {x; 1 x 1}. (proposition) A = {x; P (x)} P (x) x x a A a A Remark. (i) {2, 0, 0, 2005 4 1 1 2 2 6 3 8 4 11 5 14 6 18 7 20 8 22 9 24 10 26 11 27 http://matcmadison.edu/alehnen/weblogic/logset.htm 1 1 = = = (set) (element) a A a A a A a A a A {2, 5, (0, 1)}, [ 1, 1] = {x; 1 x 1}. (proposition)

More information

main.dvi

main.dvi Nim [email protected],.,.,,,.,,,,,,, Nim,.,,,,. Nim Nim,.,.,.,,.,.,. [1, 3],,, Nim,,., Nim. Date:. August 10-11, 1999 2 1 Nim.. Pile., Pile.,. normal case.,. reverse case.,.. Pile. N 1, N 2, N

More information

Armstrong culture Web

Armstrong culture Web 2004 5 10 M.A. Armstrong, Groups and Symmetry, Springer-Verlag, NewYork, 1988 (2000) (1989) (2001) (2002) 1 Armstrong culture Web 1 3 1.1................................. 3 1.2.................................

More information

k + (1/2) S k+(1/2) (Γ 0 (N)) N p Hecke T k+(1/2) (p 2 ) S k+1/2 (Γ 0 (N)) M > 0 2k, M S 2k (Γ 0 (M)) Hecke T 2k (p) (p M) 1.1 ( ). k 2 M N M N f S k+

k + (1/2) S k+(1/2) (Γ 0 (N)) N p Hecke T k+(1/2) (p 2 ) S k+1/2 (Γ 0 (N)) M > 0 2k, M S 2k (Γ 0 (M)) Hecke T 2k (p) (p M) 1.1 ( ). k 2 M N M N f S k+ 1 SL 2 (R) γ(z) = az + b cz + d ( ) a b z h, γ = SL c d 2 (R) h 4 N Γ 0 (N) {( ) } a b Γ 0 (N) = SL c d 2 (Z) c 0 mod N θ(z) θ(z) = q n2 q = e 2π 1z, z h n Z Γ 0 (4) j(γ, z) ( ) a b θ(γ(z)) = j(γ, z)θ(z)

More information

2016 Course Description of Undergraduate Seminars (2015 12 16 ) 2016 12 16 ( ) 13:00 15:00 12 16 ( ) 1 21 ( ) 1 13 ( ) 17:00 1 14 ( ) 12:00 1 21 ( ) 15:00 1 27 ( ) 13:00 14:00 2 1 ( ) 17:00 2 3 ( ) 12

More information

2014 F/ E 1 The arithmetic of elliptic curves from a viewpoint of computation 1 Shun ichi Yokoyama / JST CREST,.

2014 F/ E 1 The arithmetic of elliptic curves from a viewpoint of computation 1 Shun ichi Yokoyama / JST CREST,. 2014 F/ E 1 The arithmetic of elliptic curves from a viewpoint of computation 1 Shun ichi Yokoyama / JST CREST,. http://www2.math.kyushu-u.ac.jp/~s-yokoyama/yamagata2014.html. K Q, C, F p.,, f = 0.,,.,

More information

3 3.3. I 3.3.2. [ ] N(µ, σ 2 ) σ 2 (X 1,..., X n ) X := 1 n (X 1 + + X n ): µ X N(µ, σ 2 /n) 1.8.4 Z = X µ σ/ n N(, 1) 1.8.2 < α < 1/2 Φ(z) =.5 α z α

3 3.3. I 3.3.2. [ ] N(µ, σ 2 ) σ 2 (X 1,..., X n ) X := 1 n (X 1 + + X n ): µ X N(µ, σ 2 /n) 1.8.4 Z = X µ σ/ n N(, 1) 1.8.2 < α < 1/2 Φ(z) =.5 α z α 2 2.1. : : 2 : ( ): : ( ): : : : ( ) ( ) ( ) : ( pp.53 6 2.3 2.4 ) : 2.2. ( ). i X i (i = 1, 2,..., n) X 1, X 2,..., X n X i (X 1, X 2,..., X n ) ( ) n (x 1, x 2,..., x n ) (X 1, X 2,..., X n ) : X 1,

More information

2.1 H f 3, SL(2, Z) Γ k (1) f H (2) γ Γ f k γ = f (3) f Γ \H cusp γ SL(2, Z) f k γ Fourier f k γ = a γ (n)e 2πinz/N n=0 (3) γ SL(2, Z) a γ (0) = 0 f c

2.1 H f 3, SL(2, Z) Γ k (1) f H (2) γ Γ f k γ = f (3) f Γ \H cusp γ SL(2, Z) f k γ Fourier f k γ = a γ (n)e 2πinz/N n=0 (3) γ SL(2, Z) a γ (0) = 0 f c GL 2 1 Lie SL(2, R) GL(2, A) Gelbart [Ge] 1 3 [Ge] Jacquet-Langlands [JL] Bump [Bu] Borel([Bo]) ([Ko]) ([Mo]) [Mo] 2 2.1 H = {z C Im(z) > 0} Γ SL(2, Z) Γ N N Γ (N) = {γ SL(2, Z) γ = 1 2 mod N} g SL(2,

More information

5.. z = f(x, y) y y = b f x x g(x) f(x, b) g x ( ) A = lim h g(a + h) g(a) h g(x) a A = g (a) = f x (a, b)............................................

5.. z = f(x, y) y y = b f x x g(x) f(x, b) g x ( ) A = lim h g(a + h) g(a) h g(x) a A = g (a) = f x (a, b)............................................ 5 partial differentiation (total) differentiation 5. z = f(x, y) (a, b) A = lim h f(a + h, b) f(a, b) h........................................................... ( ) f(x, y) (a, b) x A (a, b) x (a, b)

More information

... 3... 3... 3... 3... 4... 7... 10... 10... 11... 12... 12... 13... 14... 15... 18... 19... 20... 22... 22... 23 2

... 3... 3... 3... 3... 4... 7... 10... 10... 11... 12... 12... 13... 14... 15... 18... 19... 20... 22... 22... 23 2 1 ... 3... 3... 3... 3... 4... 7... 10... 10... 11... 12... 12... 13... 14... 15... 18... 19... 20... 22... 22... 23 2 3 4 5 6 7 8 9 Excel2007 10 Excel2007 11 12 13 - 14 15 16 17 18 19 20 21 22 Excel2007

More information

2 probably 3 probability theory probability theory (gàil`ü) 2 1960 2.1, 1: 583 589 666 1853 2967 2 3 1973

2 probably 3 probability theory probability theory (gàil`ü) 2 1960 2.1, 1: 583 589 666 1853 2967 2 3 1973 ( ) ( C) 1 probability probability probable probable probable probable probably probably maybe perhaps possibly likely possibly

More information

1 2 1.1............................................ 3 1.2.................................... 7 1.3........................................... 9 1.4..

1 2 1.1............................................ 3 1.2.................................... 7 1.3........................................... 9 1.4.. 2010 8 3 ( ) 1 2 1.1............................................ 3 1.2.................................... 7 1.3........................................... 9 1.4........................................

More information

1 n =3, 2 n 3 x n + y n = z n x, y, z 3 a, b b = aq q a b a b b a b a a b a, b a 0 b 0 a, b 2

1 n =3, 2 n 3 x n + y n = z n x, y, z 3 a, b b = aq q a b a b b a b a a b a, b a 0 b 0 a, b 2 n =3, 200 2 10 1 1 n =3, 2 n 3 x n + y n = z n x, y, z 3 a, b b = aq q a b a b b a b a a b a, b a 0 b 0 a, b 2 a, b (a, b) =1a b 1 x 2 + y 2 = z 2, (x, y) =1, x 0 (mod 2) (1.1) x =2ab, y = a 2 b 2, z =

More information

Fig. 3 Flow diagram of image processing. Black rectangle in the photo indicates the processing area (128 x 32 pixels).

Fig. 3 Flow diagram of image processing. Black rectangle in the photo indicates the processing area (128 x 32 pixels). Fig. 1 The scheme of glottal area as a function of time Fig. 3 Flow diagram of image processing. Black rectangle in the photo indicates the processing area (128 x 32 pixels). Fig, 4 Parametric representation

More information

一般演題(ポスター)

一般演題(ポスター) 6 5 13 : 00 14 : 00 A μ 13 : 00 14 : 00 A β β β 13 : 00 14 : 00 A 13 : 00 14 : 00 A 13 : 00 14 : 00 A β 13 : 00 14 : 00 A β 13 : 00 14 : 00 A 13 : 00 14 : 00 A β 13 : 00 14 : 00 A 13 : 00 14 : 00 A

More information

Sage for Mathematics : a Primer ‚æ1Łfl - Sage ‡ð™m‡é

Sage for Mathematics : a Primer   ‚æ1Łfl - Sage ‡ð™m‡é .. / JST CREST. [email protected] 2013 1 15 / IIJ Shun ichi Yokoyama (IMI/JST CREST) at IIJ January 15th, 2013 1 / 23 1 Sage Sage Sage Notebook Sage Salvus 2 Sage Sage Cryptosystem Sage Shun

More information

untitled

untitled 20 7 1 22 7 1 1 2 3 7 8 9 10 11 13 14 15 17 18 19 21 22 - 1 - - 2 - - 3 - - 4 - 50 200 50 200-5 - 50 200 50 200 50 200 - 6 - - 7 - () - 8 - (XY) - 9 - 112-10 - - 11 - - 12 - - 13 - - 14 - - 15 - - 16 -

More information

untitled

untitled 19 1 19 19 3 8 1 19 1 61 2 479 1965 64 1237 148 1272 58 183 X 1 X 2 12 2 15 A B 5 18 B 29 X 1 12 10 31 A 1 58 Y B 14 1 25 3 31 1 5 5 15 Y B 1 232 Y B 1 4235 14 11 8 5350 2409 X 1 15 10 10 B Y Y 2 X 1 X

More information

http://banso.cocolog-nifty.com/ 100 100 250 5 1 1 http://www.banso.com/ 2009 5 2 10 http://www.banso.com/ 2009 5 2 http://www.banso.com/ 2009 5 2 http://www.banso.com/ < /> < /> / http://www.banso.com/

More information

R R P N (C) 7 C Riemann R K ( ) C R C K 8 (R ) R C K 9 Riemann /C /C Riemann 10 C k 11 k C/k 12 Riemann k Riemann C/k k(c)/k R k F q Riemann 15

R R P N (C) 7 C Riemann R K ( ) C R C K 8 (R ) R C K 9 Riemann /C /C Riemann 10 C k 11 k C/k 12 Riemann k Riemann C/k k(c)/k R k F q Riemann 15 (Gen KUROKI) 1 1 : Riemann Spec Z 2? 3 : 4 2 Riemann Riemann Riemann 1 C 5 Riemann Riemann R compact R K C ( C(x) ) K C(R) Riemann R 6 (E-mail address: [email protected]) 1 1 ( 5 ) 2 ( Q ) Spec

More information

x 3 a (mod p) ( ). a, b, m Z a b m a b (mod m) a b m 2.2 (Z/mZ). a = {x x a (mod m)} a Z m 0, 1... m 1 Z/mZ = {0, 1... m 1} a + b = a +

x 3 a (mod p) ( ). a, b, m Z a b m a b (mod m) a b m 2.2 (Z/mZ). a = {x x a (mod m)} a Z m 0, 1... m 1 Z/mZ = {0, 1... m 1} a + b = a + 1 1 22 1 x 3 (mod ) 2 2.1 ( )., b, m Z b m b (mod m) b m 2.2 (Z/mZ). = {x x (mod m)} Z m 0, 1... m 1 Z/mZ = {0, 1... m 1} + b = + b, b = b Z/mZ 1 1 Z Q R Z/Z 2.3 ( ). m {x 0, x 1,..., x m 1 } modm 2.4

More information

n=360 28.6% 34.4% 36.9% n=360 2.5% 17.8% 19.2% n=64 0.8% 0.3% n=69 1.7% 3.6% 0.6% 1.4% 1.9% < > n=218 1.4% 5.6% 3.1% 60.6% 0.6% 6.9% 10.8% 6.4% 10.3% 33.1% 1.4% 3.6% 1.1% 0.0% 3.1% n=360 0% 50%

More information

2 Riemann Im(s) > 0 ζ(s) s R(s) = 2 Riemann [Riemann]) ζ(s) ζ(2) = π2 6 *3 Kummer s = 2n, n N ζ( 2) = 2 2, ζ( 4) =.3 2 3, ζ( 6) = ζ( 8)

2 Riemann Im(s) > 0 ζ(s) s R(s) = 2 Riemann [Riemann]) ζ(s) ζ(2) = π2 6 *3 Kummer s = 2n, n N ζ( 2) = 2 2, ζ( 4) =.3 2 3, ζ( 6) = ζ( 8) (Florian Sprung) p 2 p * 9 3 p ζ Mazur Wiles 4 5 6 2 3 5 2006 http://www.icm2006.org/video/ eighth session [ ] Coates [Coates] 2 735 Euler n n 2 = p p 2 p 2 = π2 6 859 Riemann ζ(s) = n n s = p p s s ζ(s)

More information

P-12 P-13 3 4 28 16 00 17 30 P-14 P-15 P-16 4 14 29 17 00 18 30 P-17 P-18 P-19 P-20 P-21 P-22

P-12 P-13 3 4 28 16 00 17 30 P-14 P-15 P-16 4 14 29 17 00 18 30 P-17 P-18 P-19 P-20 P-21 P-22 1 14 28 16 00 17 30 P-1 P-2 P-3 P-4 P-5 2 24 29 17 00 18 30 P-6 P-7 P-8 P-9 P-10 P-11 P-12 P-13 3 4 28 16 00 17 30 P-14 P-15 P-16 4 14 29 17 00 18 30 P-17 P-18 P-19 P-20 P-21 P-22 5 24 28 16 00 17 30 P-23

More information

項 目

項 目 1 1 2 3 11 4 6 5 7,000 2 120 1.3 4,000 04 450 < > 5 3 6 7 8 9 4 10 11 5 12 45 6 13 E. 7 B. C. 14 15 16 17 18 19 20 21 22 23 8 24 25 9 27 2 26 6 27 3 1 3 3 28 29 30 9 31 32 33 500 1 4000 0 2~3 10 10 34

More information

32 1 7 1 20 ( ) [18 30] [21 00] 2 3 ( ( ) ) ( ) 4 1 2 95 ( 7 3 2 ) 1 2 3 2 a b

32 1 7 1 20 ( ) [18 30] [21 00] 2 3 ( ( ) ) ( ) 4 1 2 95 ( 7 3 2 ) 1 2 3 2 a b 31 1 7 1 ( ) [18 35] [20 40] 2 3 ( ( ) ) ( ) 4 1 2 1 22 1 2 a b T T A c d 32 1 7 1 20 ( ) [18 30] [21 00] 2 3 ( ( ) ) ( ) 4 1 2 95 ( 7 3 2 ) 1 2 3 2 a b 7 1 ( 34 ) 1 7 2 13 ( ) [18 20] [20 00] 2 4 7 3

More information

( ) 1., ([SU] ): F K k., Z p -, (cf. [Iw2], [Iw3], [Iw6]). K F F/K Z p - k /k., Weil., K., K F F p- ( 4.1).,, Z p -,., Weil..,,. Weil., F, F projectiv

( ) 1., ([SU] ): F K k., Z p -, (cf. [Iw2], [Iw3], [Iw6]). K F F/K Z p - k /k., Weil., K., K F F p- ( 4.1).,, Z p -,., Weil..,,. Weil., F, F projectiv ( ) 1 ([SU] ): F K k Z p - (cf [Iw2] [Iw3] [Iw6]) K F F/K Z p - k /k Weil K K F F p- ( 41) Z p - Weil Weil F F projective smooth C C Jac(C)/F ( ) : 2 3 4 5 Tate Weil 6 7 Z p - 2 [Iw1] 2 21 K k k 1 k K

More information

サービス付き高齢者向け住宅賠償責任保険.indd

サービス付き高齢者向け住宅賠償責任保険.indd 1 2 1 CASE 1 1 2 CASE 2 CASE 3 CASE 4 3 CASE 5 4 3 4 5 6 2 CASE 1 CASE 2 CASE 3 7 8 3 9 10 CASE 1 CASE 2 CASE 3 CASE 4 11 12 13 14 1 1 2 FAX:03-3375-8470 2 3 3 4 4 3 15 16 FAX:03-3375-8470 1 2 0570-022808

More information

Pari-gp 2006/7/12 2 1. Pari-gp 2. Microsoft Windows 3. 4. Pari-gp 5. 2 6. 7. Galois 8.

Pari-gp 2006/7/12 2 1. Pari-gp 2. Microsoft Windows 3. 4. Pari-gp 5. 2 6. 7. Galois 8. Pari-gp 2006/7/12 1 Pari-GP 2006 7 12 Pari-gp 2006/7/12 2 1. Pari-gp 2. Microsoft Windows 3. 4. Pari-gp 5. 2 6. 7. Galois 8. Pari-gp 2006/7/12 2-1 +ε. Pari-gp 2006/7/12 3 Pari-gp Pari C gp GP gp K. Belabas

More information