crypto.sagews

E = EllipticCurve(GF(2^7-1),[0,1])

P = E.random_element()
Q = E.random_element()

g = E.plot()
g += points([P[:2],Q[:2]],size=90,color="red")
g += points([(E(P)+E(Q))[:2]],size=90,color="green")
g.show(title="$E/\\mathbb{F}_{%d}: %s = 0$" % (E.base_field().cardinality(),latex(E.defining_polynomial()(z=1))))