小蜗牛的博客

guass消元A=input('输入系数矩阵A：');B=input('输入b向量（按行向量）：');n=length(B);x=zeros(n,1);c=zeros(1,n);d1=0;for i=1:n-1 max=abs(A(i,i)); m=i; for j=i+1:n if max<abs(A(j,j)) max=abs(A(j,j)); m=j; end end

一、MATLAB中实现1、代码function LU(a,n)m = zeros(n,n); for i = 1:n; m(i,i) = 1; end for j = 1 : n-1 if abs(a(j,j))<eps; error('zero pivot encountered'); end for i = j+1 : n mult = a(i,j)/a(j,j);