forsub.m

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

function X = forsub( A, B )
n = length( B );
X = zeros( n, 1 );
X(1) = B(1) / A(1,1);
for k = 2:1:n
    X(k) = ( B(k) - A(k,1:k-1) * X(1:k-1) ) / A(k,k);
end