
format long

A = [ 2.222 16.710 9.612; -1.5611 5.1792 -1.6855; 3.3330 15920 10.333]

b = [.965; 2.714; 7953]

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

P'*L*U

A \ b

AS = sparse(A);

[L,U,P] = lu(AS,0);

full(L)
full(U)
full(P)


