backsub.m

数值计算课本中第三章程序源代码

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