-->A = [1, 2, 3; 4, 5, 6]; C = A'*A;
-->[L, U, E] = lu(C)
E =
0. 0. 1.
1. 0. 0.
0. 1. 0.
U =
27. 36. 45.
0. - 0.6666667 - 1.3333333
0. 0. 0.
L =
1. 0. 0.
0.6296296 1. 0.
0.8148148 0.5 1.
-->L*U
ans =
27. 36. 45.
17. 22. 27.
22. 29. 36.
-->E*C
ans =
27. 36. 45.
17. 22. 27.
22. 29. 36.
-->E'*L*U
ans =
17. 22. 27.
22. 29. 36.
27. 36. 45.
-->C
C =
17. 22. 27.
22. 29. 36.
27. 36. 45.
-->E'*L*U - C
ans =
1.0D-14 *
- 0.3552714 0. 0.
0. 0. 0.
0. 0. 0.
-->E'*L*U - C == zeros(2, 3)
ans =
F
-->E'*L*U - C == zeros(3, 3)
ans =
F T T
T T T
T T T
-->E'*L*U == C
ans =
F T T
T T T
T T T