backsub.m

来自「数值计算课本中第三章程序源代码」· M 代码 · 共 7 行

function X = backsub( A, 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