A = [
  1   3   5
  2   4   7
  1   1   0];

  [L,U,P] = lu(A)