data LJX;

do i=1 to 1000;

count=0;

m=sqrt(i);

do j=1 to m;

if mod(i,j)=0 then count+1;

end;

if count=1;

output;

drop count j;

run;

proc print;

run;