sage: alpha=sqrt(-17/2 + sqrt(17^2-4(17))/2)
sage: f=alpha.minpoly()
sage: K.<a>=NumberField(x^4+17*x^2+17)
sage: R<y>=PolynomialRing(K)
sage: R(f).factor()
(y-a)(y+a)(y^2+a^2+17)
sage: g=y^2+a^2+17
sage: L.<b>=K.extension(g)
sage: L.relative_discriminant()
Fractional ideal (1)