# SageMath code for working with elliptic curve isogeny class 279414.bq # Define the isogeny class: E = EllipticCurve([1, -1, 1, -4001, 98385]) 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