Criba de Eratóstenes Matlab

Escribir un script de Matlab que calcule números primos existentes entre los primeros 1000 número naturales. Para ello se implementará aplicando lo que se conoce como “criba de Eratóstenes”. Este procedimiento consiste en la eliminación sucesiva de los múltiplos de un determinado número, a partir de 2, considerando dicho número inicial como primo si no ha sido eliminado.



function primos=practica1_8()
n=1000;

enteros=2:n;  
primos =[1];  
    while length(enteros)>0  
        sicPrim=enteros(1) ; 
        primos=[primos,sicPrim]; 
        enteros1=[];
        for i=1:length(enteros) 
            if(mod(enteros(i),sicPrim)~=0) 
                enteros1=[enteros1, enteros(i)] 
            end % if  
        end % for 
        enteros=enteros1;  
    end % while   
end % function