gauss.m

来自「这是线性议程组的直接解法」· M 代码 · 共 20 行

function [x,q]=Gauss(A,b)
B=[A,b];
if find(diag(A)==0)
    disp('不满足高斯迭代条件！');
    return
end
n=length(A);
for k=2:n
   for kk=k:n
    B(kk,:)=-B(kk,k-1)*B(k-1,:)/B(k-1,k-1)+B(kk,:)
   end
end
x(n,1)=B(n,n+1)/B(n,n);
for k=n-1:-1:1
    sx=B(k,k+1:n)*x(k+1:n,1);
    x(k,1)=(B(k,n+1)-sx)/B(k,k);
end
q=(n^3)/3+n^2-1/3*n;
end