2、
程序
a=input('请输入线性方程组的系数矩阵a=');
b=input('请输入b=');
[m,n]=size(a);
if m~=n
error('此矩阵非方阵,运行截止')
end
for t=1:n-1
i=t;h=a(t,t);
for s=t+1:n
if abs(h)
h=a(s,t);i=s;
end
end
if i~=t
for j=1:n
c=a(i,j);a(i,j)=a(t,j);a(t,j)=c;
end
c=b(i);b(i)=b(t);b(t)=c;
end
for i=t+1:n
k=a(i,t)/a(t,t);
for j=1:n
a(i,j)=a(i,j)-k*a(t,j);
end
b(i)=b(i)-k*b(t);
end
end
fprintf('经%1d次消元后的增广矩阵为\n',(n-1));
disp([a,b])
x(n)=b(n)/a(n,n);
for i=n-1:-1:1
x(i)=b(i);
for j=i+1:n
x(i)=x(i)-x(j)*a(i,j);
end
x(i)=x(i)/a(i,i);
end
fprintf('线性方程组的解向量是');
x