e0325.m

来自「matlab的一些学习资料和matlab的部分经典源码」· M 代码 · 共 23 行

%1
m=input('m=');
p=1:m;p(1)=0;
for i=2:sqrt(m)
    for j=2*i:i:m
        p(j)=0;
    end
end
n=find(p~=0);
p(n)
%2
m=input('m=');
p=2:m;
for i=2:sqrt(m)
    n=find(rem(p,i)==0&p~=i);
    p(n)=[];
end
p