Processing math: 100%
Ask Your Question

Revision history [back]

click to hide/show revision 1
initial version

asked 5 years ago

Lee gravatar image

How to find a CM point with the image in the elliptic curve under modular parametrization given

everyone! Let E:y2+y=x361 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 (333,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?

click to hide/show revision 2
retagged

updated 5 years ago

FrédéricC gravatar image

How to find a CM point with the image in the elliptic curve under modular parametrization given

everyone! Let E:y2+y=x361 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 (333,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?