zhongyao7.m

该程序用来解方程组的,非线性的和线性的,带有参数的方程组

syms x1 x2 x3 x4 x5 x6
a=2*sqrt(55)*sqrt(x1)/sqrt(9 - x1);
e=sqrt(x1*(10*(13+x1)/(9-x1)-1));
a1=x2;
b1=(-x2-x1*x2)/10;
c1=-x2*(9-x1)*e;
a2=((-9+x1)*(x1*(11+x1)/(-x1+9))^(1/2)*(-11*(1+x1)+90-x1*(22+x1))*x2^2)/(4*sqrt(11)*(-11979-2541*x1+11*x1^2+24*x1^3+x1^4));
b2=((x1*(11+x1)/(-x1+9))^(1/2)*(-54-3*x1+x1^2)*(-11*(x1+1)+90-x1*(22+x1))*x2^2)/(20*sqrt(11)*(-11979-2541*x1+11*x1^2+24*x1^3+x1^4));
c2=((1+x1)*(-99-(7+90-x1)*x1+2*x1^2)*x2^2)/(20*(-1089-132*x1+13*x1^2+x1^3));
d1=x3;
e1=(-a*x4+10*x3)/10;
g1=e*(-a*x4*x1+20*x3*x1+a^2*x3+20*x4*a)/10;
f1='e+a1+a2-x5';
f2='a1+a2-(m+d1)*e^(-1)';
f3='e+b1+b2-x5';
f4='b1+b2-2*e1*e^(-1)';
f5='e+c1+c2-x6';
f6='c1+c2-2*g1*e^(-1)';
[x1,x2,x3,x4,x5,x6]=fsolve(f1,f1,f1,f1,f1,f1);
solution=[x1,x2,x3,x4,x5,x6]