求解20000以内的所有质数，并做出在等距区间内质数的个数的直方图（以1000为间距）

b=1;
s=zeros(1,20000);%建立一个储存空间
for a=1:20000;
if a<=20000
     x=isprime(a);%isprime判断素数函数,如果是a是素数，则输出为1
     if x==1%如果x等于1，则说明a是素数，则将a存入s
         s(b)=a;
         b=b+1;
     else 
     end
end
end
s(s==0)=[];
hist(s,20);%将ｓ中所有数据画柱形图，分成20个柱形图