3723.sagews

A = 2811483865363529644328185773079724833947000934911780977075879578375099762840
B = 250446441104874779541582568787395262709240144181428082354092935623339223922413
p = 745716355583897279336091334816026338459826430674460685108127047713267472241559
ec = EllipticCurve(GF(p), [A, B])

ec.order()

745716355583897279336091334816026338459826098367461738879158821761502402155416

ec.j_invariant()

376738357847214901164140665055058431119917844710765744434849749843642790706

x0 = 332306998946228968225951765070086144
y0 = 60890339900407113676354129625720714119895503851481483716090698020108489580573
ecp = ec.point([x0, y0])

ecp.order()

186429088895974319834022833704006584614956524591865434719789705440375600538854

ecp2 = ecp + ecp

ecp2.order()

93214544447987159917011416852003292307478262295932717359894852720187800269427