matlab用求特征值求方程,MATLAB如何求五阶符号函数特征值和特征向量

我希望求一个五阶符号函数的特征值和特征向量，将符号函数对角化。但得到的解非常长。请问错在哪里，如何计算特征值的精确解呢？谢谢解答！

syms r13 r1 r14 r15 r23 r24 r25 f1 f2 f3 f4 f5 s q;

A=[f1^2*s^2 0 r13*f1*f3*s*q r14*f1*f4*s*q r15*f1*f5*s*q

0 f2^2*s^2 r23*f2*f3*s*q r24*f2*f4*s*q r25*f2*f5*s*q

r13*f1*f3*s*q r23*f2*f3*s*q f3^2*q^2 0 0

r14*f1*f4*s*q r24*f2*f4*s*q 0 f4^2*q^2 0

r15*f1*f5*s*q r25*f2*f5*s*q 0 0 f5^2*q^2];

e=eig(A)复制代码

得到的解:

e =

root(z1^5 - z1^4*(f3^2*q^2 + f4^2*q^2 + f5^2*q^2 + f1^2*s^2 + f2^2*s^2) + z1^3*(f3^2*f4^2*q^4 + f3^2*f5^2*q^4 + f4^2*f5^2*q^4 + f1^2*f2^2*s^4 + f1^2*f3^2*q^2*s^2 + f1^2*f4^2*q^2*s^2 + f2^2*f3^2*q^2*s^2 + f1^2*f5^2*q^2*s^2 + f2^2*f4^2*q^2*s^2 + f2^2*f5^2*q^2*s^2 - f1^2*f3^2*q^2*r13^2*s^2 - f1^2*f4^2*q^2*r14^2*s^2 - f1^2*f5^2*q^2*r15^2*s^2 - f2^2*f3^2*q^2*r23^2*s^2 - f2^2*f4^2*q^2*r24^2*s^2 - f2^2*f5^2*q^2*r25^2*s^2) + z1^2*(- f3^2*f4^2*f5^2*q^6 - f1^2*f2^2*f3^2*q^2*s^4 - f1^2*f2^2*f4^2*q^2*s^4 - f1^2*f2^2*f5^2*q^2*s^4 - f1^2*f3^2*f4^2*q^4*s^2 - f1^2*f3^2*f5^2*q^4*s^2 - f2^2*f3^2*f4^2*q^4*s^2 - f1^2*f4^2*f5^2*q^4*s^2 - f2^2*f3^2*f5^2*q^4*s^2 - f2^2*f4^2*f5^2*q^4*s^2 + f1^2*f2^2*f3^2*q^2*r13^2*s^4 + f1^2*f2^2*f4^2*q^2*r14^2*s^4 + f1^2*f3^2*f4^2*q^4*r13^2*s^2 + f1^2*f3^2*f4^2*q^4*r14^2*s^2 + f1^2*f3^2*f5^2*q^4*r13^2*s^2 + f1^2*f2^2*f5^2*q^2*r15^2*s^4 + f1^2*f3^2*f5^2*q^4*r15^2*s^2 + f1^2*f4^2*f5^2*q^4*r14^2*s^2 + f1^2*f4^2*f5^2*q^4*r15^2*s^2 + f1^2*f2^2*f3^2*q^2*r23^2*s^4 + f1^2*f2^2*f4^2*q^2*r24^2*s^4 + f2^2*f3^2*f4^2*q^4*r23^2*s^2 + f1^2*f2^2*f5^2*q^2*r25^2*s^4 + f2^2*f3^2*f4^2*q^4*r24^2*s^2 + f2^2*f3^2*f5^2*q^4*r23^2*s^2 + f2^2*f3^2*f5^2*q^4*r25^2*s^2 + f2^2*f4^2*f5^2*q^4*r24^2*s^2 + f2^2*f4^2*f5^2*q^4*r25^2*s^2) - z1*(- f1^2*f2^2*f3^2*f4^2*q^4*s^4 - f1^2*f2^2*f3^2*f5^2*q^4*s^4 - f1^2*f2^2*f4^2*f5^2*q^4*s^4 - f1^2*f3^2*f4^2*f5^2*q^6*s^2 - f2^2*f3^2*f4^2*f5^2*q^6*s^2 + f1^2*f2^2*f3^2*f4^2*q^4*r13^2*s^4 + f1^2*f2^2*f3^2*f4^2*q^4*r14^2*s^4 + f1^2*f2^2*f3^2*f5^2*q^4*r13^2*s^4 + f1^2*f2^2*f3^2*f5^2*q^4*r15^2*s^4 + f1^2*f2^2*f4^2*f5^2*q^4*r14^2*s^4 + f1^2*f3^2*f4^2*f5^2*q^6*r13^2*s^2 + f1^2*f2^2*f4^2*f5^2*q^4*r15^2*s^4 + f1^2*f3^2*f4^2*f5^2*q^6*r14^2*s^2 + f1^2*f3^2*f4^2*f5^2*q^6*r15^2*s^2 + f1^2*f2^2*f3^2*f4^2*q^4*r23^2*s^4 + f1^2*f2^2*f3^2*f4^2*q^4*r24^2*s^4 + f1^2*f2^2*f3^2*f5^2*q^4*r23^2*s^4 + f1^2*f2^2*f3^2*f5^2*q^4*r25^2*s^4 + f1^2*f2^2*f4^2*f5^2*q^4*r24^2*s^4 + f1^2*f2^2*f4^2*f5^2*q^4*r25^2*s^4 + f2^2*f3^2*f4^2*f5^2*q^6*r23^2*s^2 + f2^2*f3^2*f4^2*f5^2*q^6*r24^2*s^2 + f2^2*f3^2*f4^2*f5^2*q^6*r25^2*s^2 - f1^2*f2^2*f3^2*f4^2*q^4*r13^2*r24^2*s^4 - f1^2*f2^2*f3^2*f4^2*q^4*r14^2*r23^2*s^4 - f1^2*f2^2*f3^2*f5^2*q^4*r13^2*r25^2*s^4 - f1^2*f2^2*f3^2*f5^2*q^4*r15^2*r23^2*s^4 - f1^2*f2^2*f4^2*f5^2*q^4*r14^2*r25^2*s^4 - f1^2*f2^2*f4^2*f5^2*q^4*r15^2*r24^2*s^4 + 2*f1^2*f2^2*f3^2*f4^2*q^4*r13*r14*r23*r24*s^4 + 2*f1^2*f2^2*f3^2*f5^2*q^4*r13*r15*r23*r25*s^4 + 2*f1^2*f2^2*f4^2*f5^2*q^4*r14*r15*r24*r25*s^4) + 2*f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r14*r15*r24*r25*s^4 + 2*f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r13*r15*r23*r25*s^4 + 2*f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r13*r14*r23*r24*s^4 - f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r15^2*r24^2*s^4 - f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r14^2*r25^2*s^4 - f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r15^2*r23^2*s^4 - f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r13^2*r25^2*s^4 - f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r14^2*r23^2*s^4 - f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r13^2*r24^2*s^4 - f1^2*f2^2*f3^2*f4^2*f5^2*q^6*s^4 + f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r25^2*s^4 + f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r24^2*s^4 + f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r23^2*s^4 + f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r15^2*s^4 + f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r14^2*s^4 + f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r13^2*s^4, z1, 1)

root(z1^5 - z1^4*(f3^2*q^2 + f4^2*q^2 + f5^2*q^2 + f1^2*s^2 + f2^2*s^2) + z1^3*(f3^2*f4^2*q^4 + f3^2*f5^2*q^4 + f4^2*f5^2*q^4 + f1^2*f2^2*s^4 + f1^2*f3^2*q^2*s^2 + f1^2*f4^2*q^2*s^2 + f2^2*f3^2*q^2*s^2 + f1^2*f5^2*q^2*s^2 + f2^2*f4^2*q^2*s^2 + f2^2*f5^2*q^2*s^2 - f1^2*f3^2*q^2*r13^2*s^2 - f1^2*f4^2*q^2*r14^2*s^2 - f1^2*f5^2*q^2*r15^2*s^2 - f2^2*f3^2*q^2*r23^2*s^2 - f2^2*f4^2*q^2*r24^2*s^2 - f2^2*f5^2*q^2*r25^2*s^2) + z1^2*(- f3^2*f4^2*f5^2*q^6 - f1^2*f2^2*f3^2*q^2*s^4 - f1^2*f2^2*f4^2*q^2*s^4 - f1^2*f2^2*f5^2*q^2*s^4 - f1^2*f3^2*f4^2*q^4*s^2 - f1^2*f3^2*f5^2*q^4*s^2 - f2^2*f3^2*f4^2*q^4*s^2 - f1^2*f4^2*f5^2*q^4*s^2 - f2^2*f3^2*f5^2*q^4*s^2 - f2^2*f4^2*f5^2*q^4*s^2 + f1^2*f2^2*f3^2*q^2*r13^2*s^4 + f1^2*f2^2*f4^2*q^2*r14^2*s^4 + f1^2*f3^2*f4^2*q^4*r13^2*s^2 + f1^2*f3^2*f4^2*q^4*r14^2*s^2 + f1^2*f3^2*f5^2*q^4*r13^2*s^2 + f1^2*f2^2*f5^2*q^2*r15^2*s^4 + f1^2*f3^2*f5^2*q^4*r15^2*s^2 + f1^2*f4^2*f5^2*q^4*r14^2*s^2 + f1^2*f4^2*f5^2*q^4*r15^2*s^2 + f1^2*f2^2*f3^2*q^2*r23^2*s^4 + f1^2*f2^2*f4^2*q^2*r24^2*s^4 + f2^2*f3^2*f4^2*q^4*r23^2*s^2 + f1^2*f2^2*f5^2*q^2*r25^2*s^4 + f2^2*f3^2*f4^2*q^4*r24^2*s^2 + f2^2*f3^2*f5^2*q^4*r23^2*s^2 + f2^2*f3^2*f5^2*q^4*r25^2*s^2 + f2^2*f4^2*f5^2*q^4*r24^2*s^2 + f2^2*f4^2*f5^2*q^4*r25^2*s^2) - z1*(- f1^2*f2^2*f3^2*f4^2*q^4*s^4 - f1^2*f2^2*f3^2*f5^2*q^4*s^4 - f1^2*f2^2*f4^2*f5^2*q^4*s^4 - f1^2*f3^2*f4^2*f5^2*q^6*s^2 - f2^2*f3^2*f4^2*f5^2*q^6*s^2 + f1^2*f2^2*f3^2*f4^2*q^4*r13^2*s^4 + f1^2*f2^2*f3^2*f4^2*q^4*r14^2*s^4 + f1^2*f2^2*f3^2*f5^2*q^4*r13^2*s^4 + f1^2*f2^2*f3^2*f5^2*q^4*r15^2*s^4 + f1^2*f2^2*f4^2*f5^2*q^4*r14^2*s^4 + f1^2*f3^2*f4^2*f5^2*q^6*r13^2*s^2 + f1^2*f2^2*f4^2*f5^2*q^4*r15^2*s^4 + f1^2*f3^2*f4^2*f5^2*q^6*r14^2*s^2 + f1^2*f3^2*f4^2*f5^2*q^6*r15^2*s^2 + f1^2*f2^2*f3^2*f4^2*q^4*r23^2*s^4 + f1^2*f2^2*f3^2*f4^2*q^4*r24^2*s^4 + f1^2*f2^2*f3^2*f5^2*q^4*r23^2*s^4 + f1^2*f2^2*f3^2*f5^2*q^4*r25^2*s^4 + f1^2*f2^2*f4^2*f5^2*q^4*r24^2*s^4 + f1^2*f2^2*f4^2*f5^2*q^4*r25^2*s^4 + f2^2*f3^2*f4^2*f5^2*q^6*r23^2*s^2 + f2^2*f3^2*f4^2*f5^2*q^6*r24^2*s^2 + f2^2*f3^2*f4^2*f5^2*q^6*r25^2*s^2 - f1^2*f2^2*f3^2*f4^2*q^4*r13^2*r24^2*s^4 - f1^2*f2^2*f3^2*f4^2*q^4*r14^2*r23^2*s^4 - f1^2*f2^2*f3^2*f5^2*q^4*r13^2*r25^2*s^4 - f1^2*f2^2*f3^2*f5^2*q^4*r15^2*r23^2*s^4 - f1^2*f2^2*f4^2*f5^2*q^4*r14^2*r25^2*s^4 - f1^2*f2^2*f4^2*f5^2*q^4*r15^2*r24^2*s^4 + 2*f1^2*f2^2*f3^2*f4^2*q^4*r13*r14*r23*r24*s^4 + 2*f1^2*f2^2*f3^2*f5^2*q^4*r13*r15*r23*r25*s^4 + 2*f1^2*f2^2*f4^2*f5^2*q^4*r14*r15*r24*r25*s^4) + 2*f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r14*r15*r24*r25*s^4 + 2*f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r13*r15*r23*r25*s^4 + 2*f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r13*r14*r23*r24*s^4 - f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r15^2*r24^2*s^4 - f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r14^2*r25^2*s^4 - f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r15^2*r23^2*s^4 - f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r13^2*r25^2*s^4 - f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r14^2*r23^2*s^4 - f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r13^2*r24^2*s^4 - f1^2*f2^2*f3^2*f4^2*f5^2*q^6*s^4 + f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r25^2*s^4 + f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r24^2*s^4 + f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r23^2*s^4 + f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r15^2*s^4 + f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r14^2*s^4 + f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r13^2*s^4, z1, 2)

root(z1^5 - z1^4*(f3^2*q^2 + f4^2*q^2 + f5^2*q^2 + f1^2*s^2 + f2^2*s^2) + z1^3*(f3^2*f4^2*q^4 + f3^2*f5^2*q^4 + f4^2*f5^2*q^4 + f1^2*f2^2*s^4 + f1^2*f3^2*q^2*s^2 + f1^2*f4^2*q^2*s^2 + f2^2*f3^2*q^2*s^2 + f1^2*f5^2*q^2*s^2 + f2^2*f4^2*q^2*s^2 + f2^2*f5^2*q^2*s^2 - f1^2*f3^2*q^2*r13^2*s^2 - f1^2*f4^2*q^2*r14^2*s^2 - f1^2*f5^2*q^2*r15^2*s^2 - f2^2*f3^2*q^2*r23^2*s^2 - f2^2*f4^2*q^2*r24^2*s^2 - f2^2*f5^2*q^2*r25^2*s^2) + z1^2*(- f3^2*f4^2*f5^2*q^6 - f1^2*f2^2*f3^2*q^2*s^4 - f1^2*f2^2*f4^2*q^2*s^4 - f1^2*f2^2*f5^2*q^2*s^4 - f1^2*f3^2*f4^2*q^4*s^2 - f1^2*f3^2*f5^2*q^4*s^2 - f2^2*f3^2*f4^2*q^4*s^2 - f1^2*f4^2*f5^2*q^4*s^2 - f2^2*f3^2*f5^2*q^4*s^2 - f2^2*f4^2*f5^2*q^4*s^2 + f1^2*f2^2*f3^2*q^2*r13^2*s^4 + f1^2*f2^2*f4^2*q^2*r14^2*s^4 + f1^2*f3^2*f4^2*q^4*r13^2*s^2 + f1^2*f3^2*f4^2*q^4*r14^2*s^2 + f1^2*f3^2*f5^2*q^4*r13^2*s^2 + f1^2*f2^2*f5^2*q^2*r15^2*s^4 + f1^2*f3^2*f5^2*q^4*r15^2*s^2 + f1^2*f4^2*f5^2*q^4*r14^2*s^2 + f1^2*f4^2*f5^2*q^4*r15^2*s^2 + f1^2*f2^2*f3^2*q^2*r23^2*s^4 + f1^2*f2^2*f4^2*q^2*r24^2*s^4 + f2^2*f3^2*f4^2*q^4*r23^2*s^2 + f1^2*f2^2*f5^2*q^2*r25^2*s^4 + f2^2*f3^2*f4^2*q^4*r24^2*s^2 + f2^2*f3^2*f5^2*q^4*r23^2*s^2 + f2^2*f3^2*f5^2*q^4*r25^2*s^2 + f2^2*f4^2*f5^2*q^4*r24^2*s^2 + f2^2*f4^2*f5^2*q^4*r25^2*s^2) - z1*(- f1^2*f2^2*f3^2*f4^2*q^4*s^4 - f1^2*f2^2*f3^2*f5^2*q^4*s^4 - f1^2*f2^2*f4^2*f5^2*q^4*s^4 - f1^2*f3^2*f4^2*f5^2*q^6*s^2 - f2^2*f3^2*f4^2*f5^2*q^6*s^2 + f1^2*f2^2*f3^2*f4^2*q^4*r13^2*s^4 + f1^2*f2^2*f3^2*f4^2*q^4*r14^2*s^4 + f1^2*f2^2*f3^2*f5^2*q^4*r13^2*s^4 + f1^2*f2^2*f3^2*f5^2*q^4*r15^2*s^4 + f1^2*f2^2*f4^2*f5^2*q^4*r14^2*s^4 + f1^2*f3^2*f4^2*f5^2*q^6*r13^2*s^2 + f1^2*f2^2*f4^2*f5^2*q^4*r15^2*s^4 + f1^2*f3^2*f4^2*f5^2*q^6*r14^2*s^2 + f1^2*f3^2*f4^2*f5^2*q^6*r15^2*s^2 + f1^2*f2^2*f3^2*f4^2*q^4*r23^2*s^4 + f1^2*f2^2*f3^2*f4^2*q^4*r24^2*s^4 + f1^2*f2^2*f3^2*f5^2*q^4*r23^2*s^4 + f1^2*f2^2*f3^2*f5^2*q^4*r25^2*s^4 + f1^2*f2^2*f4^2*f5^2*q^4*r24^2*s^4 + f1^2*f2^2*f4^2*f5^2*q^4*r25^2*s^4 + f2^2*f3^2*f4^2*f5^2*q^6*r23^2*s^2 + f2^2*f3^2*f4^2*f5^2*q^6*r24^2*s^2 + f2^2*f3^2*f4^2*f5^2*q^6*r25^2*s^2 - f1^2*f2^2*f3^2*f4^2*q^4*r13^2*r24^2*s^4 - f1^2*f2^2*f3^2*f4^2*q^4*r14^2*r23^2*s^4 - f1^2*f2^2*f3^2*f5^2*q^4*r13^2*r25^2*s^4 - f1^2*f2^2*f3^2*f5^2*q^4*r15^2*r23^2*s^4 - f1^2*f2^2*f4^2*f5^2*q^4*r14^2*r25^2*s^4 - f1^2*f2^2*f4^2*f5^2*q^4*r15^2*r24^2*s^4 + 2*f1^2*f2^2*f3^2*f4^2*q^4*r13*r14*r23*r24*s^4 + 2*f1^2*f2^2*f3^2*f5^2*q^4*r13*r15*r23*r25*s^4 + 2*f1^2*f2^2*f4^2*f5^2*q^4*r14*r15*r24*r25*s^4) + 2*f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r14*r15*r24*r25*s^4 + 2*f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r13*r15*r23*r25*s^4 + 2*f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r13*r14*r23*r24*s^4 - f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r15^2*r24^2*s^4 - f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r14^2*r25^2*s^4 - f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r15^2*r23^2*s^4 - f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r13^2*r25^2*s^4 - f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r14^2*r23^2*s^4 - f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r13^2*r24^2*s^4 - f1^2*f2^2*f3^2*f4^2*f5^2*q^6*s^4 + f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r25^2*s^4 + f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r24^2*s^4 + f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r23^2*s^4 + f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r15^2*s^4 + f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r14^2*s^4 + f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r13^2*s^4, z1, 3)

root(z1^5 - z1^4*(f3^2*q^2 + f4^2*q^2 + f5^2*q^2 + f1^2*s^2 + f2^2*s^2) + z1^3*(f3^2*f4^2*q^4 + f3^2*f5^2*q^4 + f4^2*f5^2*q^4 + f1^2*f2^2*s^4 + f1^2*f3^2*q^2*s^2 + f1^2*f4^2*q^2*s^2 + f2^2*f3^2*q^2*s^2 + f1^2*f5^2*q^2*s^2 + f2^2*f4^2*q^2*s^2 + f2^2*f5^2*q^2*s^2 - f1^2*f3^2*q^2*r13^2*s^2 - f1^2*f4^2*q^2*r14^2*s^2 - f1^2*f5^2*q^2*r15^2*s^2 - f2^2*f3^2*q^2*r23^2*s^2 - f2^2*f4^2*q^2*r24^2*s^2 - f2^2*f5^2*q^2*r25^2*s^2) + z1^2*(- f3^2*f4^2*f5^2*q^6 - f1^2*f2^2*f3^2*q^2*s^4 - f1^2*f2^2*f4^2*q^2*s^4 - f1^2*f2^2*f5^2*q^2*s^4 - f1^2*f3^2*f4^2*q^4*s^2 - f1^2*f3^2*f5^2*q^4*s^2 - f2^2*f3^2*f4^2*q^4*s^2 - f1^2*f4^2*f5^2*q^4*s^2 - f2^2*f3^2*f5^2*q^4*s^2 - f2^2*f4^2*f5^2*q^4*s^2 + f1^2*f2^2*f3^2*q^2*r13^2*s^4 + f1^2*f2^2*f4^2*q^2*r14^2*s^4 + f1^2*f3^2*f4^2*q^4*r13^2*s^2 + f1^2*f3^2*f4^2*q^4*r14^2*s^2 + f1^2*f3^2*f5^2*q^4*r13^2*s^2 + f1^2*f2^2*f5^2*q^2*r15^2*s^4 + f1^2*f3^2*f5^2*q^4*r15^2*s^2 + f1^2*f4^2*f5^2*q^4*r14^2*s^2 + f1^2*f4^2*f5^2*q^4*r15^2*s^2 + f1^2*f2^2*f3^2*q^2*r23^2*s^4 + f1^2*f2^2*f4^2*q^2*r24^2*s^4 + f2^2*f3^2*f4^2*q^4*r23^2*s^2 + f1^2*f2^2*f5^2*q^2*r25^2*s^4 + f2^2*f3^2*f4^2*q^4*r24^2*s^2 + f2^2*f3^2*f5^2*q^4*r23^2*s^2 + f2^2*f3^2*f5^2*q^4*r25^2*s^2 + f2^2*f4^2*f5^2*q^4*r24^2*s^2 + f2^2*f4^2*f5^2*q^4*r25^2*s^2) - z1*(- f1^2*f2^2*f3^2*f4^2*q^4*s^4 - f1^2*f2^2*f3^2*f5^2*q^4*s^4 - f1^2*f2^2*f4^2*f5^2*q^4*s^4 - f1^2*f3^2*f4^2*f5^2*q^6*s^2 - f2^2*f3^2*f4^2*f5^2*q^6*s^2 + f1^2*f2^2*f3^2*f4^2*q^4*r13^2*s^4 + f1^2*f2^2*f3^2*f4^2*q^4*r14^2*s^4 + f1^2*f2^2*f3^2*f5^2*q^4*r13^2*s^4 + f1^2*f2^2*f3^2*f5^2*q^4*r15^2*s^4 + f1^2*f2^2*f4^2*f5^2*q^4*r14^2*s^4 + f1^2*f3^2*f4^2*f5^2*q^6*r13^2*s^2 + f1^2*f2^2*f4^2*f5^2*q^4*r15^2*s^4 + f1^2*f3^2*f4^2*f5^2*q^6*r14^2*s^2 + f1^2*f3^2*f4^2*f5^2*q^6*r15^2*s^2 + f1^2*f2^2*f3^2*f4^2*q^4*r23^2*s^4 + f1^2*f2^2*f3^2*f4^2*q^4*r24^2*s^4 + f1^2*f2^2*f3^2*f5^2*q^4*r23^2*s^4 + f1^2*f2^2*f3^2*f5^2*q^4*r25^2*s^4 + f1^2*f2^2*f4^2*f5^2*q^4*r24^2*s^4 + f1^2*f2^2*f4^2*f5^2*q^4*r25^2*s^4 + f2^2*f3^2*f4^2*f5^2*q^6*r23^2*s^2 + f2^2*f3^2*f4^2*f5^2*q^6*r24^2*s^2 + f2^2*f3^2*f4^2*f5^2*q^6*r25^2*s^2 - f1^2*f2^2*f3^2*f4^2*q^4*r13^2*r24^2*s^4 - f1^2*f2^2*f3^2*f4^2*q^4*r14^2*r23^2*s^4 - f1^2*f2^2*f3^2*f5^2*q^4*r13^2*r25^2*s^4 - f1^2*f2^2*f3^2*f5^2*q^4*r15^2*r23^2*s^4 - f1^2*f2^2*f4^2*f5^2*q^4*r14^2*r25^2*s^4 - f1^2*f2^2*f4^2*f5^2*q^4*r15^2*r24^2*s^4 + 2*f1^2*f2^2*f3^2*f4^2*q^4*r13*r14*r23*r24*s^4 + 2*f1^2*f2^2*f3^2*f5^2*q^4*r13*r15*r23*r25*s^4 + 2*f1^2*f2^2*f4^2*f5^2*q^4*r14*r15*r24*r25*s^4) + 2*f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r14*r15*r24*r25*s^4 + 2*f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r13*r15*r23*r25*s^4 + 2*f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r13*r14*r23*r24*s^4 - f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r15^2*r24^2*s^4 - f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r14^2*r25^2*s^4 - f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r15^2*r23^2*s^4 - f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r13^2*r25^2*s^4 - f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r14^2*r23^2*s^4 - f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r13^2*r24^2*s^4 - f1^2*f2^2*f3^2*f4^2*f5^2*q^6*s^4 + f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r25^2*s^4 + f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r24^2*s^4 + f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r23^2*s^4 + f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r15^2*s^4 + f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r14^2*s^4 + f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r13^2*s^4, z1, 4)

root(z1^5 - z1^4*(f3^2*q^2 + f4^2*q^2 + f5^2*q^2 + f1^2*s^2 + f2^2*s^2) + z1^3*(f3^2*f4^2*q^4 + f3^2*f5^2*q^4 + f4^2*f5^2*q^4 + f1^2*f2^2*s^4 + f1^2*f3^2*q^2*s^2 + f1^2*f4^2*q^2*s^2 + f2^2*f3^2*q^2*s^2 + f1^2*f5^2*q^2*s^2 + f2^2*f4^2*q^2*s^2 + f2^2*f5^2*q^2*s^2 - f1^2*f3^2*q^2*r13^2*s^2 - f1^2*f4^2*q^2*r14^2*s^2 - f1^2*f5^2*q^2*r15^2*s^2 - f2^2*f3^2*q^2*r23^2*s^2 - f2^2*f4^2*q^2*r24^2*s^2 - f2^2*f5^2*q^2*r25^2*s^2) + z1^2*(- f3^2*f4^2*f5^2*q^6 - f1^2*f2^2*f3^2*q^2*s^4 - f1^2*f2^2*f4^2*q^2*s^4 - f1^2*f2^2*f5^2*q^2*s^4 - f1^2*f3^2*f4^2*q^4*s^2 - f1^2*f3^2*f5^2*q^4*s^2 - f2^2*f3^2*f4^2*q^4*s^2 - f1^2*f4^2*f5^2*q^4*s^2 - f2^2*f3^2*f5^2*q^4*s^2 - f2^2*f4^2*f5^2*q^4*s^2 + f1^2*f2^2*f3^2*q^2*r13^2*s^4 + f1^2*f2^2*f4^2*q^2*r14^2*s^4 + f1^2*f3^2*f4^2*q^4*r13^2*s^2 + f1^2*f3^2*f4^2*q^4*r14^2*s^2 + f1^2*f3^2*f5^2*q^4*r13^2*s^2 + f1^2*f2^2*f5^2*q^2*r15^2*s^4 + f1^2*f3^2*f5^2*q^4*r15^2*s^2 + f1^2*f4^2*f5^2*q^4*r14^2*s^2 + f1^2*f4^2*f5^2*q^4*r15^2*s^2 + f1^2*f2^2*f3^2*q^2*r23^2*s^4 + f1^2*f2^2*f4^2*q^2*r24^2*s^4 + f2^2*f3^2*f4^2*q^4*r23^2*s^2 + f1^2*f2^2*f5^2*q^2*r25^2*s^4 + f2^2*f3^2*f4^2*q^4*r24^2*s^2 + f2^2*f3^2*f5^2*q^4*r23^2*s^2 + f2^2*f3^2*f5^2*q^4*r25^2*s^2 + f2^2*f4^2*f5^2*q^4*r24^2*s^2 + f2^2*f4^2*f5^2*q^4*r25^2*s^2) - z1*(- f1^2*f2^2*f3^2*f4^2*q^4*s^4 - f1^2*f2^2*f3^2*f5^2*q^4*s^4 - f1^2*f2^2*f4^2*f5^2*q^4*s^4 - f1^2*f3^2*f4^2*f5^2*q^6*s^2 - f2^2*f3^2*f4^2*f5^2*q^6*s^2 + f1^2*f2^2*f3^2*f4^2*q^4*r13^2*s^4 + f1^2*f2^2*f3^2*f4^2*q^4*r14^2*s^4 + f1^2*f2^2*f3^2*f5^2*q^4*r13^2*s^4 + f1^2*f2^2*f3^2*f5^2*q^4*r15^2*s^4 + f1^2*f2^2*f4^2*f5^2*q^4*r14^2*s^4 + f1^2*f3^2*f4^2*f5^2*q^6*r13^2*s^2 + f1^2*f2^2*f4^2*f5^2*q^4*r15^2*s^4 + f1^2*f3^2*f4^2*f5^2*q^6*r14^2*s^2 + f1^2*f3^2*f4^2*f5^2*q^6*r15^2*s^2 + f1^2*f2^2*f3^2*f4^2*q^4*r23^2*s^4 + f1^2*f2^2*f3^2*f4^2*q^4*r24^2*s^4 + f1^2*f2^2*f3^2*f5^2*q^4*r23^2*s^4 + f1^2*f2^2*f3^2*f5^2*q^4*r25^2*s^4 + f1^2*f2^2*f4^2*f5^2*q^4*r24^2*s^4 + f1^2*f2^2*f4^2*f5^2*q^4*r25^2*s^4 + f2^2*f3^2*f4^2*f5^2*q^6*r23^2*s^2 + f2^2*f3^2*f4^2*f5^2*q^6*r24^2*s^2 + f2^2*f3^2*f4^2*f5^2*q^6*r25^2*s^2 - f1^2*f2^2*f3^2*f4^2*q^4*r13^2*r24^2*s^4 - f1^2*f2^2*f3^2*f4^2*q^4*r14^2*r23^2*s^4 - f1^2*f2^2*f3^2*f5^2*q^4*r13^2*r25^2*s^4 - f1^2*f2^2*f3^2*f5^2*q^4*r15^2*r23^2*s^4 - f1^2*f2^2*f4^2*f5^2*q^4*r14^2*r25^2*s^4 - f1^2*f2^2*f4^2*f5^2*q^4*r15^2*r24^2*s^4 + 2*f1^2*f2^2*f3^2*f4^2*q^4*r13*r14*r23*r24*s^4 + 2*f1^2*f2^2*f3^2*f5^2*q^4*r13*r15*r23*r25*s^4 + 2*f1^2*f2^2*f4^2*f5^2*q^4*r14*r15*r24*r25*s^4) + 2*f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r14*r15*r24*r25*s^4 + 2*f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r13*r15*r23*r25*s^4 + 2*f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r13*r14*r23*r24*s^4 - f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r15^2*r24^2*s^4 - f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r14^2*r25^2*s^4 - f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r15^2*r23^2*s^4 - f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r13^2*r25^2*s^4 - f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r14^2*r23^2*s^4 - f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r13^2*r24^2*s^4 - f1^2*f2^2*f3^2*f4^2*f5^2*q^6*s^4 + f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r25^2*s^4 + f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r24^2*s^4 + f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r23^2*s^4 + f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r15^2*s^4 + f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r14^2*s^4 + f1^2*f2^2*f3^2*f4^2*f5^2*q^6*r13^2*s^4, z1, 5)