2348.sagews

EllipticCurve([0,-1]).integral_points(both_signs=True)

[(1 : 0 : 1)]

EllipticCurve([0,1]).integral_points(both_signs=True)

[(-1 : 0 : 1), (0 : -1 : 1), (0 : 1 : 1), (2 : -3 : 1), (2 : 3 : 1)]