k=5000;i=2;h=3:k;while(h[1]<sqrt(k)){i=c(i,h[1]);h=h[apply(outer(h,i,'%%'),1,prod)!=0]};c(i,h)