LU分解版本1

对系数矩阵A进行LU分解，没有采用原地存储方式，没有进行选主元

function ex2_7_1
%对系数矩阵A进行LU分解【自己写的版本，没有采用原地存储方式，没有进行选主元】
%将L矩阵和U矩阵输出，

 A = [10 -7 0;
      -3 2 6;
      5 -1 5];
 %对系数矩阵A进行LU分解，最终得到L 和 U矩阵
 %初始化L和U矩阵.
 L = [0 0 0;
      0 0 0;
      0 0 0]
 U = [0 0 0;
      0 0 0;
      0 0 0]
  [n,n] = size(A);
  for k = 1:1:n
      i = k+1:n;
      %计算乘子，将乘子填入到L矩阵的第（2、3）行，第1列
      L(i,k) = A(i,k)/A(k,k);
      %更新系数矩阵A，就是将A的第一行乘以相应的乘子的相反数，加到第二行和第三行
      j = k:n;
      A(i,j) = A(i,j) - L(i,k)*A(k,j);
  end
  U = A;
  L 
  U