例1.52   矩陣的三角分解。
>> X=magic(4);
>> X(2,3)=7
X =
    	16     2     3    13
    	 5    11     7     8
    	 9     7     6    12
    	 4    14    15     1

>> [L,U]=lu(X)
% 產(chǎn)生一個(gè)下三角矩陣L和一個(gè)上三角矩陣U，并滿足X=LU。

L =
    	1.0000         0         0          0
    	0.3125    0.7685    1.0000        0
    	0.5625    0.4352    0.3864    1.0000
    	0.2500    1.0000         0         0
U =
   		16.0000    2.0000    3.0000   13.0000
  		     0     13.5000   14.2500   -2.2500
         	 0         0      -4.8889        5.6667
         	 0         0          0           3.4773

>> [L,U,P]=lu(X)
% 產(chǎn)生一個(gè)下三角矩陣L、一個(gè)上三角矩陣U和交換矩陣P，并滿足PX=LU。

L =
    	1.0000         0          0          0
    	0.2500    1.0000         0          0
    	0.3125    0.7685    1.0000         0
    	0.5625    0.4352    0.3864    1.0000
U =
   	16.0000    2.0000    3.0000   13.0000
      0     13.5000   14.2500   -2.2500
      0          0     -4.8889    5.6667
       0          0         0       3.4773
P =
     1     0     0     0
     0     0     0     1
     0     1     0     0
     0     0     1     0