MATLABで素数定理を確認する

n_min = 2;
n_max = 1000;

x = n_min:1:n_max;
p_count=zeros(size(x));

pi_n = x ./ log(x);
for i=1:1:length(x)
    p_count(1,i) = length(primes(i));
end

% グラフに描画
figure;
plot(x, pi_n,"k");
hold on;
plot(x, p_count,"k--");
legend('素数定理', '素数の個数');