Question

有人会解matlab这道题么，后天就得交作业，求大神解答ORZ





Answer 1

% 提示用户输入
N=input('Input a number: ');

% 计算估计密度
N=(2:N);
pi_n=N./log(N);

% 求2-N的所有整数n，求素数密度
% 采用的方法是，看能否被2-sqrt(n)之间的某个整数整除，如果可以整除
% 则不是素数
pi_n_exact = ones( size(N) );
for n=N
if n>2, pi_n_exact(n-1) = pi_n_exact(n-2)+1; end
for i=2:sqrt(n)
if rem(n,i) == 0
pi_n_exact(n-1) = pi_n_exact(n-1) - 1;
break
end
end
end

% 绘图
clf
area(N,pi_n_exact,'facec','c')
hold on
plot(N,pi_n)
xlabel('n')
ylabel('Density of primes')
legend('Estimated \pi(n)','\pi(n)',0)