2002.sagews

prime_pi(1000000)

78498

E  = EllipticCurve('389a')

P, Q = E.gens()

7*P + Q

(10074900482821/4690473062500 : 27802875709438987131/10158392035109375000 : 1)