# SageMath code for working with elliptic curve isogeny class 3.3.148.1-400.2-a
# (Note that not all these functions may be available, and some may take a long time to execute.)


# Define the base number field: 
R.<x> = PolynomialRing(QQ);  K.<a> = NumberField(R([1, -3, -1, 1]))

# Define the isogeny class: 
E = EllipticCurve([K([-1,0,1]),K([1,0,-1]),K([1,1,0]),K([-2052,5213,4455]),K([16900486,-42913527,-36675517])])
E.isogeny_class()

# Compute the Mordell-Weil rank: 
E.rank()

# Isogeny matrix: 
E.isogeny_class().matrix()

# Isogeny graph: 
E.isogeny_class().graph().plot(edge_labels=True)

# List of curves in the isogeny class: 
E.isogeny_class().curves