%例35-22：LU分解法求解下列方程組。
A=[3,1,-4,1;1,-5,0,7;0,2,1,-1;1,6,-1,-4]; %系數矩陣A
b=[13,-9,6,0]';  %矩陣B
[L,U]=lu(A);  %LU分解
x=U\(L\b)   % 求解
[L,U ,P]=lu(A); %  采用LU的第2種格式分解
x=U\(L\P*b)   % 求解