// maior_pivo.sci

// dada uma coluna j da matriz A
// devolve o indice i da linha
// com maior pivo

function i = maior_pivo (A, j)
  [nlin, ncol] = size(A);
  
  // observe que comeco a procurar pivos
  // abaixo da diagonal principal
  i = j;
  for k = j:nlin do
    if abs(A(i,j)) < abs(A(k,j))
      i = k;
    end
  end
endfunction
