|
MATHEMATICA
|
Discriminant[p_?PolynomialQ, x_] := With[{n = Exponent[p, x], k = Exponent[D[p, x], x]}, Cancel[((-1)^(n(n - 1)/2)Resultant[ p, D[p, x], x]) Coefficient[p, x, n]^(n - k - 2)]] ; Table[Discriminant[x^5 - 10p x^3 + 45p^2 x - p^2, x], {p, 1, 20}]
|