primes.m

来自「关于质数的一些程序」· M 代码 · 共 10 行

function [y,z]=primes(n)
m=1;
for k=1:n
    if(isprime(k))
        y(m)=k;m=m+1;
    end
    z(k)=m-1;
end
k=1:n;
steps(k,z)