// seidel_decomp.sci

function [N, P] = seidel_decomp(A)
  [nl, nc] = size(A);
  n = nl;
  
  // construindo P do Metodo de Gauss-Seidel
  P = zeros (n, n);
  
  for i = 1:n do
    for j = (i + 1):n do
      P(i,j) = -A(i,j);
      A(i,j) = 0;
    end
  end
  
  // o que sobra fica na N
  N = A;
endfunction

