Revision history [back]

this would get you close:

E.torsion_polynomial(n).splitting_field()

You're just some quadratic extensions away from the desired field (you just need to make sure that roots of the division polynomial are indeed x-coordinates of points).

Be careful, though: generally such explicit representations of splitting fields are entirely unworkable because they tend to have very high degrees.