discriminant(t^3 - (12m^2n^2-(m^2-n^2)^2)t + 4 m n(m^2+n^2)^2, t)