function value = sigma ( n )

%% SIGMA computes the sum of the divisors of N.
%
  value = 0.0;

  for i = 1 : n
    if ( mod ( n, i ) == 0)
      value = value + i;
    end
  end

  return
end

