以下是我的程序,参数b(1)b(2)b(3)b(4)一共18组,最后可以得到18组的解,但是现在只能显示出最后一组的解,想问问大家怎么把18组的结果都显示出来呢? syms x1 x2 x3 x4 filename1='m2.xls'; M=xlsread(filename1,'sheet2'); for i=1:18 b(1)=M(i,1);b(2)=M(i,2);b(3)=M(i,3);b(4)=M(i,10);b(5)=M(i,11); eq1=(3.55e-10)-x1.*x2.*exp(2.*b(2).*b(5)+2.*b(1).*b(4)).*exp(-0.5547.*((b(4)+b(5)+4.*x4+x2+x3)./2).^0.5./(1+1.2.*((b(4)+b(5)+4.*x4+x2+x3)./2).^0.5))./(b(4).*exp(-0.5547.*((b(4)+b(5)+4.*x4+x2+x3)./2).^0.5./(1+1.2.*((b(4)+b(5)+4.*x4+x2+x3)./2).^0.5)+2.*b(3).*b(5)+2.*b(1).*x1)); eq2=(4.60e-07)-x2.*b(5).*exp(-0.5547.*((b(4)+b(5)+4.*x4+x2+x3)./2).^0.5./(1+1.2.*((b(4)+b(5)+4.*x4+x2+x3)./2).^0.5)).*exp(-0.5547.*((b(4)+b(5)+4.*x4+x2+x3)./2).^0.5./(1+1.2.*((b(4)+b(5)+4.*x4+x2+x3)./2).^0.5))./(0.004941889.*1); eq3=(4.66e-11)-x2.*x4.*exp(-0.5547.*((b(4)+b(5)+4.*x4+x2+x3)./2).^0.5./(1+1.2.*((b(4)+b(5)+4.*x4+x2+x3)./2).^0.5)).*exp(-0.5547.*((b(4)+b(5)+4.*x4+x2+x3)./2).^0.5./(1+1.2.*((b(4)+b(5)+4.*x4+x2+x3)./2).^0.5))./(b(5).*exp(-0.5547.*((b(4)+b(5)+4.*x4+x2+x3)./2).^0.5./(1+1.2.*((b(4)+b(5)+4.*x4+x2+x3)./2).^0.5))); eq4=(1.0577e-14)-x2.*x3.*exp(-0.5547.*((b(4)+b(5)+4.*x4+x2+x3)./2).^0.5./(1+1.2.*((b(4)+b(5)+4.*x4+x2+x3)./2).^0.5)).*exp(-0.5547.*((b(4)+b(5)+4.*x4+x2+x3)./2).^0.5./(1+1.2.*((b(4)+b(5)+4.*x4+x2+x3)./2).^0.5)); end S=vpasolve(eq1,eq2,eq3,eq4,x1,x2,x3,x4); x1=double(S.x1) x2=double(S.x2) x3=double(S.x3) x4=double(S.x4) |