第一种：三角形方程组的解

a=[1 1 1;0 4 -1;2 -2 1];
b=[6;5;1];
n=3;
x=zeros(n,1);
for i=n:-1:1
%    s=0;
%    for j=i+1:n
%        s=s+a(i,j)*x(j);
%    end
    s=a(i,i+1:n)*x(i+1:n);
    x(i)=(b(i)-s)/a(i,i);
end
x

---

第二种：高斯消去法求解

A=[1 1 1;0 4 -1;2 -2 1];
b=[6;5;1];
n=3;
A(:,n+1)=b;
for k=1:n-1
    for i=k+1:n
        lik=A(i,k)/A(k,k);
        for j=k+1:n+1
            A(i,j)=A(i,j)-lik*A(k,j);
            
        end
        A(i,k)=0;
    end
end
A(i,j)