# SageMath code for working with elliptic curve isogeny class 66066.cq


# Define the isogeny class: 
E = EllipticCurve([1, 0, 0, -6987813, 7102721385])
E.isogeny_class()

# Mordell-Weil rank: 
E.rank()

# q-expansion of modular form: 
E.q_eigenform(20)

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

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

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