列主消元法（Matlab实现）

目录

1、概述

2、代码

3、结果

---

1、概述

在前面一章知识点我们已经梳理清晰啦，这里直接上代码和运行结果。

2、代码

function [x]=abc(A,b,flag) 
if nargin<3,flag=0;end 
n=length(b); 
for k=1:(n-1) 
	[ap,p]=max(abs(A(k:n,k))); 
	p=p+k-1; 
	if p>k 
		A([k p],:)=A([p k],:); 
		b([k p],:)=b([p k],:); 
	end 
	m=A(k+1:n,k)/A(k,k); 
	A(k+1:n,k+1:n)=A(k+1:n,k+1:n)-m*A(k,k+1:n); 
	b(k+1:n)=b(k+1:n)-m*b(k); 
	A(k+1:n,k)=zeros(n-k,1); 
	if flag~=0, Ab=[A,b],end 
end 
x=zeros(n,1); 
x(n)=b(n)/A(n,n); 
for k=n-1:-1:1 
	x(k)=(b(k)-A(k,k+1:n)*x(k+1:n))/A(k,k); 
end

3、结果