2.2.m

来自「学习使用各种数值方法计算数值代数问题。其中包括计算矩阵的特征值问题」· M 代码 · 共 8 行

for k=1:n
    for i=l+1:n
        p(i,k)=A(i,k)/A(k,k);
        for j=k:m
            A(i,j)=A(i,j)-A(k,j)*p(i,k);
        end
    end
end