getf ("gausspv.sci");

// declarando as matrizes do problema
A = [6, -2, 2, 4; 12, -8, 6, 10; 3, -13, 9, 3; -6, 4, 1, -18]
b = [12; 34; 27; -38]

// resolvendo por MEG com pivotamento
[p, B, x] = gausspv (A, b)

// soh por curiosidade vamos obter 
// as matrizes P, L e U.

getf ("troca_linhas.sci");

P = eye (4, 4);
for i = 1:4 do
  P = troca_linhas(P, i, p(i));
end

L = eye (4, 4);
for i = 2:4 do
  L(i, 1:(i - 1)) = B(i, 1:(i - 1));
end

U = zeros (4,4);
for i = 1:4 do
  U(i, i:4) = B(i, i:4);
end

C1 = P * A
C2 = L * U
  
