\\package:seadata
ellmodulareqn(5)
do(i,v)=
{
  E = ellinit([0,0,0,v[2],v[3]]*Mod(1,v[1]));
  print(i, ": ", ellap(E,v[1]));
}
{
v=[[202600005603433095160409308644759862837,25496852782325453225973142890909600552,129550610797481291887769966647995045232],
[173327739907566197112155895875385467119,52716988591102938437323369716512206005,43087597392844950895070462564402654315],
[523583591335747530615071369664554118036421993253,25679429559575246581833628827363203226862930934,78220741356817481602535950765825830003112603824],
[1319450668936329467137913739322239157303860926441,807652438980115949692649657326438677571309575087,1017125626316888896817395440127041355136940446205],
[439581010348913995032270658729785287035964480323270935583,323922016281172901245123590380881598241426088528431020005,54496426275749371996644207660602248980615186517525561222],
[2979720374579183569554262247145622188470249961843364603751,428654869348535206084607029945439317783967748844874233571,1182279475380088064870625220629639405548336474256523329003],
[6243380271698146227966925307851825694742847655729810693741,4068721281680536125235363885580194460678653324971583338307,2519148351962491328666249705249360758373031631978108875818],
[1606938044258990275541962092341162602522202993782792835304761,1,252199199707645577897249048746397012330572101453777389069968],
[1267650600228229401496703205953,1,417990942431022911086532367249],
[590295810358705651741,1,3],
[18446744073709551629,1,42]];
}
for(i=1,#v, do(i,v[i]));
g=ffgen(nextprime(2^32)^2);E=ellinit([g,3]);ellap(E)
