matlab编写M文件用while循环找出从100到200之间的素数

a=[];
i=100;
while i<=200,
j=1;
while j<i,
if gcd(j,i)~=1,
a=[a,i]; %i为合数
end
j=j+1;
end
i=i+1;
end
b=100:200; %b为全集
disp(setdiff(b,a)); %打印质数集

primes命令

i=100:200;
i(isprime(i))