everyone! Let E:y2+y=x3−61 be the minimal model of the elliptic curve 243b. How can I find the CM point τ in X0(243) such that τ maps to the point (33√3,4) under the modular parametrization? Can anyone tell me the answer or how to use sagemath to find it?
I use the sagemath code
EllipticCurve([0,0,1,0,-61])
phi = EllipticCurve([0,0,1,0,-61]).modular_parametrization()
f=phi.power_series(prec = 10000)[1]
f.truncate(20000)
to get the parametrization of y coordinate, then I use
q=var('q')
f(q)=
df=diff(f,q)
NewtonIt(q)=q-(f/df)(q)
xn=e^(2*pi*I*a/20.031)
for i in range(1000):
xn=N(NewtonIt(xn),digits=2000)
print xn
to get the numerical e2πiτ. After taking log and dividing by 2πi, I get the numerical τ. But if I use
z=
p=z.algebraic_dependency(100)
I get the wrong polynomial. Why?