gedebahe_num.m

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

function y=gedebahe_num(n)
for p=2:2:n
    m=0;
    for k=1:p/2
        if(isprime(k)&&isprime(p-k)) m=m+1;
        end
    end
    y(p/2)=m;
end

p=2:2:n
stem(p,y)