e242.m

来自「matlab算法集 matlab算法集」· M 代码 · 共 15 行

%-----------------------------------------------------------------------
% Example 2.4.2: LU Decomposition
%-----------------------------------------------------------------------
   clc                 
   clear
   fprintf ('Example 2.4.2: LU Decomposition\n')
   A = [0 1 -1 ; 2 -2 1 ; 1 2 0];
   b = [-4 ; 9 ; 0];
   show ('A',A)
   show ('b',b)
   [L,U,P,Delta] = lufac (A);
   [x,d] = lusub (L,U,P,b);
   show ('x',x)
%-----------------------------------------------------------------------