result.m

来自「用Crout分解法分解系数矩阵」· M 代码 · 共 15 行

function result(A,B,n)
[L U]=crout(A);
%令Y=UX，先求解Y
Y(1)=B(1)/L(1,1);
for i=2:n
    M=Y(i-1)*L(i,i-1);
    Y(i)=(B(i)-M)/L(i,i);
end
%得到Y，求解X
X(n)=Y(n);
for i=n-1:-1:1
    X(i)=Y(i)-U(i,i+1)*X(i+1);
end
Y
X