matlab含参方程组,MATLAB如何求解含参非线性方程组

x =

(- k1^2*l + k2*k1^2 + 2*k1*l^2 - l^3)/(k1*(- k1^2 + 2*k1*l + k2*k1 - l^2)) - (b^2*k3*root(z1^5 - (2*k1*k2*z1^3)/b^2 - (z1^2*(a*b^2*k1*k2 - a*b^2*k2*l))/(b^4*k3) + (z1*(k1^2*k2^2*k3 + a*b*k1*k2^2 + a*b*k1^2*k2 + a*b*k2*l^2 - 2*a*b*k1*k2*l))/(b^4*k3) - (a*k1^2*k2^2 - a*k1*k2^2*l)/(b^4*k3), z1, 1)^3)/(a*k1^2 - 2*a*k1*l - a*k2*k1 + a*l^2) - (root(z1^5 - (2*k1*k2*z1^3)/b^2 - (z1^2*(a*b^2*k1*k2 - a*b^2*k2*l))/(b^4*k3) + (z1*(k1^2*k2^2*k3 + a*b*k1*k2^2 + a*b*k1^2*k2 + a*b*k2*l^2 - 2*a*b*k1*k2*l))/(b^4*k3) - (a*k1^2*k2^2 - a*k1*k2^2*l)/(b^4*k3), z1, 1)*(a*b*k1^2 + a*b*l^2 + k1^2*k2*k3 - 2*a*b*k1*l))/(a*k1*(- k1^2 + 2*k1*l + k2*k1 - l^2)) + (b^3*k3*(k1 - l)*root(z1^5 - (2*k1*k2*z1^3)/b^2 - (z1^2*(a*b^2*k1*k2 - a*b^2*k2*l))/(b^4*k3) + (z1*(k1^2*k2^2*k3 + a*b*k1*k2^2 + a*b*k1^2*k2 + a*b*k2*l^2 - 2*a*b*k1*k2*l))/(b^4*k3) - (a*k1^2*k2^2 - a*k1*k2^2*l)/(b^4*k3), z1, 1)^4)/(k1*(- a*k1^2*k2 + a*k1*k2^2 + 2*a*k1*k2*l - a*k2*l^2)) - (b*k3*(k1 - l)*root(z1^5 - (2*k1*k2*z1^3)/b^2 - (z1^2*(a*b^2*k1*k2 - a*b^2*k2*l))/(b^4*k3) + (z1*(k1^2*k2^2*k3 + a*b*k1*k2^2 + a*b*k1^2*k2 + a*b*k2*l^2 - 2*a*b*k1*k2*l))/(b^4*k3) - (a*k1^2*k2^2 - a*k1*k2^2*l)/(b^4*k3), z1, 1)^2)/(a*(- k1^2 + 2*k1*l + k2*k1 - l^2))

(- k1^2*l + k2*k1^2 + 2*k1*l^2 - l^3)/(k1*(- k1^2 + 2*k1*l + k2*k1 - l^2)) - (b^2*k3*root(z1^5 - (2*k1*k2*z1^3)/b^2 - (z1^2*(a*b^2*k1*k2 - a*b^2*k2*l))/(b^4*k3) + (z1*(k1^2*k2^2*k3 + a*b*k1*k2^2 + a*b*k1^2*k2 + a*b*k2*l^2 - 2*a*b*k1*k2*l))/(b^4*k3) - (a*k1^2*k2^2 - a*k1*k2^2*l)/(b^4*k3), z1, 2)^3)/(a*k1^2 - 2*a*k1*l - a*k2*k1 + a*l^2) - (root(z1^5 - (2*k1*k2*z1^3)/b^2 - (z1^2*(a*b^2*k1*k2 - a*b^2*k2*l))/(b^4*k3) + (z1*(k1^2*k2^2*k3 + a*b*k1*k2^2 + a*b*k1^2*k2 + a*b*k2*l^2 - 2*a*b*k1*k2*l))/(b^4*k3) - (a*k1^2*k2^2 - a*k1*k2^2*l)/(b^4*k3), z1, 2)*(a*b*k1^2 + a*b*l^2 + k1^2*k2*k3 - 2*a*b*k1*l))/(a*k1*(- k1^2 + 2*k1*l + k2*k1 - l^2)) + (b^3*k3*(k1 - l)*root(z1^5 - (2*k1*k2*z1^3)/b^2 - (z1^2*(a*b^2*k1*k2 - a*b^2*k2*l))/(b^4*k3) + (z1*(k1^2*k2^2*k3 + a*b*k1*k2^2 + a*b*k1^2*k2 + a*b*k2*l^2 - 2*a*b*k1*k2*l))/(b^4*k3) - (a*k1^2*k2^2 - a*k1*k2^2*l)/(b^4*k3), z1, 2)^4)/(k1*(- a*k1^2*k2 + a*k1*k2^2 + 2*a*k1*k2*l - a*k2*l^2)) - (b*k3*(k1 - l)*root(z1^5 - (2*k1*k2*z1^3)/b^2 - (z1^2*(a*b^2*k1*k2 - a*b^2*k2*l))/(b^4*k3) + (z1*(k1^2*k2^2*k3 + a*b*k1*k2^2 + a*b*k1^2*k2 + a*b*k2*l^2 - 2*a*b*k1*k2*l))/(b^4*k3) - (a*k1^2*k2^2 - a*k1*k2^2*l)/(b^4*k3), z1, 2)^2)/(a*(- k1^2 + 2*k1*l + k2*k1 - l^2))

(- k1^2*l + k2*k1^2 + 2*k1*l^2 - l^3)/(k1*(- k1^2 + 2*k1*l + k2*k1 - l^2)) - (b^2*k3*root(z1^5 - (2*k1*k2*z1^3)/b^