# SageMath code for working with elliptic curve isogeny class 33282.d # Define the isogeny class: E = EllipticCurve([1, -1, 0, -12408, 535094]) 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