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
数值方法 回代程序(matlab)
最新推荐文章于 2024-09-21 10:03:18 发布