ge.m

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

function [l,u]=GE(A)
r=rank(A);
u=A;
l=eye(r);
   for p=1:(r-1)
       for m=(p+1):r
           l(m,p)=u(m,p)/u(p,p);
           for n=1:r
               y=u(m,n)-l(m,p)*u(p,n);
               u(m,n)=y;
           end
       end
  end