# SageMath code for working with elliptic curve isogeny class 47040.fd


# Define the isogeny class: 
E = EllipticCurve([0, 1, 0, 159, -225])
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
