-
.
高斯选列主元消元法
1. 程序:
clear
format rat
A=input('输入增广矩阵 A=')
[m,n]=size(A);
for i=1:(m-1)
numb=int2str(i);
disp(['第',numb,'次选列主元后的增广矩阵'])
temp=max(abs(A(i:m,i)));
[a,b]=find(abs(A(i:m,i))==temp);
tempo=A(a(1)+i-1,:);
A(a(1)+i-1,:)=A(i,:);
A(i,:)=tempo
disp(['第',numb,'次消元后的增广矩阵'])
for j=(i+1):m
A(j,:)=A(j,:)-A(i,:)*A(j,i)/A(i,i);
end
A
end
%回代过程
disp('回代求解')
- .word.zl.