syms t a b c d f_0 f_1
eq1=t*((a*t+b)^2+f_0^2*(c*t+d)^2)^2-(a*d-b*c)*(1+f_1^2*t^2)^2*(a*t+b)*(c*t+d);
eq2=expand(eq1);
A=poly_coef(eq2,t)
A =
[ a*f_1^4*c*(-a*d+b*c), c^4*f_0^4+b^2*c^2*f_1^4+a^4-a^2*d^2*f_1^4+2*a^2*f_0^2*c^2, -2*a^2*d*f_1^2*c+2*b*c^2*f_1^2*a+4*c*f_0^2*d*a^2+4*f_0^4*c^3*d+4*c^2*f_0^2*b*a+4*b*a^3+b^2*c*f_1^4*d-a*d^2*f_1^4*b, 2*b^2*f_0^2*c^2+6*b^2*a^2+8*c*f_0^2*d*b*a+6*f_0^4*c^2*d^2+2*b^2*c^2*f_1^2+2*a^2*f_0^2*d^2-2*a^2*d^2*f_1^2, c^2*b*a-a^2*d*c+2*b^2*c*f_1^2*d+4*c*f_0^2*d*b^2+4*a*b^3+4*f_0^2*d^2*b*a-2*a*d^2*f_1^2*b+4*c*f_0^4*d^3, 2*b^2*f