% factorization_lu.m
% 矩陣的LU分解示例
A = [ 1    2    3
      4    5    6
      7    8    0 ];
disp('lu分解方法一');
[L1,U]=lu(A)
X=L1*U
disp('lu分解方法二');
[P,L2,U2]=lu(A)
PL2=P*L2
disp('lu分解求線性方程組');
b=[2 3 5]'
xb0=A\b
y1=L1\b;
xb1=U\y1