eratosten

function [result_vector] = eratosten(max_number)
stop_variable =  max_number;

if (max_number==1)
    disp('1 is not allowed');
   return;
end

if(max_number<4)

for l=2:max_number
    result_vector(l-1)=l;

end

return;

end

if(max_number==4)
result_vector(1)=2;
result_vector(2)=3;
return;
end


for i=2:max_number
        test_vector(i-1)=i;
end

p=test_vector(1);
j=1;
    while(p^2 < stop_variable)
       k=1;
       result_vector(j)=p;
       j=j+1;


        for i=1:length(test_vector)

            if( mod( test_vector(i), p )~=0)
            test_internal(k)=test_vector(i);
            k=k+1;
            end

        end


        for i=1:k-1
            test_vector(i)=test_internal(i);

        end
          %%izbegavanje upisa iz originalnog test vektora

        for i=k:length(test_vector)
            test_vector(i)=0;

        end


         p = test_vector(1);


    end


     for i=1:k-1

         result_vector(j)=test_vector(i);
         j=j+1;
       end



end