%例2  求[1，100]区间的所有素数。
x=1:100;
k=isprime(x);
k1=find(k);
p=x(k1)