lusovlve.m

来自「数学工具： LU分解法解线性方程组」· M 代码 · 共 17 行

function x=LUsovlve(l,u,b,n)
y(1)=b(1);
for i=2:n
    m=0;
    for j=1:i-1
        m=m+l(i,j)*y(j);
    end
    y(i)=b(i)-m;
end
x(n)=y(n)/u(n,n);
for i=n-1:-1:1
    m=0;
    for j=n:-1:i+1
        m=m+u(i,j)*x(j);
    end
    x(i)=(y(i)-m)/u(i,i);
end