bslashtx.m

Based on Matlab,Gauss Iteration Method

function x = bslashtx(A,b)
%   被交互实验trylu调用
%   BSLASHTX  Solve linear system (backslash)
%   x = bslashtx(A,b) solves A*x = b

[n,n] = size(A);

% Triangular factorization
[L,U,p] = lutx(A);

% Permutation and forward elimination
y = zeros(n,1);
for k = 1:n
   j = 1:k-1;
   y(k) = b(p(k)) - L(k,j)*y(j);
end

% Back substitution
x = zeros(n,1);
for k = n:-1:1
   j = k+1:n;
   x(k) = (y(k) - U(k,j)*x(j))/U(k,k);
end