![](https://csdnimg.cn/release/download_crawler_static/86132197/bg8.jpg)
dif=RB-RA;
if dif>0
disp('此方程组无解');
return
end
if RA==RB
if RA==n
format long;
disp('此方程组有唯一解');
for p=1:n-1
for k=p+1:n
m=B(k,p)/B(p,p);
B(k,p:n+1)=B(k,p:n+1)-m*B(p,p:n+1);
end
end %顺序消元形成上三角矩阵
b=B(1:n,n+1);
A=B(1:n,1:n);
b(n)=b(n)/A(n,n);
for q=n-1:-1:1
b(q)=(b(q)-sum(A(q,q+1:n)*b(q+1:n)))/A(q,q);
end %回代求解
else
disp('此方程组有无数组解');
end
end
上机运行结果为
>>
A=[1.1348,3.8326,1.1651,3.4017;0.5301,1.7875,2.5330,1.5435;
3.4129,4.9317,8.7643,1.3142;1.2371,4.9998,10.6721,0.0147];
b=[9.5342;6.3941;18.4231;16.9237];
>> X=gaus(A,b)
此方程组有唯一解