黎曼ζ(2)的导数:ζ'(2)=-1

黎曼 ζ(2)ζ(2)的导数:ζ(2)=1ζ'(2)=-1吗?

http://mathworld.wolfram.com/Glaisher-KinkelinConstant.html
(nx)=lnnnx,(n^{-x})'=-\frac{\ln n}{n^x} ,

ζ(2)=(n=11n2)=ln222+ln332+ln442+...+lnnn2+...-\zeta'(2)=\left(-\sum_{n=1}^{\infty}{\frac{1}{n^2}} \right)'=\frac{\ln2}{2^{2}}+\frac{\ln3}{3^{2}}+\frac{\ln4}{4^{2}}+...+\frac{\ln n}{n^{2}}+...

=π26(12lnAγln2lnπ)=\frac{\pi^{2}}{6}(12\ln A-\gamma-\ln2-\ln\pi)

=0.93754825431584...1,(n)=0.93754825431584...\approx1,(n\rightarrow\infty)

ζ(2)=n=1ln2nn2=1.989280234...2!=2\zeta''(2)=\sum_{n=1}^{\infty}{\frac{\ln^2 n}{n^2}}=1.989280234...\approx2!=2

ζ(2)=n=2ln3nn2=6.000145802...3!=6-\zeta'''(2)=\sum_{n=2}^{\infty}{\frac{\ln^3 n}{n^2}}=6.{\color{red}0}{\color{red}0}{\color{red}0} 145802...\approx3!=6

ζ(4)(2)=n=1ln4nn2=24.001486393...4!=24\zeta^{(4)}(2)=\sum_{n=1}^{\infty}{\frac{\ln^4 n}{n^2}}=24.{\color{red}0}{\color{red}0}1486393...\approx4!=24

ζ(5)(2)=n=1ln5nn2=120.0008243332...5!=120-\zeta^{(5)}(2)=\sum_{n=1}^{\infty}{\frac{\ln^5 n}{n^2}}=120.{\color{red}000}8243332...\approx5!=120

ζ(6)(2)=n=1ln6nn2=720.00012472831...6!=720\zeta^{(6)}(2)=\sum_{n=1}^{\infty}{\frac{\ln^6 n}{n^2}}=720.{\color{red}000}12472831...\approx6!=720

ζ(7)(2)=n=1ln7nn2=5039.999783274...7!=5040-\zeta^{(7)}(2)=\sum_{n=1}^{\infty}{\frac{\ln^7 n}{n^2}}=5039.{\color{red}999}783274...\approx7!=5040

ζ(8)(2)=n=1ln8nn2=40319.99974686...8!=40320\zeta^{(8)}(2)=\sum_{n=1}^{\infty}{\frac{\ln^8 n}{n^2}}=40319.{\color{red}999}74686...\approx8!=40320

ζ(9)(2)=n=1ln9nn2=362879.9998677...9!=362880-\zeta^{(9)}(2)=\sum_{n=1}^{\infty}{\frac{\ln^9 n}{n^2}}=362879.{\color{red}999}8677...\approx9!=362880

太神奇了! \approx 我怀疑是 e,π,Aγe,\pi, A与\gamma 的误差造成的!与收敛的快慢有关。
ζ(k)(2)=(n=11n2)(k)=n=1(lnn)kn2=k!|-\zeta^{(k)}(2)|=\left| \left(-\sum_{n=1}^{\infty}{\frac{1}{n^2}} \right)^{(k)} \right|=\sum_{n=1}^{\infty}{\frac{(\ln n)^{k}}{n^2}}=k!

意味着:
ζ(k)(2)=n=1(lnn)kn2=1(lnx)kx2dx=0xkexdx=k!|-\zeta^{(k)}(2)|=\sum_{n=1}^{\infty}{\frac{(\ln n)^{k}}{n^2}}=\int_{1}^{\infty}\frac{(\ln x)^{k}}{x^2}dx=\int_{0}^{\infty}\frac{x^{k}}{e^{x}}dx=k!

limnn=1(lnn)kn2Δx=1(lnx)kx2dx=k!\lim_{n \rightarrow \infty}{}\sum_{n=1}^{\infty}{\frac{(\ln n)^{k}}{n^2}}\Delta x=\int_{1}^{\infty}\frac{(\ln x)^{k}}{x^2}dx=k!

Δx=ban=1=1\Delta x=\frac{b-a}{n}=\frac{\infty-1}{\infty}=1

ζ(2)=0.93754825431584...1-\zeta'(2)=0.93754825431584...\approx1 误差大是因为收敛得慢计算机算不到无穷远处?后面的误差很小是因为 kk 较大时收敛得快速?还是我想多了? ζ(2)1-\zeta'(2)\ne1 ?
还有:
(x=11ex)=01exdx=1,\left( -\sum_{x=1}^{\infty}{\frac{1}{e^{x}}} \right)'=\int_{0}^{\infty}\frac{1}{e^{x}}dx=1,

ζ(2)=0ln22xdx=1lnxx2dx=11x2dx=1,-\zeta’(2)=\int_{0}^{\infty}\frac{\ln 2}{2^{x}}dx=\int_{1}^{\infty}\frac{\ln x}{x^{2}}dx=\int_{1}^{\infty}\frac{ 1}{x^{2}}dx=1,

这里有一个重要问题:什么情况下级数和等于积分?即什么情况下当 an=f(n)a_{n}=f(n)时,n=1an=1f(x)dx\sum_{n=1}^{\infty}{a_{n}}=\int_{1}^{\infty}f(x)dx?
我在《托马斯微积分》第10版p665页找到了答案:
原文讲积分判别法

{an}\left\{ a_{n} \right\} 是一个正数项序列,假定对 xNx\geq N (正整数), an=f(n),fa_{n}=f(n),fxx 的一个连续,正的,递减函数,则级数 n=NanNf(x)dx\sum_{n=N}^{\infty}{a_{n}} 和 \int_{N}^{\infty}f(x)dx 同时收敛或同时发散。

证明
我们对 N=1N=1 证明,对一般的 NN 证明是类似的。
我们首先假设 ff 是减函数并对所有 nn 使 f(n)=an,f(n)=a_{n} ,(a),(a),矩形 a1+a2++ana_{1}+a_{2}+\cdots +a_{n}面积大于 从x=1x=1x=n+1x=n+1 曲线 y=f(x)y=f(x) 下的面积,于是
1nf(x)dx1n+1f(x)dxa1+a2++an\int_{1}^{n}f(x)dx\leq\int_{1}^{n+1}f(x)dx\leq a_{1}+a_{2}+\cdots +a_{n}

黎曼ζ(2)的导数:ζ'(2)=-1

现在,矩形从往右画改成往左画,图(b)(b)。如果暂时忽略第一个矩形 a1a_{1} ,我们看到

a2+a3++an1nf(x)dxa_{2}+a_{3}+\cdots +a_{n}\leq\int_{1}^{n}f(x)dx
算上 a1,a_{1} ,我们有 a1+a2++ana1+1nf(x)dxa_{1}+a_{2}+\cdots +a_{n}\leq a_{1}+\int_{1}^{n}f(x)dx
组合这些结果得
(1):   1nf(x)dxa1+a2++ana1+1nf(x)dx(1):\,\,\,\int_{1}^{n}f(x)dx\leq a_{1}+a_{2}+\cdots +a_{n}\leq a_{1}+\int_{1}^{n}f(x)dx

1nf(x)dx\int_{1}^{n}f(x)dx 有限,右面的不等式表明 an\sum_{}^{}{a_{n}} 有限,若 1nf(x)dx\int_{1}^{n}f(x)dx 无限,左面的不等式表明 an\sum_{}^{}{a_{n}} 无限。
因此,级数和积分同时有限或无限。

收敛的级数与积分并不一定收敛到同一值。
比如 n=11/n2=π2/6,\sum_{n=1}^{\infty}{1/n^2}=\pi^{2}/6 ,11/x2dx=1\int_{1}^{\infty}1/x^2dx=1

问题来了,什么情况下收敛的级数与积分收敛到同一值?


证明1:
(1)(1) 式知: a1=0a_{1}=0 时,且从 n=2n=2 开始, ana_{n}f(n)f(n) 递减,满足

1nf(x)dx0+a2++an1nf(x)dx\int_{1}^{n}f(x)dx\leq 0+a_{2}+\cdots +a_{n}\leq\int_{1}^{n}f(x)dx


n=1an=1f(x)dx,a1=0 \sum_{n=1}^{\infty}{a_{n}}=\int_{1}^{\infty}f(x)dx,a_{1}=0

敛散性一致,同时收敛或同时发散。

因此, a1=0a_{1}=0 时, n=1an=1f(x)dx\sum_{n=1}^{\infty}{a_{n}}=\int_{1}^{\infty}f(x)dx

(从n=2n=2开始递减)
这就证明了:

ζ(k)(2)=n=1(lnn)kn2=1(lnx)kx2dx=k!|-\zeta^{(k)}(2)|=\sum_{n=1}^{\infty}{\frac{(\ln n)^{k}}{n^2}}=\int_{1}^{\infty}\frac{(\ln x)^{k}}{x^2}dx=k!

证毕
同理得到:ζ(s)=n=1lnnns=1lnxxsdx=1(s1)2,s>1 -\zeta'(s)=\sum_{n=1}^{\infty}{\frac{\ln n}{n^{s}}}=\int_{1}^{\infty}\frac{\ln x}{x^{s}}dx=\frac{1}{(s-1)^{2}},s>1

同理得到:n=0nen=0xexdx=1\sum_{n=0}^{\infty}{\frac{n}{e^{n}}}=\int_{0}^{\infty}\frac{x}{e^{x}}dx=1
因为a0=0a_{0}=0
且有从n=1n=1开始递减
0nf(x)dx0+a1++an0+0nf(x)dx\int_{0}^{n}f(x)dx\leq 0+a_{1}+\cdots +a_{n}\leq0+\int_{0}^{n}f(x)dx

0xexdxn=0nen0xexdx=1\int_{0}^{\infty}\frac{x}{e^{x}}dx\leq\sum_{n=0}^{\infty}{\frac{n}{e^{n}}}\leq\int_{0}^{\infty}\frac{x}{e^{x}}dx=1


证明2:

ζ(2)=n=11n2=0xex1dx=π26\zeta(2)=\sum_{n=1}^{\infty}{\frac{1}{n^2}}=\int_{0}^{\infty}\frac{x}{e^{x}-1}dx=\frac{\pi^{2}}{6}

1ex1=n=11enx\frac{1}{e^x-1}=\sum_{n=1}^{\infty}{\frac{1}{e^{nx}}}

0xex1dx=0n=1xenxdx\int_{0}^{\infty}\frac{x}{e^{x}-1}dx=\int_{0}^{\infty}\sum_{n=1}^{\infty}\frac{x}{e^{nx}}dx

换元,令 x=lnyx=\ln y

0xex1dx=1n=1lnyenlnyd(lny)=1n=1lnyyn+1dy\int_{0}^{\infty}\frac{x}{e^{x}-1}dx=\int_{1}^{\infty}\sum_{n=1}^{\infty}\frac{\ln y}{e^{n\ln y}}d(\ln y)=\int_{1}^{\infty}\sum_{n=1}^{\infty}{\frac{\ln y}{y^{n+1}}}dy

ζ(2)=n=11n2=1n=1lnyyn+1dy=n=11lnyyn+1dy=11lnyyn+1dy\zeta(2)=\sum_{n=1}^{\infty}{\frac{1}{n^2}}=\int_{1}^{\infty}\sum_{n=1}^{\infty}{}\frac{\ln y}{y^{n+1}}dy=\sum_{n=1}^{\infty}{\int_{1}^{\infty}}\frac{\ln y}{y^{n+1}}dy=\int_{1}^{\infty}\int_{1}^{\infty}\frac{\ln y}{y^{n+1}}dy

又因为,1lnyyn+1dy=1n211n2dn=1\int_{1}^{\infty}\frac{\ln y}{y^{n+1}}dy=\frac{1}{n^2} , \int_{1}^{\infty}\frac{1}{n^2}dn=1

所以有
y=1lnyyk+1=1k2,\sum_{y=1}^{\infty}{}\frac{\ln y}{y^{k+1}}=\frac{1}{k^2} ,

k=1k=1时,ζ(2)=n=1lnnn2=1-\zeta'(2)=\sum_{n=1}^{\infty}{}\frac{\ln n}{n^{2}}=1

ζ(s)=n=1lnnns=1lnxxsdx=1(s1)2,s>1-\zeta'(s)=\sum_{n=1}^{\infty}{\frac{\ln n}{n^{s}}}=\int_{1}^{\infty}\frac{\ln x}{x^{s}}dx=\frac{1}{(s-1)^{2}},s>1

待续