# solve for Z - >[L][Z]=[C]

Z(1) = b(1);
for i = 2 : n
  sum = 0;
  for j = i - 1:-1:1
    sum = sum + (L(i,j) * Z(j));
  endfor
  Z(i) = b(i) - sum;
endfor

##disp(Z);

# solve for X -> [U][X]=[Z]

X(n) = Z(n) / U(n,n);
for i = n - 1: -1: 1
  sum = 0;
  for j = i + 1:1:n
    sum = sum + (U(i,j)*X(j));
  endfor
  X(i) = (Z(i) - sum) / U(i,i);
endfor
disp(X);