euler.m

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

function y=euler(n)
m=1;%Euler primes
p=0;%potential incorrect Euler primes
for k=2:n
    if(isprime(k))
        for x=0:k-2
            z=x^2+x+k;
            if(isprime(z))
                y(m)=z;
                m=m+1;
            else
                p=p+1;
                q(p)=z;
            end
        end
    end
end
p,q