li.m

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

function [x,y]=li(n)
y(1)=1;
y(2)=1/log(2);
for k=3:n
    y(k)=y(k-1)+1/(log(k));
end
x=y(n);

k=1:n;
plot(k,y(k));