O enunciado nos dá um sistema linear 4\times 4 e pede para resolvê-lo via decomposição LU (método de Doolittle). Vamos fatorar a matriz dos coeficientes em A = LU, onde L é unitária (1 na diagonal) e U é triangular superior, e depois resolver Ly=b e Ux=y.