k7=7*r, k6=-21*r^2+7*b*h*m+7*c*g*m+7*d*e*m, k5=35*r^3+m*(e*(14*b*c-35*d*r)+g*(7*b^2-35*c*r)-35*b*h*r+7*b*d^2+7*c^2*d)+(7*c*h^2+(14*d*g+7*e^2)*h+7*e*g^2)*m^2, k4=-35*r^4+m*(e*(70*d*r^2-56*b*c*r+7*b^3)+g*(70*c*r^2-28*b^2*r)+70*b*h*r^2+d*(21*b^2*c-28*c^2*r)-28*b*d^2*r+7*b*c^3)+m^2*(h*(g*(-56*d*r-7*b*c)- 28*e^2*r+(21*c^2-7*b*d)*e+21*c*d^2)+g^2*(-28*e*r+21*b*d-14*c^2)+h^2*(-28*c*r-14*b^2)+(21*b*e^2-7*c*d*e+7*d^3)*g+7*c*e^3-14*d^2*e^2)+ (7*d*h^3+21*e*g*h^2+7*g^3*h)*m^3, k3=21*r^5+m*(e*(-70*d*r^3+84*b*c*r^2-21*b^3*r)+g*(42*b^2*r^2-70*c*r^3)-70*b*h*r^3+d*(42*c^2*r^2-63*b^2*c*r+7*b^4)+42*b*d^2*r^2-21*b*c^3*r+14*b^3*c^2)+ m^2*(h*(g*(84*d*r^2+21*b*c*r-21*b^3)+42*e^2*r^2+e*(21*b*d*r-63*c^2*r-14*b^2*c)+d^2*(-63*c*r-7*b^2)+35*b*c^2*d+7*c^4)+g^2*(42*e*r^2-63*b*d*r+ 42*c^2*r-7*b^2*c)+h^2*(42*c*r^2+42*b^2*r)+g*(e*(d*(21*c*r+35*b^2)-14*b*c^2)-63*b*e^2*r-21*d^3*r+35*b*c*d^2-21*c^3*d)+e^2*(42*d^2*r-14*b*c*d+ 14*c^3)+e^3*(14*b^2-21*c*r)+(-21*b*d^3-7*c^2*d^2)*e+7*c*d^4)+m^3*(h^2*(g*(-63*e*r-14*b*d-7*c^2)-7*b*e^2+35*c*d*e+14*d^3)+h*(-21*g^3*r+(35*b*e- 14*c*d)*g^2+(35*c*e^2-14*d^2*e)*g-21*d*e^3)+h^3*(-21*d*r-21*b*c)+7*b*g^4+(14*d^2-21*c*e)*g^3-7*d*e^2*g^2+7*e^4*g)+(7*e*h^4+14*g^2*h^3)*m^4, k2=-7*r^6+m*(e*(35*d*r^4-56*b*c*r^3+21*b^3*r^2)+g*(35*c*r^4-28*b^2*r^3)+35*b*h*r^4+d*(-28*c^2*r^3+63*b^2*c*r^2-14*b^4*r)-28*b*d^2*r^3+21*b*c^3*r^2- 28*b^3*c^2*r+7*b^5*c)+m^2*(h*(g*(-56*d*r^3-21*b*c*r^2+42*b^3*r)-28*e^2*r^3+e*(-21*b*d*r^2+63*c^2*r^2+28*b^2*c*r-14*b^4)+d^2*(63*c*r^2+14*b^2*r)+ d*(-70*b*c^2*r-7*b^3*c)-14*c^4*r+21*b^2*c^3)+g^2*(-28*e*r^3+63*b*d*r^2-42*c^2*r^2+14*b^2*c*r-7*b^4)+h^2*(-28*c*r^3-42*b^2*r^2)+g*(e*(d*(-21*c*r^2- 70*b^2*r)+28*b*c^2*r-7*b^3*c)+63*b*e^2*r^2+21*d^3*r^2+d*(42*c^3*r+14*b^2*c^2)+d^2*(21*b^3-70*b*c*r)-14*b*c^4)+e^2*(-42*d^2*r^2+d*(28*b*c*r+21*b^3)- 28*c^3*r+7*b^2*c^2)+e^3*(21*c*r^2-28*b^2*r)+e*(d^2*(14*c^2*r-35*b^2*c)+42*b*d^3*r-7*b*c^3*d+7*c^5)+d^4*(-14*c*r-7*b^2)+21*b*c^2*d^3-7*c^4*d^2)+ m^3*(h^2*(g*(63*e*r^2+28*b*d*r+14*c^2*r-35*b^2*c)+e*(d*(14*b^2-70*c*r)-35*b*c^2)+14*b*e^2*r-28*d^3*r+14*b*c*d^2+21*c^3*d)+h*(21*g^3*r^2+g* (e*(28*d^2*r+105*b*c*d-7*c^3)+e^2*(14*b^2-70*c*r)-7*b*d^3-35*c^2*d^2)+g^2*(d*(28*c*r-35*b^2)-70*b*e*r+14*b*c^2)+e^3*(42*d*r-7*b*c)+(14*c^2*d- 35*b*d^2)*e^2-7*c*d^3*e+7*d^5)+h^3*(21*d*r^2+42*b*c*r+7*b^3)+g^2*(e^2*(14*d*r-35*b*c)+(14*b*d^2-35*c^2*d)*e+21*c*d^3)+g^3*(e*(42*c*r+21*b^2)- 28*d^2*r-7*b*c*d+7*c^3)+g*(-14*e^4*r+(21*c^2-7*b*d)*e^3+14*c*d^2*e^2-14*d^4*e)-14*b*g^4*r+7*b*e^5-14*c*d*e^4+7*d^3*e^3)+m^4*(h^3*(-28*g^2*r+ (-7*b*e-7*c*d)*g+21*c*e^2+21*d^2*e)+h^4*(-14*e*r-14*b*d-7*c^2)+(21*b*g^3+(14*c*e+7*d^2)*g^2-35*d*e^2*g-7*e^4)*h^2+(-14*c*g^4-7*d*e*g^3+ 21*e^3*g^2)*h+7*d*g^5-7*e^2*g^4)+7*g*h^5*m^5, k1=r^7+m*(e*(-7*d*r^5+14*b*c*r^4-7*b^3*r^3)+g*(7*b^2*r^4-7*c*r^5)-7*b*h*r^5+d*(7*c^2*r^4-21*b^2*c*r^3+7*b^4*r^2)+7*b*d^2*r^4-7*b*c^3*r^3+14*b^3*c^2*r^2- 7*b^5*c*r+b^7)+m^2*(h*(g*(14*d*r^4+7*b*c*r^3-21*b^3*r^2)+7*e^2*r^4+e*(7*b*d*r^3-21*c^2*r^3-14*b^2*c*r^2+14*b^4*r)+d^2*(-21*c*r^3-7*b^2*r^2)+ d*(35*b*c^2*r^2+7*b^3*c*r-7*b^5)+7*c^4*r^2-21*b^2*c^3*r+7*b^4*c^2)+g^2*(7*e*r^4-21*b*d*r^3+14*c^2*r^3-7*b^2*c*r^2+7*b^4*r)+h^2*(7*c*r^4+14*b^2*r^3)+ g*(e*(d*(7*c*r^3+35*b^2*r^2)-14*b*c^2*r^2+7*b^3*c*r-7*b^5)-21*b*e^2*r^3-7*d^3*r^3+d*(-21*c^3*r^2-14*b^2*c^2*r+14*b^4*c)+d^2*(35*b*c*r^2-21*b^3*r)+ 14*b*c^4*r-7*b^3*c^3)+e^2*(14*d^2*r^3+d*(-14*b*c*r^2-21*b^3*r)+14*c^3*r^2-7*b^2*c^2*r+7*b^4*c)+e^3*(14*b^2*r^2-7*c*r^3)+e*(d^2*(-7*c^2*r^2+ 35*b^2*c*r+7*b^4)-21*b*d^3*r^2+d*(7*b*c^3*r-21*b^3*c^2)-7*c^5*r+7*b^2*c^4)+d^4*(7*c*r^2+7*b^2*r)+d^2*(7*c^4*r+14*b^2*c^3)+d^3*(-21*b*c^2*r-7*b^3*c)- 7*b*c^5*d+c^7)+m^3*(h^2*(g*(-21*e*r^3-14*b*d*r^2-7*c^2*r^2+35*b^2*c*r+7*b^4)+e*(d*(35*c*r^2-14*b^2*r)+35*b*c^2*r-21*b^3*c)-7*b*e^2*r^2+14*d^3*r^2+ d*(-21*c^3*r-7*b^2*c^2)+d^2*(14*b^3-14*b*c*r)+7*b*c^4)+h*(-7*g^3*r^3+g*(e*(-14*d^2*r^2+d*(7*b^3-105*b*c*r)+7*c^3*r+35*b^2*c^2)+e^2*(35*c*r^2- 14*b^2*r)+d^2*(35*c^2*r-14*b^2*c)+7*b*d^3*r+7*b*c^3*d-7*c^5)+g^2*(d*(35*b^2*r-14*c*r^2)+35*b*e*r^2-14*b*c^2*r-21*b^3*c)+e^3*(-21*d*r^2+7*b*c*r- 7*b^3)+e^2*(d*(35*b^2*c-14*c^2*r)+35*b*d^2*r-21*b*c^3)+e*(d^3*(7*c*r-21*b^2)-14*b*c^2*d^2+14*c^4*d)-7*d^5*r+14*b*c*d^4-7*c^3*d^3)+h^3*(-7*d*r^3- 21*b*c*r^2-7*b^3*r)+g^2*(e^2*(-7*d*r^2+35*b*c*r+14*b^3)+e*(d*(35*c^2*r-14*b^2*c)-14*b*d^2*r-21*b*c^3)+d^3*(14*b^2-21*c*r)-7*b*c^2*d^2+7*c^4*d)+g^3* (e*(-21*c*r^2-21*b^2*r)+14*d^2*r^2+d*(7*b*c*r-7*b^3)-7*c^3*r+14*b^2*c^2)+g*(7*e^4*r^2+e^2*(d^2*(-14*c*r-7*b^2)+35*b*c^2*d+7*c^4)+e*(14*d^4*r+ 7*b*c*d^3-21*c^3*d^2)+e^3*(7*b*d*r-21*c^2*r-21*b^2*c)-7*b*d^5+7*c^2*d^4)+7*b*g^4*r^2+e^4*(d*(14*c*r+7*b^2)+7*b*c^2)+e^3*(-7*d^3*r-21*b*c*d^2- 7*c^3*d)-7*b*e^5*r+(7*b*d^4+14*c^2*d^3)*e^2-7*c*d^5*e+d^7)+m^4*(h^3*(14*g^2*r^2+g*(d*(7*c*r-21*b^2)+7*b*e*r-21*b*c^2)+e*(-21*d^2*r+7*b*c*d-7*c^3)+ e^2*(14*b^2-21*c*r)-7*b*d^3+14*c^2*d^2)+h^4*(7*e*r^2+14*b*d*r+7*c^2*r+7*b^2*c)+h^2*(g*(e^2*(35*d*r-14*b*c)+(35*b*d^2-14*c^2*d)*e-21*c*d^3)+ g^2*(e*(-14*c*r-7*b^2)-7*d^2*r+35*b*c*d+14*c^3)-21*b*g^3*r+7*e^4*r+(14*c^2-21*b*d)*e^3-7*c*d^2*e^2+7*d^4*e)+h*(g^3*(e*(7*d*r+7*b*c)-21*b*d^2- 21*c^2*d)+g^2*(-21*e^3*r+(-14*b*d-7*c^2)*e^2+35*c*d^2*e+7*d^4)+g^4*(14*c*r+7*b^2)+(14*b*e^4+7*c*d*e^3-21*d^3*e^2)*g-7*c*e^5+7*d^2*e^4)+g^4*(7*e^2*r+ (14*b*d+7*c^2)*e+7*c*d^2)+g^5*(-7*d*r-7*b*c)+(-7*b*e^3-21*c*d*e^2-7*d^3*e)*g^3+(7*c*e^4+14*d^2*e^3)*g^2-7*d*e^5*g+e^7)+m^5*(h^5*(-7*g*r-7*b*e- 7*c*d)+(7*b*g^2+(14*c*e+7*d^2)*g+7*d*e^2)*h^4+(-7*c*g^3-21*d*e*g^2-7*e^3*g)*h^3+(7*d*g^4+14*e^2*g^3)*h^2-7*e*g^5*h+g^7)+h^7*m^6, for the polynomial x^7-k7*x^6-k6*x^5-k5*x^4-k4*x^3-k3*x^2-k2*x-k1, roots: x1=r+b*m^(1/7)+c*m^(2/7)+d*m^(3/7)+e*m^(4/7)+g*m^(5/7)+h*m^(6/7), x2=r+%e^((2*%i*%pi)/7)*h*m^(6/7)+%e^((4*%i*%pi)/7)*g*m^(5/7)+%e^((6*%i*%pi)/7)*e*m^(4/7)+%e^(-(6*%i*%pi)/7)*d*m^(3/7)+%e^(-(4*%i*%pi)/7)*c*m^(2/7)+%e^(-(2*%i*%pi)/7)*b*m^(1/7), x3=r+%e^(-(2*%i*%pi)/7)*h*m^(6/7)+%e^(-(4*%i*%pi)/7)*g*m^(5/7)+%e^(-(6*%i*%pi)/7)*e*m^(4/7)+%e^((6*%i*%pi)/7)*d*m^(3/7)+%e^((4*%i*%pi)/7)*c*m^(2/7)+%e^((2*%i*%pi)/7)*b*m^(1/7), x4=r+%e^((4*%i*%pi)/7)*h*m^(6/7)+%e^(-(6*%i*%pi)/7)*g*m^(5/7)+%e^(-(2*%i*%pi)/7)*e*m^(4/7)+%e^((2*%i*%pi)/7)*d*m^(3/7)+%e^((6*%i*%pi)/7)*c*m^(2/7)+%e^(-(4*%i*%pi)/7)*b*m^(1/7), x5=r+%e^(-(4*%i*%pi)/7)*h*m^(6/7)+%e^((6*%i*%pi)/7)*g*m^(5/7)+%e^((2*%i*%pi)/7)*e*m^(4/7)+%e^(-(2*%i*%pi)/7)*d*m^(3/7)+%e^(-(6*%i*%pi)/7)*c*m^(2/7)+%e^((4*%i*%pi)/7)*b*m^(1/7), x6=r+%e^(-(6*%i*%pi)/7)*h*m^(6/7)+%e^((2*%i*%pi)/7)*g*m^(5/7)+%e^(-(4*%i*%pi)/7)*e*m^(4/7)+%e^((4*%i*%pi)/7)*d*m^(3/7)+%e^(-(2*%i*%pi)/7)*c*m^(2/7)+%e^((6*%i*%pi)/7)*b*m^(1/7), x7=r+%e^((6*%i*%pi)/7)*h*m^(6/7)+%e^(-(2*%i*%pi)/7)*g*m^(5/7)+%e^((4*%i*%pi)/7)*e*m^(4/7)+%e^(-(4*%i*%pi)/7)*d*m^(3/7)+%e^((2*%i*%pi)/7)*c*m^(2/7)+%e^(-(6*%i*%pi)/7)*b*m^(1/7).