lu.m

CG共轭梯度法迭代解矩阵方程

function[L,U]=lu(A%不选主元函数
[m,n]=size(A);
U=A;
L=eye(m,n);
for k=1:m-1
  for j=k+1:m
      L[j,k]=U[j,k]/U[k,k];
      U[j,k:m]=U[j,k:m]-L[j,k]*U[k,k:m];
  end
end