default(realprecision,38);
rk(x)=x=ellinit(x);ellanalyticrank(x);
rk([0, -1, 1, -10, -20])
rk([0, 0, 1, -1, 0])
rk([0, 1, 1, -2, 0])
rk([0, 0, 1, -7, 6])
rk([-5187, 176830])

he(x)=x=ellinit(x);ellheegner(x);
he([1, 1, 0, -1297, -18530])
he([0, -1, 1, -33, 93])
he([-157^2,0])
he([0,0,-9/484,0,-27/234256])
getheap()[1]

E=ellinit([0,-1437004800,0,458885065605120000,0]);
ellglobalred(E); ellheegner(E)

L = ellinit([0,0,1,7,6]);
ellL1(L)
ellL1(L,1)
ellL1(L,3)
ellL1(ellinit([0,1,1,-2,0]),2)
default(realprecision,115);
ellL1(L)
ellL1(L,1)
ellL1(L,3)
ellL1(ellinit([0,1,1,-2,0]),2)

default(realprecision,19);
ellmoddegree(ellinit([0,1,0,-4,-4]))
ellmoddegree(ellinit([-4,0]))
ellmoddegree(ellinit([0,-1,0,4,-4]))
ellmoddegree(ellinit([0,-1,0,-8,-16]))
ellmoddegree(ellinit([0,0,1,0,-7]))
ellmoddegree(ellinit([1,-1,0,0,-5]))
ellmoddegree(ellinit([1,-1,0,-69,-208]))
localbitprec(9); ellmoddegree(ellinit([0, -1, 1, 0, 0]))
