func2p2.asv

来自「数值分析最常用的四十种算法」· ASV 代码 · 共 24 行

a=[20,2,3;1,8,1;2,-3,15];
[L,U]=lu(a)
b=eye(3,3)
n=3
X=[0,0,0]'
Y=[0,0,0]'
for i=1:n
   Y=inv(L)*b(:,i)      %计算结果，运用LU分解结果
   X=inv(U)*Y
   b(:,i)=X
end
inv(a)

%b =
%    0.0517   -0.0164   -0.0093
%   -0.0055    0.1237   -0.0072
%   -0.0080    0.0269    0.0665

%直接算逆矩阵的结果：