Q49-49a1-49a2.gp (846B)
1 elltorsfield_deg( E, n, deg ) = { 2 my( f, eqn, res ); 3 f = elldivpol( E, n ); 4 eqn = y^2 + E[1]*x*y + E[3]*y - x^3 - E[2]*x^2 - E[4]*x - E[5]; 5 res = substpol( polresultant( f, eqn ), y, x ); 6 printf( "Computing %d-torsion field...\n", n ); 7 return( nfinit( nfsplitting( f*res, deg ) ) ); 8 } 9 10 print("Computing 49-torsion field for 49a1, expecting maximal degree 42*49\n"); 11 elltorsfield_deg( ellinit( "49a1" ), 49, 42*49 ); 12 13 print("Computing 49-torsion field for 49a1, expecting maximal degree 42*49\n"); 14 elltorsfield_deg( ellinit( "49a1" ), 49, 42*49 );