Matlab实现——Forward-substitution

欢迎前往个人博客 驽马点滴 和视频空间 哔哩哔哩-《挨踢日志》

%Program 3.1.1 (Forward-substitution)

function X=forwardsub(A,B)

%Input  - A is an n x n lower-triangular nonsingular  matrix

%        - B is an n x 1 matrix

%Output - X is the solution to the linear system AX = B

%Find the dimension of B and initialize X

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

Untitled.m

 

A=[3,0,0,0;5,6,0,0;-4,7,-2,0;3,2,-1,4;];

B=[6;4;-7;20;];

X=forwardsub(A,B)

欢迎前往个人博客 驽马点滴 和视频空间 哔哩哔哩-《挨踢日志》