Factor A = LU using the Doolittle algorithm. L is lower triangular with 1s on the diagonal; U is upper triangular. Then solve Ly = b and Ux = y.