primes.m

关于质数的一些程序

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)