\\ Pari/GP code for working with elliptic curve isogeny class 63870.p \\ Define the isogeny class: E = ellinit([1, 0, 0, -40, -400]) ellisomat(E) \\ Mordell-Weil rank: [lower,upper] = ellrank(E) \\ q-expansion of modular form: Ser(ellan(E,20),q)*q \\ Isogeny matrix: ellisomat(E)