backsub.m

来自「余向阳博士Matlab教程的源代码。 包括：1.非线性方程的解法；2.线性方程」· M 代码 · 共 11 行

function X=backsub(A,b)
%  A是一个n阶上三角非奇异阵。
%  b是一个n维向量。
%  X是线性方程组AX=b的解。
n=length(b);
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