load west0479;                  %加載MATLAB 7.0自帶的一個稀疏矩陣
S = west0479;
[L,U,P] = luinc(S,'0');             %做不完全LU分解    
figure(1);                       %新開一個圖形顯示窗口
subplot(2,3,1);                   %把該窗口分成2*3個字窗口，并在第1個字窗口上畫圖
spy(S);                         %顯示稀疏矩陣S
title('S');                        %給圖添上標題
subplot(2,3,2);                   %在第2個字窗口上畫圖
spy(L);                         %顯示稀疏矩陣L
title('L');                        %給圖添上標題
subplot(2,3,3);                   %在第2個字窗口上畫圖
spy(U);                         %顯示稀疏矩陣U
title('U');                        %給圖添上標題
subplot(2,3,4);                   %在第2個字窗口上畫圖
spy(P*S);                       %顯示稀疏矩陣P*S
title('P*S');                      %給圖添上標題
subplot(2,3,5);                   %在第2個字窗口上畫圖
spy(L*U);                       %顯示稀疏矩陣L*U
title('L*U');                     %給圖添上標題