a =. ; b =. ; c =.

AL = {{1, 0, 0}, {a, 1, 0}, {0, b, 1}}

AU = {{3, 1, 0}, {0, c, 1}, {0, 0, d}}

B = AL . AU

a = 1; b = -1; c = 3;  d = 4

B

LinearSolve[AL, {4, 8, 0}]

LinearSolve[AU, {4, 4, 4}]

M = {{3, 0, 0}, {3, 4, 0}, {0, -3, 3}}

NN = {{0, 1, 0}, {0, 0, 1}, {0, 0, 0}}

(* Logo,  
      a matriz de iteração C é dada por *)

CC = -Inverse[M] . NN

g = Inverse[M] . {0, 0, 1}