| Hosted by CoCalc | Download
K.<a> = NumberField(x^2+1)
R.<y> = K[] L.<b> = K.extension(y^2 + 2*a*y + 1) L
Number Field in b with defining polynomial y^2 + 2*a*y + 1 over its base field
b.charpoly()
x^2 + 2*a*x + 1