\p 38
\e
ellinit([0,0,0,-1,0])
ellinit([0,0,0,-17,0],1)
ellj(I)
ellsub(ellinit([0,0,0,-17,0]),[-1,4],[-4,2])
\\
acurve=ellinit([0,0,1,-1,0])
apoint=[2,2]
elladd(acurve,apoint,apoint)
ellak(acurve,1000000007)
ellan(acurve,100)
ellap(acurve,10007)
ellap(acurve,10007,1)
deu=direuler(p=2,100,1/(1-ellap(acurve,p)*x+if(acurve[12]%p,p,0)*x^2))
ellan(acurve,100)==deu
ellisoncurve(acurve,apoint)
acurve=ellchangecurve(acurve,[-1,1,2,3])
apoint=ellchangepoint(apoint,[-1,1,2,3])
ellisoncurve(acurve,apoint)
ellglobalred(acurve)
ellheight(acurve,apoint)
ellheight(acurve,apoint,1)
ellordinate(acurve,1)
ellztopoint(acurve,ellpointtoz(acurve,apoint))
ellpow(acurve,apoint,10)
ellwp(acurve)
ellpointtoz(acurve,apoint)
q*Ser(ellan(acurve,100),q)
\\
bcurve=ellinit([0,0,0,-3,0])
elllocalred(bcurve,2)
elltaniyama(bcurve)
\\
ccurve=ellinit([0,0,-1,-1,0])
l=elllseries(ccurve,2)
elllseries(ccurve,2,1.2)-l
\\
tcurve=ellinit([1,0,1,-19,26]);
ellorder(tcurve,[1,2])
elltors(tcurve)
\\
mcurve=ellinit([0,0,0,-17,0])
mpoints=[[-1,4],[-4,2]]~
mhbi=ellbil(mcurve,mpoints,[9,24])
ma=ellheightmatrix(mcurve,mpoints)
matsolve(ma,mhbi)
\\
cmcurve=ellinit([0,-3/4,0,-2,-1])
ellpow(cmcurve,[x,y],quadgen(-7))
\p 96
precision(cmcurve)
getheap
print("Total time spent: ",gettime);
\q