eratos(m)

function x = eratos(m)
%eratos  nombres premiers <= m
%        eratos est 5 fois plus
%        rapide que primes0
%        mais utilisez plutt
%        la  version Matlab (primes)

x = 1:m; x(1)=0;
for i = 2:fix(m/2)
  if x(i) ~= 0
    x(2*i:i:m) = zeros(size(2*i:i:m));
  end;
end;
x = x(x>0); % on retire tous les
            % zros du tableau ;