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

edit

everyone! Let $E:y^2+y=x^3-61$ be the minimal model of the elliptic curve 243b. How can I find the CM point $\tau$ in $X_0(243)$ such that $\tau$ maps to the point $(3\sqrt[3]{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 $e^{2\pi i \tau}$. After taking log and dividing by $2 \pi i$, I get the numerical $\tau$. But if I use

z=
p=z.algebraic_dependency(100)

I get the wrong polynomial. Why?