gauss.m

来自「不使用matlab内部调用程序的按照高斯消元法本身原理来编写。」· M 代码 · 共 23 行

function Gauss(A,B,n);
A=input('Please input an array A=');
B=input('Please input an array B=');
n=input('Please input number of x=');
x=zeros(1,n);
for k=1:n-1
    for i=k+1:n
        m=A(i,k)/A(k,k);
        for j=k+1:n
            A(i,j)=A(i,j)-m*A(k,j);
        end
        B(i)=B(i)-m*B(k);
    end
end
x(n)=B(n)/A(n,n);
for i=n-1:-1:1
    sum=0;
    for j=i+1:n
        sum=sum+A(i,j)*x(j);
    end
    x(i)=(B(i)-sum)/A(i,i);
end
x=x'