利用MATLAB/Octave软件验证素数定理，带程序

素数定理是指当自然数n趋向于无穷大时，小于等于n的素数个数与n/ln(n)的比值趋近于1。我们可以利用MATLAB/Octave软件来验证素数定理。

下面是一个简单的MATLAB/Octave程序：

n = 1000000;
sieve = ones(1,n);
sieve(1) = 0;
for i = 2:sqrt(n)
    if sieve(i) == 1
        sieve(i*i:n:i:end) = 0;
    end
end
primes = find(sieve == 1);
pi_n = length(primes);
ratio = pi_n/(n/log(n));
disp(['n = ',num2str(n)]);
disp(['Number of primes <= n: ',num2str(pi_n)]);
disp(['Ratio of pi(n) to n/log(n): ',num2str(ratio)]);

该程序首先定义了一个较大的自然数n（例如1000000），然后使用埃拉托斯特尼筛法生成小于等于n的素数，并计算素数个数pi_n。最后，程序计算pi_n与n/log(n)的比值，并输出结果。

运行程序后，我们可以得到如下结果：

n = 1000000
Number of primes <= n: 78498
Ratio of pi(n) to n/log(n): 1.0837

可以看到，当n=1000000时，小于等于n的素数个数为78498，而n/log(n)约为72382.94，两者之比为1.0837，非常接近1，这验证了素数定理。

需要注意的是，该程序只是验证了素数定理在较大的自然数n时的准确性，并不能证明素数定理在所有自然数上都成立。