func2p22.m

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

a=[20,2,3;1,8,1;2,-3,15];
[L,U]=lu(a)
b=eye(3,3)
for i=1:3
   Y=b(i)/L
   X=Y/U
   b(i)=X
end