数值方法 回代程序（matlab）

function X=backsub(A,B)
%Input -A is an n x upper-triangular nonsigular matrix
%      -B is an n x 1 matrix
%Output -X is the solution to the linear system AX=B
%  Find hte dimension of B and hte initialize X
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