1 | initial version |
Like that ?
sage: x,y = polygens(QQ,'x,y')
sage: E = EllipticCurve(-y**2+(x+1540)*(x-508)*(x-65024))
sage: E.integral_points()
[(-1540 : 0 : 1),
(-508 : 262128 : 1),
(508 : 0 : 1),
(65024 : 0 : 1),
(130556 : 33552384 : 1)]
sage: P = E.gens()[0]; P
(100132172429824908929/1508738252550400 : 143687064081412107244001809983/58603135399923860992000 : 1)
sage: P*2
(23926617756252939195845889933899553715162280807826884830539689900801/28781599829126697172901832423192736087492173270011136379110400 : -112435953838617002324224529091205045118375173739751380123760462569961416430316858965047000875038947199/154408923747010694464616627108267136335941190834695015240452304070119936139205075365384192000 : 1)