conjuncture_compare.m

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

function a=conjuncture_compare(n)

m=0;
for k=1:n
    if(isprime(k))
        m=m+1;
    end
    y1(k)=m;
end
k=1:n;
plot(k,y1(k),'r');hold on; %red

y2(1)=1;
for k=2:n
    y2(k)=k/log(k);
end
k=1:n;
plot(k,y2(k),'g');hold on; %green

y3(1)=1;
y3(2)=1/log(2);
for k=3:n
    y3(k)=y3(k-1)+1/(log(k));
end
k=1:n;
plot(k,y3(k),'b');    %blue

legend('actual prime numbers','Gauss conjuncture','li')