3233.sagews

1+1

2

E = EllipticCurve([0,1])

E.weierstrass_p()

z^-2 - 1/7*z^4 + 1/637*z^10 - 1/84721*z^16 + O(z^20)