Matlab方程求解问题.怎么解下面的方程!
function F =Komulainen(x,F_aq,F1_aq,F_org,V,A,C,D,a,a1,alfa,beta,x_pls,x_el)
% B = (beta*(35.15*x_pls+27.15*x(6)))
% Ba = ((beta*(35.15*x_pls+27.15*x(6)))*a)
% b = (x(6)-a*x_pls)
F = [(F_aq/V)*(x(2)-x(1))-K*(alfa*a*(-((beta*(35.15*x_pls+27.15*x(6)))*a-A+(x(6)-a*x(2)))-sqrt(((beta*(35.15*x_pls+27.15*x(6)))*a-A+(x(6)-a*x(2)))^2-4*a*(beta*(35.15*x_pls+27.15*x(6)))*(x(6)-a*x(2))))/(2*a)+alfa*(x(6)-a*x(2))+(1-alfa)*x(6)-x(4));
(F_org/V)*(x(6)-x(4))+K*(alfa*a*(-((beta*(35.15*x_pls+27.15*x(6)))*a-A+(x(6)-a*x(2)))-sqrt(((beta*(35.15*x_pls+27.15*x(6)))*a-A+(x(6)-a*x(2)))^2-4*a*(beta*(35.15*x_pls+27.15*x(6)))*(x(6)-a*x(2))))/(2*a)+alfa*(x(6)-a*x(2))+(1-alfa)*x(6)-x(4));
(F_aq/V)*(x_pls-x(2))-K*(alfa*a*(-((beta*(35.15*x_pls+27.15*x(6)))*a-A+(x(4)-a*x_pls))-sqrt(((beta*(35.15*x_pls+27.15*x(6)))*a-A+(x(4)-a*x_pls))^2-4*a*(beta*(35.15*x_pls+27.15*x(6)))*(x(4)-a*x_pls)))/(2*a)+alfa*(x(4)-a*x_pls)+(1-alfa)*x(4)-x(5));
(F_org/V)*(x(4)-x(5))+K*(alfa*a*(-((beta*(35.15*x_pls+27.15*x(6)))*a-A+(x(4)-a*x_pls))-sqrt(((beta*(35.15*x_pls+27.15*x(6)))*a-A+(x(4)-a*x_pls))^2-4*a*(beta*(35.15*x_pls+27.15*x(6)))*(x(4)-a*x_pls)))/(2*a)+alfa*(x(4)-a*x_pls)+(1-alfa)*x(4)-x(5));
(F1_aq/V)*(x_el-x(3))+K(alfa*a1*(((x(5)-a1*x_el)-D)/(C-a1))+alfa*(x(5)-a1*x_el)+(1-alfa)*x(5)-x(6));
(F_org/V)*(x(5)-x(6))-K(alfa*a1*(((x(5)-a1*x_el)-D)/(C-a1))+alfa*(x(5)-a1*x_el)+(1-alfa)*x(5)-x(6))];