\\ Pari/GP code for working with elliptic curve isogeny class 3249.b


\\ Define the isogeny class: 
E = ellinit([0, 0, 1, -7581, 279504])
ellisomat(E)

\\ Mordell-Weil rank: 
[lower,upper] = ellrank(E)

\\ q-expansion of modular form: 
Ser(ellan(E,20),q)*q

\\ Isogeny matrix: 
ellisomat(E)