MATLABで素数定理を確認する

primesを使えば簡単。

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
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('素数定理', '素数の個数');

コメント

タイトルとURLをコピーしました