Das Polynom
P
n
(
x
)
:=
1
n
!
⋅
d
n
d
x
n
[
x
n
(
1
−
x
)
n
]
=:
∑
ν
=
0
n
a
ν
x
ν
{\displaystyle P_{n}(x):={\frac {1}{n!}}\cdot {\frac {{\rm {d}}^{n}}{{\rm {d}}x^{n}}}[x^{n}(1-x)^{n}]=:\sum _{\nu =0}^{n}a_{\nu }x^{\nu }}
hat ausschließlich ganzzahlige Koeffizienten
a
ν
∈
Z
{\displaystyle a_{\nu }\in \mathbb {Z} }
.
Nach der Leibniz'schen Formel gilt
P
n
(
x
)
=
1
n
!
∑
ν
=
0
n
(
n
ν
)
n
(
n
−
1
)
⋯
(
n
−
ν
+
1
)
x
n
−
ν
n
⋯
(
n
−
(
n
−
ν
)
+
1
)
(
1
−
x
)
n
−
(
n
−
ν
)
(
−
1
)
n
−
ν
=
1
n
!
∑
ν
=
0
n
(
n
ν
)
n
!
(
n
−
ν
)
!
x
n
−
ν
n
!
ν
!
(
1
−
x
)
ν
(
−
1
)
n
−
ν
=
∑
ν
=
0
n
(
n
ν
)
2
(
−
1
)
n
−
ν
x
n
−
ν
(
1
−
x
)
ν
,
{\displaystyle {\begin{array}{lll}P_{n}(x)&=&{\frac {1}{n!}}\sum \limits _{\nu =0}^{n}{n \choose \nu }n(n-1)\cdots (n-\nu +1)x^{n-\nu }n\cdots (n-(n-\nu )+1)(1-x)^{n-(n-\nu )}(-1)^{n-\nu }\\&=&{\frac {1}{n!}}\sum \limits _{\nu =0}^{n}{n \choose \nu }{\frac {n!}{(n-\nu )!}}x^{n-\nu }{\frac {n!}{\nu !}}(1-x)^{\nu }(-1)^{n-\nu }=\sum \limits _{\nu =0}^{n}{n \choose \nu }^{2}(-1)^{n-\nu }x^{n-\nu }(1-x)^{\nu }\ ,\\\end{array}}}
d.h.,
P
n
{\displaystyle P_{n}}
hat nur ganzzahlige Koeffizienten.
◻
{\displaystyle \Box }
Ist
f
:
Q
:=
[
0
,
1
]
2
→
R
{\displaystyle f:Q:=[0,1]^{2}\rightarrow \mathbb {R} }
stetig, so existiert das Integral
∫
Q
f
(
x
,
y
)
1
−
x
y
d
x
d
y
<
∞
{\displaystyle \int _{Q}{\frac {f(x,y)}{1-xy}}{\rm {d}}x{\rm {d}}y<\infty }
.
Da
f
{\displaystyle f}
auf dem Kompaktum
Q
{\displaystyle Q}
stetig, also beschränkt ist, muss
∫
Q
1
1
−
x
y
d
x
d
y
<
∞
{\displaystyle \int _{Q}{\frac {1}{1-xy}}{\rm {d}}x{\rm {d}}y<\infty }
gezeigt werden. Dies ist offenbar äquivalent zu
∫
Q
1
1
−
(
1
−
u
)
(
1
−
v
)
d
u
d
v
<
∞
{\displaystyle \int _{Q}{\frac {1}{1-(1-u)(1-v)}}{\rm {d}}u{\rm {d}}v<\infty }
mit
x
=
1
−
u
{\displaystyle x=1-u}
,
y
=
1
−
v
{\displaystyle y=1-v}
. Weil der Integrand nur in
u
=
v
=
0
{\displaystyle u=v=0}
singulär wird, brauchen wir die Endlichkeit nur auf dem Viertelkreis
V
:=
{
(
r
cos
φ
,
r
sin
φ
)
|
0
<
r
≤
1
,
0
≤
φ
≤
π
2
}
{\displaystyle V:=\{(r\cos \varphi ,r\sin \varphi )\ |\ 0<r\leq 1,\ 0\leq \varphi \leq {\frac {\pi }{2}}\}}
zu zeigen. Hier können wir Polarkoordinaten einführen und erhalten
∫
V
1
1
−
(
1
−
u
)
(
1
−
v
)
d
u
d
v
=
∫
0
π
/
2
∫
0
1
r
1
−
(
1
−
r
cos
φ
)
(
1
−
r
sin
φ
)
d
r
d
φ
=
∫
0
π
/
2
∫
0
1
1
sin
φ
+
cos
φ
−
r
sin
φ
cos
φ
d
r
d
φ
=
∫
0
π
/
4
∫
0
1
1
sin
φ
(
1
−
r
cos
φ
)
+
cos
φ
d
r
d
φ
+
∫
π
/
4
π
/
2
∫
0
1
1
cos
φ
(
1
−
r
sin
φ
)
+
sin
φ
d
r
d
φ
≤
∫
0
π
/
4
∫
0
1
1
cos
φ
d
r
d
φ
+
∫
π
/
4
π
/
2
∫
0
1
1
sin
φ
d
r
d
φ
≤
2
π
4
+
2
π
4
<
∞
,
{\displaystyle {\begin{array}{lll}\int _{V}{\frac {1}{1-(1-u)(1-v)}}{\rm {d}}u{\rm {d}}v&=&\int _{0}^{\pi /2}\int _{0}^{1}{\frac {r}{1-(1-r\cos \varphi )(1-r\sin \varphi )}}{\rm {d}}r{\rm {d}}\varphi =\int _{0}^{\pi /2}\int _{0}^{1}{\frac {1}{\sin \varphi +\cos \varphi -r\sin \varphi \cos \varphi }}{\rm {d}}r{\rm {d}}\varphi \\&=&\int _{0}^{\pi /4}\int _{0}^{1}{\frac {1}{\sin \varphi (1-r\cos \varphi )+\cos \varphi }}{\rm {d}}r{\rm {d}}\varphi +\int _{\pi /4}^{\pi /2}\int _{0}^{1}{\frac {1}{\cos \varphi (1-r\sin \varphi )+\sin \varphi }}{\rm {d}}r{\rm {d}}\varphi \\&\leq &\int _{0}^{\pi /4}\int _{0}^{1}{\frac {1}{\cos \varphi }}{\rm {d}}r{\rm {d}}\varphi +\int _{\pi /4}^{\pi /2}\int _{0}^{1}{\frac {1}{\sin \varphi }}{\rm {d}}r{\rm {d}}\varphi \leq {\sqrt {2}}{\frac {\pi }{4}}+{\sqrt {2}}{\frac {\pi }{4}}<\infty \ ,\end{array}}}
was zu zeigen war.
◻
{\displaystyle \Box }
Wir untersuchen im Folgenden das Integral
I
n
:=
∫
0
1
∫
0
1
−
ln
x
y
1
−
x
y
P
n
(
x
)
P
n
(
y
)
d
x
d
y
=
∑
r
,
s
=
0
n
a
r
a
s
∫
0
1
∫
0
1
x
r
y
s
1
−
x
y
(
−
ln
x
y
)
d
x
d
y
.
{\displaystyle I_{n}:=\int _{0}^{1}\int _{0}^{1}{\frac {-\ln xy}{1-xy}}P_{n}(x)P_{n}(y){\rm {d}}x{\rm {d}}y=\sum _{r,s=0}^{n}a_{r}a_{s}\int _{0}^{1}\int _{0}^{1}{\frac {x^{r}y^{s}}{1-xy}}(-\ln xy){\rm {d}}x{\rm {d}}y\ .}
Die Wohldefiniertheit wird im folgenden Lemma formuliert:
Für alle
r
,
s
∈
N
0
{\displaystyle r,s\in \mathbb {N} _{0}}
gilt
∫
0
1
∫
0
1
x
r
y
s
1
−
x
y
(
−
ln
x
y
)
d
x
d
y
=
{
1
r
−
s
∑
ν
=
s
+
1
r
1
ν
2
,
falls
r
>
s
,
2
[
ζ
(
3
)
−
∑
ν
=
1
r
1
ν
3
]
,
falls
r
=
s
,
1
s
−
r
∑
ν
=
r
+
1
s
1
ν
2
,
falls
r
<
s
.
{\displaystyle \int _{0}^{1}\int _{0}^{1}{\frac {x^{r}y^{s}}{1-xy}}(-\ln xy){\rm {d}}x{\rm {d}}y=\left\{{\begin{array}{lll}{\frac {1}{r-s}}\sum _{\nu =s+1}^{r}{\frac {1}{\nu ^{2}}}\ ,&{\text{falls}}&r>s\ ,\\2[\zeta (3)-\sum _{\nu =1}^{r}{\frac {1}{\nu ^{3}}}]\ ,&{\text{falls}}&r=s\ ,\\{\frac {1}{s-r}}\sum _{\nu =r+1}^{s}{\frac {1}{\nu ^{2}}}\ ,&{\text{falls}}&r<s\ .\\\end{array}}\right.}
Insbesondere ist
I
n
{\displaystyle I_{n}}
wohldefiniert.
Die Wohldefiniertheit ist lediglich für
r
=
s
=
0
{\displaystyle r=s=0}
zu zeigen. Sei
m
∈
N
{\displaystyle m\in \mathbb {N} }
und
V
:=
{
(
R
cos
φ
,
R
sin
φ
)
|
R
∈
(
0
,
2
]
,
φ
∈
(
0
,
π
2
)
}
{\displaystyle V:=\{(R\cos \varphi ,R\sin \varphi )\ |\ R\in (0,{\sqrt {2}}],\ \varphi \in (0,{\frac {\pi }{2}})\}}
. Wegen
|
ln
x
|
m
≤
c
x
{\displaystyle |\ln x|^{m}\leq {\frac {c}{\sqrt {x}}}}
für
x
∈
(
0
,
2
]
{\displaystyle x\in (0,2]}
ist
∫
V
|
ln
x
y
|
m
d
x
d
y
≤
c
∫
V
1
x
y
d
x
d
y
=
c
∫
0
π
/
2
∫
0
2
R
R
cos
φ
R
sin
φ
d
R
d
φ
=
2
c
∫
0
π
/
2
1
cos
φ
sin
φ
d
φ
=
2
c
∫
0
π
/
4
φ
sin
φ
⋅
1
φ
cos
φ
d
φ
+
2
c
∫
π
/
4
π
/
2
π
2
−
φ
cos
φ
⋅
1
π
2
−
φ
sin
φ
d
φ
≤
C
[
∫
0
π
/
4
1
φ
d
φ
+
∫
π
/
4
π
/
2
1
π
2
−
φ
d
φ
]
<
∞
.
{\displaystyle {\begin{array}{lll}\int _{V}|\ln xy|^{m}{\rm {d}}x{\rm {d}}y&\leq &c\int _{V}{\frac {1}{\sqrt {xy}}}{\rm {d}}x{\rm {d}}y=c\int _{0}^{\pi /2}\int _{0}^{\sqrt {2}}{\frac {R}{{\sqrt {R\cos \varphi }}{\sqrt {R\sin \varphi }}}}{\rm {d}}R{\rm {d}}\varphi ={\sqrt {2}}c\int _{0}^{\pi /2}{\frac {1}{\sqrt {\cos \varphi \sin \varphi }}}{\rm {d}}\varphi \\&=&{\sqrt {2}}c\int _{0}^{\pi /4}{\sqrt {\frac {\varphi }{\sin \varphi }}}\cdot {\frac {1}{{\sqrt {\varphi }}{\sqrt {\cos \varphi }}}}{\rm {d}}\varphi +{\sqrt {2}}c\int _{\pi /4}^{\pi /2}{\sqrt {\frac {{\frac {\pi }{2}}-\varphi }{\cos \varphi }}}\cdot {\frac {1}{{\sqrt {{\frac {\pi }{2}}-\varphi }}{\sqrt {\sin \varphi }}}}{\rm {d}}\varphi \\&\leq &C[\int _{0}^{\pi /4}{\frac {1}{\sqrt {\varphi }}}{\rm {d}}\varphi +\int _{\pi /4}^{\pi /2}{\frac {1}{\sqrt {{\frac {\pi }{2}}-\varphi }}}{\rm {d}}\varphi ]<\infty \ .\\\end{array}}}
Nach Lemma Nummer 2 ist jedes Integral
∫
0
1
∫
0
1
x
r
y
s
1
−
x
y
|
ln
x
y
|
m
d
x
d
y
{\displaystyle \int _{0}^{1}\int _{0}^{1}{\frac {x^{r}y^{s}}{1-xy}}|\ln xy|^{m}{\rm {d}}x{\rm {d}}y}
und somit
I
n
{\displaystyle I_{n}}
wohldefiniert. Um den Wert dieses Integrals auszurechnen, müssen wir die Funktion
I
(
t
)
:=
∫
0
1
∫
0
1
x
r
+
t
y
s
+
t
1
−
x
y
d
x
d
y
,
t
≥
0
,
{\displaystyle I(t):=\int _{0}^{1}\int _{0}^{1}{\frac {x^{r+t}y^{s+t}}{1-xy}}{\rm {d}}x{\rm {d}}y\ ,\ t\geq 0\ ,}
welche nach Lemma Nummer 2 wohldefiniert ist, genauer untersuchen. Und zwar ist
I
(
t
)
=
∫
0
1
∫
0
1
∑
k
=
0
∞
x
r
+
t
+
k
y
s
+
t
+
k
d
x
d
y
nach geometrischer Reihe
=
∑
k
=
0
∞
∫
0
1
∫
0
1
x
r
+
t
+
k
y
s
+
t
+
k
d
x
d
y
nach Satz von Beppo Levi
=
∑
k
=
0
∞
(
∫
0
1
x
r
+
t
+
k
d
x
)
(
∫
0
1
y
s
+
t
+
k
d
y
)
=
∑
k
=
0
∞
1
r
+
t
+
k
+
1
⋅
1
s
+
t
+
k
+
1
.
{\displaystyle {\begin{array}{lll}I(t)&=&\int _{0}^{1}\int _{0}^{1}\sum _{k=0}^{\infty }x^{r+t+k}y^{s+t+k}{\rm {d}}x{\rm {d}}y{\text{ nach geometrischer Reihe}}\\&=&\sum _{k=0}^{\infty }\int _{0}^{1}\int _{0}^{1}x^{r+t+k}y^{s+t+k}{\rm {d}}x{\rm {d}}y{\text{ nach Satz von Beppo Levi}}\\&=&\sum _{k=0}^{\infty }(\int _{0}^{1}x^{r+t+k}{\rm {d}}x)(\int _{0}^{1}y^{s+t+k}{\rm {d}}y)\\&=&\sum _{k=0}^{\infty }{\frac {1}{r+t+k+1}}\cdot {\frac {1}{s+t+k+1}}\ .\\\end{array}}}
Der Trick besteht nun darin, dass man
I
{\displaystyle I}
im Nullpunkt (rechtsseitig) differenzieren darf und daraus die entsprechenden Formeln erhält. Es ist nämlich für
h
>
0
{\displaystyle h>0}
|
I
(
h
)
−
I
(
0
)
h
−
∫
0
1
∫
0
1
x
r
y
s
ln
x
y
1
−
x
y
d
x
d
y
|
≤
∫
0
1
∫
0
1
|
x
r
+
h
y
s
+
h
−
x
r
y
s
(
1
−
x
y
)
h
−
x
r
y
s
ln
x
y
1
−
x
y
|
d
x
d
y
=
∫
0
1
∫
0
1
x
r
y
s
1
−
x
y
|
(
x
y
)
h
−
1
h
−
ln
x
y
|
d
x
d
y
=
∫
0
1
∫
0
1
x
r
y
s
1
−
x
y
|
h
(
ln
x
y
)
2
2
(
x
y
)
ξ
|
d
x
d
y
mit
ξ
=
ξ
(
h
,
x
,
y
)
∈
(
0
,
h
)
≤
h
2
∫
0
1
∫
0
1
(
ln
x
y
)
2
1
−
x
y
x
r
y
s
d
x
d
y
→
0
f
u
¨
r
h
↘
0
{\displaystyle {\begin{array}{lll}\left|{\frac {I(h)-I(0)}{h}}-\int _{0}^{1}\int _{0}^{1}{\frac {x^{r}y^{s}\ln xy}{1-xy}}{\rm {d}}x{\rm {d}}y\right|&\leq &\int _{0}^{1}\int _{0}^{1}\left|{\frac {x^{r+h}y^{s+h}-x^{r}y^{s}}{(1-xy)h}}-{\frac {x^{r}y^{s}\ln xy}{1-xy}}\right|{\rm {d}}x{\rm {d}}y\\&=&\int _{0}^{1}\int _{0}^{1}{\frac {x^{r}y^{s}}{1-xy}}\left|{\frac {(xy)^{h}-1}{h}}-\ln xy\right|{\rm {d}}x{\rm {d}}y\\&=&\int _{0}^{1}\int _{0}^{1}{\frac {x^{r}y^{s}}{1-xy}}\left|{\frac {h(\ln xy)^{2}}{2}}(xy)^{\xi }\right|{\rm {d}}x{\rm {d}}y\ {\text{mit}}\ \xi =\xi (h,x,y)\in (0,h)\\&\leq &{\frac {h}{2}}\int _{0}^{1}\int _{0}^{1}{\frac {(\ln xy)^{2}}{1-xy}}x^{r}y^{s}{\rm {d}}x{\rm {d}}y\rightarrow 0\ \mathrm {f{\ddot {u}}r} \ h\searrow 0\\\end{array}}}
wegen
∫
0
1
∫
0
1
(
ln
x
y
)
2
1
−
x
y
x
r
y
s
d
x
d
y
<
∞
{\displaystyle \int _{0}^{1}\int _{0}^{1}{\frac {(\ln xy)^{2}}{1-xy}}x^{r}y^{s}{\rm {d}}x{\rm {d}}y<\infty }
, wie wir oben gesehen haben.
Insgesamt haben wir somit
I
′
(
0
)
=
∫
0
1
∫
0
1
x
r
y
s
ln
x
y
1
−
x
y
d
x
d
y
{\displaystyle I'(0)=\int _{0}^{1}\int _{0}^{1}{\frac {x^{r}y^{s}\ln xy}{1-xy}}{\rm {d}}x{\rm {d}}y}
. Im Falle
r
=
s
{\displaystyle r=s}
gilt
I
′
(
0
)
=
d
d
t
[
∑
k
=
0
∞
1
(
r
+
t
+
k
+
1
)
2
]
t
=
0
=
∑
k
=
0
∞
[
d
d
t
1
(
r
+
t
+
k
+
1
)
2
]
t
=
0
=
−
2
∑
k
=
0
∞
1
(
r
+
k
+
1
)
3
=
−
2
[
ζ
(
3
)
−
∑
ν
=
1
r
1
ν
3
]
,
{\displaystyle {\begin{array}{lll}I'(0)&=&{\frac {\rm {d}}{{\rm {d}}t}}\left[\sum _{k=0}^{\infty }{\frac {1}{(r+t+k+1)^{2}}}\right]_{t=0}=\sum _{k=0}^{\infty }\left[{\frac {\rm {d}}{{\rm {d}}t}}{\frac {1}{(r+t+k+1)^{2}}}\right]_{t=0}\\&=&-2\sum _{k=0}^{\infty }{\frac {1}{(r+k+1)^{3}}}=-2[\zeta (3)-\sum _{\nu =1}^{r}{\frac {1}{\nu ^{3}}}]\ ,\\\end{array}}}
weil wir wegen der normalen Konvergenz der Zeta-Funktion auf
Re
s
>
1
{\displaystyle \operatorname {Re} s>1}
gliedweise differenzieren dürfen. Im Falle
r
>
s
{\displaystyle r>s}
(und analog
r
<
s
{\displaystyle r<s}
) ist
I
(
t
)
=
1
r
−
s
∑
k
=
0
∞
[
1
s
+
1
+
t
+
k
−
1
r
+
1
+
t
+
k
]
=
1
r
−
s
∑
ν
=
s
+
1
r
1
ν
+
t
{\displaystyle I(t)={\frac {1}{r-s}}\sum _{k=0}^{\infty }\left[{\frac {1}{s+1+t+k}}-{\frac {1}{r+1+t+k}}\right]={\frac {1}{r-s}}\sum _{\nu =s+1}^{r}{\frac {1}{\nu +t}}}
und somit
I
′
(
0
)
=
d
d
t
[
1
r
−
s
∑
ν
=
s
+
1
r
1
ν
+
t
]
t
=
0
=
−
1
r
−
s
∑
ν
=
s
+
1
r
1
ν
2
,
{\displaystyle {\begin{array}{lll}I'(0)&=&{\frac {\rm {d}}{{\rm {d}}t}}\left[{\frac {1}{r-s}}\sum _{\nu =s+1}^{r}{\frac {1}{\nu +t}}\right]_{t=0}=-{\frac {1}{r-s}}\sum _{\nu =s+1}^{r}{\frac {1}{\nu ^{2}}}\ ,\\\end{array}}}
wie behauptet wurde.
◻
{\displaystyle \Box }
Es gilt
d
n
:=
kgV
(
1
,
…
,
n
)
=
∏
p
≤
n
p
⌊
ln
n
ln
p
⌋
≤
n
π
(
n
)
{\displaystyle d_{n}:=\operatorname {kgV} (1,\ldots ,n)=\prod _{p\leq n}p^{\lfloor {\frac {\ln n}{\ln p}}\rfloor }\leq n^{\pi (n)}}
.
Jedes
r
∈
{
1
,
…
,
n
}
{\displaystyle r\in \{1,\ldots ,n\}}
hat nach dem Fundamentalsatz der Arithmetik eine Primfaktorzerlegung
r
=
∏
p
≤
n
p
α
(
p
,
r
)
≤
n
{\displaystyle r=\prod _{p\leq n}p^{\alpha (p,r)}\leq n}
mit
α
(
p
,
r
)
∈
N
0
{\displaystyle \alpha (p,r)\in \mathbb {N} _{0}}
. Dann muss
p
α
(
p
,
r
)
≤
n
{\displaystyle p^{\alpha (p,r)}\leq n}
, also
α
(
p
,
r
)
≤
ln
n
ln
p
{\displaystyle \alpha (p,r)\leq {\frac {\ln n}{\ln p}}}
und somit wegen der Ganzzahligkeit
α
(
p
,
r
)
≤
⌊
ln
n
ln
p
⌋
{\displaystyle \alpha (p,r)\leq \lfloor {\frac {\ln n}{\ln p}}\rfloor }
gelten.
Daraus folgt
d
n
≤
∏
p
∈
P
p
⌊
ln
n
ln
p
⌋
=
∏
p
≤
n
p
⌊
ln
n
ln
p
⌋
{\displaystyle d_{n}\leq \prod _{p\in \mathbb {P} }p^{\lfloor {\frac {\ln n}{\ln p}}\rfloor }=\prod _{p\leq n}p^{\lfloor {\frac {\ln n}{\ln p}}\rfloor }}
. Die umgekehrte Ungleichung folgt, weil
p
⌊
ln
n
ln
p
⌋
≤
n
{\displaystyle p^{\lfloor {\frac {\ln n}{\ln p}}\rfloor }\leq n}
gilt. Also ist
d
n
=
∏
p
≤
n
p
⌊
ln
n
ln
p
⌋
≤
∏
p
≤
n
p
ln
n
ln
p
=
n
π
(
n
)
{\displaystyle d_{n}=\prod _{p\leq n}p^{\lfloor {\frac {\ln n}{\ln p}}\rfloor }\leq \prod _{p\leq n}p^{\frac {\ln n}{\ln p}}=n^{\pi (n)}}
.
◻
{\displaystyle \Box }
Es sind
(
d
r
)
3
⋅
∑
ν
=
1
r
1
ν
3
{\displaystyle (d_{r})^{3}\cdot \sum _{\nu =1}^{r}{\frac {1}{\nu ^{3}}}}
und
(
d
r
)
3
⋅
1
r
−
s
∑
ν
=
s
+
1
r
1
ν
2
{\displaystyle (d_{r})^{3}\cdot {\frac {1}{r-s}}\sum _{\nu =s+1}^{r}{\frac {1}{\nu ^{2}}}}
für
r
>
s
{\displaystyle r>s}
jeweils ganze Zahlen.
Ist
ν
≤
r
{\displaystyle \nu \leq r}
, so ist
ν
=
∏
p
≤
r
p
α
(
p
,
ν
)
{\displaystyle \nu =\prod _{p\leq r}p^{\alpha (p,\nu )}}
mit
α
(
p
,
ν
)
≤
⌊
ln
r
ln
p
⌋
{\displaystyle \alpha (p,\nu )\leq \lfloor {\frac {\ln r}{\ln p}}\rfloor }
. Dann ist
ν
3
=
∏
p
≤
r
p
3
α
(
p
,
ν
)
|
∏
p
≤
r
p
3
⌊
ln
r
ln
p
⌋
=
(
∏
p
≤
r
p
⌊
ln
r
ln
p
⌋
)
3
=
(
d
r
)
3
,
{\displaystyle \nu ^{3}=\prod _{p\leq r}p^{3\alpha (p,\nu )}\ \left|\ \prod _{p\leq r}\right.p^{3\lfloor {\frac {\ln r}{\ln p}}\rfloor }=\left(\prod _{p\leq r}p^{\lfloor {\frac {\ln r}{\ln p}}\rfloor }\right)^{3}=(d_{r})^{3}\ ,}
also
(
d
r
)
3
⋅
∑
ν
=
1
r
1
ν
3
∈
Z
{\displaystyle (d_{r})^{3}\cdot \sum _{\nu =1}^{r}{\frac {1}{\nu ^{3}}}\in \mathbb {Z} }
. Sei nun
r
>
s
{\displaystyle r>s}
. Dann ist
r
−
s
=
∏
p
≤
r
p
α
(
p
,
r
−
s
)
{\displaystyle r-s=\prod _{p\leq r}p^{\alpha (p,r-s)}}
mit
α
(
p
,
r
−
s
)
≤
⌊
ln
r
ln
p
⌋
{\displaystyle \alpha (p,r-s)\leq \lfloor {\frac {\ln r}{\ln p}}\rfloor }
. Dann folgt
(
r
−
s
)
ν
2
=
∏
p
≤
r
p
2
α
(
p
,
ν
)
+
α
(
p
,
r
−
s
)
|
∏
p
≤
r
p
3
⌊
ln
r
ln
p
⌋
=
(
∏
p
≤
r
p
⌊
ln
r
ln
p
⌋
)
3
=
(
d
r
)
3
,
{\displaystyle (r-s)\nu ^{2}=\prod _{p\leq r}p^{2\alpha (p,\nu )+\alpha (p,r-s)}\ \left|\ \prod _{p\leq r}\right.p^{3\lfloor {\frac {\ln r}{\ln p}}\rfloor }=\left(\prod _{p\leq r}p^{\lfloor {\frac {\ln r}{\ln p}}\rfloor }\right)^{3}=(d_{r})^{3}\ ,}
also
(
d
r
)
3
⋅
1
r
−
s
∑
ν
=
s
+
1
r
1
ν
2
∈
Z
{\displaystyle (d_{r})^{3}\cdot {\frac {1}{r-s}}\sum _{\nu =s+1}^{r}{\frac {1}{\nu ^{2}}}\in \mathbb {Z} }
.
◻
{\displaystyle \Box }
Es gilt
0
<
I
n
≤
2
ζ
(
3
)
(
2
−
1
)
4
n
.
{\displaystyle 0<I_{n}\leq 2\zeta (3)({\sqrt {2}}-1)^{4n}\ .}
Es ist
−
ln
x
y
=
∫
x
y
1
1
u
d
u
{\displaystyle -\ln xy=\int _{xy}^{1}{\frac {1}{u}}{\rm {d}}u}
. Mit der Transformation
u
=
1
−
(
1
−
x
y
)
v
{\displaystyle u=1-(1-xy)v}
folgt
−
ln
x
y
=
∫
0
1
d
v
1
−
(
1
−
x
y
)
v
(
1
−
x
y
)
{\displaystyle -\ln xy=\int _{0}^{1}{\frac {{\rm {d}}v}{1-(1-xy)v}}(1-xy)}
, also
−
ln
x
y
1
−
x
y
=
∫
0
1
d
v
1
−
(
1
−
x
y
)
v
{\displaystyle -{\frac {\ln xy}{1-xy}}=\int _{0}^{1}{\frac {{\rm {d}}v}{1-(1-xy)v}}}
. Es folgt
I
n
=
∫
0
1
P
n
(
y
)
∫
0
1
(
∫
0
1
1
n
!
d
n
d
x
n
[
x
n
(
1
−
x
)
n
]
1
1
−
(
1
−
x
y
)
v
d
x
)
d
v
d
y
.
{\displaystyle I_{n}=\int _{0}^{1}P_{n}(y)\int _{0}^{1}\left(\int _{0}^{1}{\frac {1}{n!}}{\frac {{\rm {d}}^{n}}{{\rm {d}}x^{n}}}[x^{n}(1-x)^{n}]{\frac {1}{1-(1-xy)v}}{\rm {d}}x\right){\rm {d}}v{\rm {d}}y\ .}
Im inneren Integral wende man partielle Integration an. Da alle Ableitungen
d
ν
d
x
ν
[
x
n
(
1
−
x
)
n
]
{\displaystyle {\frac {{\rm {d}}^{\nu }}{{\rm {d}}x^{\nu }}}[x^{n}(1-x)^{n}]}
für
ν
<
n
{\displaystyle \nu <n}
in
x
=
0
{\displaystyle x=0}
und
x
=
1
{\displaystyle x=1}
verschwinden, treten bei
n
{\displaystyle n}
-maliger partieller Integration keine Randterme auf. Wir bekommen somit
∫
0
1
1
n
!
d
n
d
x
n
[
x
n
(
1
−
x
)
n
]
1
1
−
v
+
x
y
v
d
x
=
(
−
1
)
n
∫
0
1
1
n
!
x
n
(
1
−
x
)
n
∂
n
∂
x
n
[
1
1
−
v
+
x
y
v
]
d
x
=
∫
0
1
x
n
(
1
−
x
)
n
1
(
1
−
v
+
x
y
v
)
n
+
1
(
y
v
)
n
d
x
.
{\displaystyle {\begin{array}{lll}\int _{0}^{1}{\frac {1}{n!}}{\frac {{\rm {d}}^{n}}{{\rm {d}}x^{n}}}[x^{n}(1-x)^{n}]{\frac {1}{1-v+xyv}}{\rm {d}}x&=&(-1)^{n}\int _{0}^{1}{\frac {1}{n!}}x^{n}(1-x)^{n}{\frac {\partial ^{n}}{\partial x^{n}}}[{\frac {1}{1-v+xyv}}]{\rm {d}}x\\&=&\int _{0}^{1}x^{n}(1-x)^{n}{\frac {1}{(1-v+xyv)^{n+1}}}(yv)^{n}{\rm {d}}x\ .\\\end{array}}}
Daraus folgt
I
n
=
∫
0
1
P
n
(
y
)
y
n
∫
0
1
x
n
(
1
−
x
)
n
∫
0
1
v
n
(
1
−
(
1
−
x
y
)
v
)
n
+
1
d
v
d
x
d
y
.
{\displaystyle I_{n}=\int _{0}^{1}P_{n}(y)y^{n}\int _{0}^{1}x^{n}(1-x)^{n}\int _{0}^{1}{\frac {v^{n}}{(1-(1-xy)v)^{n+1}}}{\rm {d}}v{\rm {d}}x{\rm {d}}y\ .}
Nun substituieren wir
v
=
1
−
z
1
−
z
(
1
−
x
y
)
{\displaystyle v={\frac {1-z}{1-z(1-xy)}}}
. Dann ist
d
v
d
z
=
−
x
y
(
1
−
z
(
1
−
x
y
)
)
2
<
0
{\displaystyle {\frac {{\rm {d}}v}{{\rm {d}}z}}=-{\frac {xy}{(1-z(1-xy))^{2}}}<0}
, also haben wir eine gültige Variablentransformation. Wegen
1
−
(
1
−
x
y
)
v
=
x
y
1
−
z
(
1
−
x
y
)
{\displaystyle 1-(1-xy)v={\frac {xy}{1-z(1-xy)}}}
ist
I
n
=
∫
0
1
P
n
(
y
)
y
n
∫
0
1
x
n
(
1
−
x
)
n
∫
0
1
(
1
−
z
)
n
(
1
−
z
(
1
−
x
y
)
)
n
(
1
−
z
(
1
−
x
y
)
)
n
+
1
(
x
y
)
n
+
1
x
y
(
1
−
z
(
1
−
x
y
)
)
2
d
z
d
x
d
y
=
∫
0
1
P
n
(
y
)
∫
0
1
(
1
−
x
)
n
∫
0
1
(
1
−
z
)
n
1
1
−
z
(
1
−
x
y
)
d
z
d
x
d
y
=
∫
0
1
(
1
−
z
)
n
∫
0
1
(
1
−
x
)
n
∫
0
1
1
n
!
d
n
d
y
n
[
y
n
(
1
−
y
)
n
]
1
1
−
z
+
x
y
z
d
y
d
x
d
z
.
{\displaystyle {\begin{array}{lll}I_{n}&=&\int _{0}^{1}P_{n}(y)y^{n}\int _{0}^{1}x^{n}(1-x)^{n}\int _{0}^{1}{\frac {(1-z)^{n}}{(1-z(1-xy))^{n}}}{\frac {(1-z(1-xy))^{n+1}}{(xy)^{n+1}}}{\frac {xy}{(1-z(1-xy))^{2}}}{\rm {d}}z{\rm {d}}x{\rm {d}}y\\&=&\int _{0}^{1}P_{n}(y)\int _{0}^{1}(1-x)^{n}\int _{0}^{1}(1-z)^{n}{\frac {1}{1-z(1-xy)}}{\rm {d}}z{\rm {d}}x{\rm {d}}y\\&=&\int _{0}^{1}(1-z)^{n}\int _{0}^{1}(1-x)^{n}\int _{0}^{1}{\frac {1}{n!}}{\frac {d^{n}}{dy^{n}}}[y^{n}(1-y)^{n}]{\frac {1}{1-z+xyz}}{\rm {d}}y{\rm {d}}x{\rm {d}}z\ .\end{array}}}
Wir führen wiederum
n
{\displaystyle n}
-malige partielle Integration (ebenfalls ohne Randterme) durch und erhalten
I
n
=
∫
0
1
(
1
−
z
)
n
∫
0
1
(
1
−
x
)
n
(
−
1
)
n
∫
0
1
1
n
!
y
n
(
1
−
y
)
n
∂
n
∂
y
n
[
1
1
−
z
+
x
y
z
]
d
y
d
x
d
z
=
∫
0
1
(
1
−
z
)
n
∫
0
1
(
1
−
x
)
n
∫
0
1
y
n
(
1
−
y
)
n
(
x
z
)
n
1
(
1
−
z
+
x
y
z
)
n
+
1
d
y
d
x
d
z
=
∫
0
1
∫
0
1
∫
0
1
[
F
(
x
,
y
,
z
)
n
]
1
−
(
1
−
x
y
)
z
d
x
d
y
d
z
{\displaystyle {\begin{array}{lll}I_{n}&=&\int _{0}^{1}(1-z)^{n}\int _{0}^{1}(1-x)^{n}(-1)^{n}\int _{0}^{1}{\frac {1}{n!}}y^{n}(1-y)^{n}{\frac {\partial ^{n}}{\partial y^{n}}}[{\frac {1}{1-z+xyz}}]{\rm {d}}y{\rm {d}}x{\rm {d}}z\\&=&\int _{0}^{1}(1-z)^{n}\int _{0}^{1}(1-x)^{n}\int _{0}^{1}y^{n}(1-y)^{n}(xz)^{n}{\frac {1}{(1-z+xyz)^{n+1}}}{\rm {d}}y{\rm {d}}x{\rm {d}}z\\&=&\int _{0}^{1}\int _{0}^{1}\int _{0}^{1}{\frac {[F(x,y,z)^{n}]}{1-(1-xy)z}}{\rm {d}}x{\rm {d}}y{\rm {d}}z\\\end{array}}}
mit
F
(
x
,
y
,
z
)
:=
x
(
1
−
x
)
y
(
1
−
y
)
z
(
1
−
z
)
1
−
(
1
−
x
y
)
z
{\displaystyle F(x,y,z):={\frac {x(1-x)y(1-y)z(1-z)}{1-(1-xy)z}}}
. Es ist
I
n
>
0
{\displaystyle I_{n}>0}
, da der Integrand echt positiv auf
(
0
,
1
)
3
{\displaystyle (0,1)^{3}}
ist. Das folgende Lemma Nummer 8 besagt
F
(
x
,
y
,
z
)
≤
(
2
−
1
)
4
{\displaystyle F(x,y,z)\leq ({\sqrt {2}}-1)^{4}}
auf
(
0
,
1
)
3
{\displaystyle (0,1)^{3}}
. Nutzen wir diese Abschätzung vorab aus, so folgt
0
<
I
n
≤
(
2
−
1
)
4
n
∫
0
1
∫
0
1
∫
0
1
1
1
−
(
1
−
x
y
)
z
d
z
d
x
d
y
=
(
2
−
1
)
4
n
∫
0
1
∫
0
1
−
ln
x
y
1
−
x
y
d
x
d
y
nach Rechnung zu Beginn des Beweises
=
2
(
2
−
1
)
4
n
ζ
(
3
)
nach Lemma Nummer 3 mit
r
=
s
=
0
,
{\displaystyle {\begin{array}{lll}0<I_{n}&\leq &({\sqrt {2}}-1)^{4n}\int _{0}^{1}\int _{0}^{1}\int _{0}^{1}{\frac {1}{1-(1-xy)z}}{\rm {d}}z{\rm {d}}x{\rm {d}}y\\&=&({\sqrt {2}}-1)^{4n}\int _{0}^{1}\int _{0}^{1}-{\frac {\ln xy}{1-xy}}{\rm {d}}x{\rm {d}}y{\text{ nach Rechnung zu Beginn des Beweises}}\\&=&2({\sqrt {2}}-1)^{4n}\zeta (3){\text{ nach Lemma Nummer 3 mit }}r=s=0\ ,\end{array}}}
was zu beweisen war.
◻
{\displaystyle \Box }
Es gilt
F
(
x
,
y
,
z
)
:=
x
(
1
−
x
)
y
(
1
−
y
)
z
(
1
−
z
)
1
−
(
1
−
x
y
)
z
≤
(
2
−
1
)
4
{\displaystyle F(x,y,z):={\frac {x(1-x)y(1-y)z(1-z)}{1-(1-xy)z}}\leq ({\sqrt {2}}-1)^{4}}
für alle
(
x
,
y
,
z
)
∈
(
0
,
1
)
3
{\displaystyle (x,y,z)\in (0,1)^{3}}
.
Zunächst zeigen wir, dass
F
{\displaystyle F}
auf den Rand von
(
0
,
1
)
3
{\displaystyle (0,1)^{3}}
durch null stetig fortgesetzt werden kann. Offenbar müssen wir nur den Bereich mit
z
=
1
{\displaystyle z=1}
untersuchen. Ist nun
z
n
→
1
{\displaystyle z_{n}\rightarrow 1}
, so folgt
|
F
(
x
n
,
y
n
,
z
n
)
|
=
|
1
−
x
n
|
|
1
−
y
n
|
|
z
n
|
|
x
n
y
n
(
1
−
z
n
)
1
−
(
1
−
x
n
y
n
)
z
n
|
≤
|
x
n
y
n
(
1
−
z
n
)
x
n
y
n
+
(
1
−
z
n
)
−
x
n
y
n
(
1
−
z
n
)
|
=
1
|
1
1
−
z
n
+
1
x
n
y
n
−
1
|
≤
1
|
1
1
−
z
n
−
1
|
=
1
−
z
n
1
−
(
1
−
z
n
)
→
0
.
{\displaystyle {\begin{array}{lll}|F(x_{n},y_{n},z_{n})|&=&|1-x_{n}||1-y_{n}||z_{n}|\left|{\frac {x_{n}y_{n}(1-z_{n})}{1-(1-x_{n}y_{n})z_{n}}}\right|\leq \left|{\frac {x_{n}y_{n}(1-z_{n})}{x_{n}y_{n}+(1-z_{n})-x_{n}y_{n}(1-z_{n})}}\right|\\&=&{\frac {1}{\left|{\frac {1}{1-z_{n}}}+{\frac {1}{x_{n}y_{n}}}-1\right|}}\leq {\frac {1}{\left|{\frac {1}{1-z_{n}}}-1\right|}}={\frac {1-z_{n}}{1-(1-z_{n})}}\rightarrow 0\ .\\\end{array}}}
Also besitzt
f
{\displaystyle f}
ein globales Maximum
(
x
0
,
y
0
,
z
0
)
{\displaystyle (x_{0},y_{0},z_{0})}
auf
[
0
,
1
]
3
{\displaystyle [0,1]^{3}}
, welches im Inneren liegt und somit ein kritischer Punkt von
F
{\displaystyle F}
ist.
Wir haben
∂
∂
x
F
(
x
,
y
,
z
)
=
[
1
−
(
1
−
x
y
)
z
]
(
1
−
2
x
)
(
y
−
y
2
)
(
z
−
z
2
)
−
y
z
(
x
−
x
2
)
(
y
−
y
2
)
(
z
−
z
2
)
[
1
−
(
1
−
x
y
)
z
]
2
{\displaystyle {\frac {\partial }{\partial x}}F(x,y,z)={\frac {[1-(1-xy)z](1-2x)(y-y^{2})(z-z^{2})-yz(x-x^{2})(y-y^{2})(z-z^{2})}{[1-(1-xy)z]^{2}}}}
, also
0
=
[
1
−
(
1
−
x
0
y
0
)
z
0
]
(
1
−
2
x
0
)
−
y
0
z
0
(
x
0
−
x
0
2
)
{\displaystyle 0=[1-(1-x_{0}y_{0})z_{0}](1-2x_{0})-y_{0}z_{0}(x_{0}-x_{0}^{2})}
und analog mit
∂
∂
y
F
(
x
,
y
,
z
)
=
[
1
−
(
1
−
x
y
)
z
]
(
x
−
x
2
)
(
1
−
2
y
)
(
z
−
z
2
)
−
x
z
(
x
−
x
2
)
(
y
−
y
2
)
(
z
−
z
2
)
[
1
−
(
1
−
x
y
)
z
]
2
{\displaystyle {\frac {\partial }{\partial y}}F(x,y,z)={\frac {[1-(1-xy)z](x-x^{2})(1-2y)(z-z^{2})-xz(x-x^{2})(y-y^{2})(z-z^{2})}{[1-(1-xy)z]^{2}}}}
0
=
[
1
−
(
1
−
x
0
y
0
)
z
0
]
(
1
−
2
y
0
)
−
x
0
z
0
(
y
0
−
y
0
2
)
.
{\displaystyle 0=[1-(1-x_{0}y_{0})z_{0}](1-2y_{0})-x_{0}z_{0}(y_{0}-y_{0}^{2})\ .}
Subtrahieren der Gleichungen impliziert
0
=
[
1
−
(
1
−
x
0
y
0
)
z
0
]
(
2
y
0
−
2
x
0
)
−
z
0
[
(
x
0
−
x
0
2
)
y
0
−
(
y
0
−
y
0
2
)
x
0
]
=
(
y
0
−
x
0
)
(
2
−
2
z
0
+
x
0
y
0
z
0
)
{\displaystyle 0=[1-(1-x_{0}y_{0})z_{0}](2y_{0}-2x_{0})-z_{0}[(x_{0}-x_{0}^{2})y_{0}-(y_{0}-y_{0}^{2})x_{0}]=(y_{0}-x_{0})(2-2z_{0}+x_{0}y_{0}z_{0})}
, also
x
0
=
y
0
{\displaystyle x_{0}=y_{0}}
, weil der zweite
Faktor echt positiv ist. Weiter ist
∂
∂
z
F
(
x
,
x
,
z
)
=
[
1
−
(
1
−
x
2
)
z
]
(
x
−
x
2
)
2
(
1
−
2
z
)
+
(
1
−
x
2
)
(
x
−
x
2
)
2
(
z
−
z
2
)
[
1
−
(
1
−
x
2
)
z
]
2
{\displaystyle {\frac {\partial }{\partial z}}F(x,x,z)={\frac {[1-(1-x^{2})z](x-x^{2})^{2}(1-2z)+(1-x^{2})(x-x^{2})^{2}(z-z^{2})}{[1-(1-x^{2})z]^{2}}}}
, und man erhält
0
=
1
−
2
z
0
+
z
0
2
−
x
0
2
z
0
2
.
{\displaystyle 0=1-2z_{0}+z_{0}^{2}-x_{0}^{2}z_{0}^{2}\ .}
Zieht man das
z
0
{\displaystyle z_{0}}
-fache der ersten Gleichung vom
x
0
{\displaystyle x_{0}}
-fachen der dritten Gleichung ab, liefert dies
0
=
x
0
−
z
0
+
z
0
2
−
x
0
z
0
2
{\displaystyle 0=x_{0}-z_{0}+z_{0}^{2}-x_{0}z_{0}^{2}}
, also
x
0
=
z
0
2
−
z
0
z
0
2
−
1
=
1
−
1
z
0
+
1
{\displaystyle x_{0}={\frac {z_{0}^{2}-z_{0}}{z_{0}^{2}-1}}=1-{\frac {1}{z_{0}+1}}}
. Einsetzen in das
(
z
0
+
1
)
2
{\displaystyle (z_{0}+1)^{2}}
-fache der dritten Gleichung impliziert
0
=
−
2
z
0
2
+
1
{\displaystyle 0=-2z_{0}^{2}+1}
, also
z
0
=
1
2
{\displaystyle z_{0}={\frac {1}{\sqrt {2}}}}
und
x
0
=
y
0
=
2
−
1
{\displaystyle x_{0}=y_{0}={\sqrt {2}}-1}
. Aus
F
(
x
0
,
y
0
,
z
0
)
=
F
(
2
−
1
,
2
−
1
,
1
2
)
=
(
2
−
1
)
4
{\displaystyle F(x_{0},y_{0},z_{0})=F\left({\sqrt {2}}-1,{\sqrt {2}}-1,{\tfrac {1}{\sqrt {2}}}\right)=({\sqrt {2}}-1)^{4}}
folgt die Behauptung.
◻
{\displaystyle \Box }
ζ
(
3
)
{\displaystyle \zeta (3)}
ist irrational.
Angenommen, es wäre
ζ
(
3
)
=
p
q
{\displaystyle \zeta (3)={\frac {p}{q}}}
mit
p
,
q
∈
N
{\displaystyle p,q\in \mathbb {N} }
und
ggT
(
p
,
q
)
=
1
{\displaystyle \operatorname {ggT} (p,q)=1}
. Nach dem Euklidischen Algorithmus gibt es wegen
ggT
(
p
,
q
)
=
1
{\displaystyle \operatorname {ggT} (p,q)=1}
Zahlen
m
0
,
n
0
∈
Z
{\displaystyle m_{0},n_{0}\in \mathbb {Z} }
mit
m
0
p
+
n
0
q
=
1
{\displaystyle m_{0}p+n_{0}q=1}
. Dann ist
inf
{
|
m
q
−
n
p
|
|
m
,
n
∈
Z
,
m
q
−
n
p
≠
0
}
=
1
{\displaystyle \inf\{|mq-np|\ |\ m,n\in \mathbb {Z} ,\ mq-np\neq 0\}=1}
, also
inf
{
|
m
−
n
ζ
(
3
)
|
|
m
,
n
∈
Z
,
m
−
n
ζ
(
3
)
≠
0
}
=
1
q
.
{\displaystyle \inf\{|m-n\zeta (3)|\ |\ m,n\in \mathbb {Z} ,\ m-n\zeta (3)\neq 0\}={\frac {1}{q}}\ .}
Nach Lemma Nummer 4, Korollar Nummer 6 und Lemma Nummer 7 ist somit
1
q
≤
2
ζ
(
3
)
(
2
−
1
)
4
n
(
d
n
)
3
≤
2
ζ
(
3
)
(
2
−
1
)
4
n
n
3
π
(
n
)
=
2
ζ
(
3
)
e
4
n
ln
(
2
−
1
)
+
3
π
(
n
)
ln
n
.
{\displaystyle {\frac {1}{q}}\leq 2\zeta (3)({\sqrt {2}}-1)^{4n}(d_{n})^{3}\leq 2\zeta (3)({\sqrt {2}}-1)^{4n}n^{3\pi (n)}=2\zeta (3)e^{4n\ln({\sqrt {2}}-1)+3\pi (n)\ln n}\ .}
Dies ist äquivalent zu
π
(
n
)
ln
n
n
≥
1
3
n
ln
(
1
2
q
ζ
(
3
)
)
−
4
3
ln
(
2
−
1
)
{\displaystyle \pi (n){\frac {\ln n}{n}}\geq {\frac {1}{3n}}\ln \left({\frac {1}{2q\zeta (3)}}\right)-{\frac {4}{3}}\ln \left({\sqrt {2}}-1\right)}
. Übergang zum Grenzwert impliziert gemäß dem Primzahlsatz
1
≥
−
4
3
ln
(
2
−
1
)
{\displaystyle 1\geq -{\frac {4}{3}}\ln({\sqrt {2}}-1)}
, aber es ist
−
4
3
ln
(
2
−
1
)
>
1
{\displaystyle -{\frac {4}{3}}\ln({\sqrt {2}}-1)>1}
.
◻
{\displaystyle \Box }
Roger Apéry: Irrationalité de
ζ
(
2
)
{\displaystyle \zeta (2)}
et
ζ
(
3
)
{\displaystyle \zeta (3)}
. Astérisque, 61 , 11-13, 1979.
Frits Beukers: A note on the irrationality of
ζ
(
2
)
{\displaystyle \zeta (2)}
and
ζ
(
3
)
{\displaystyle \zeta (3)}
. Bulletin of the London Mathematical Society, 11 , 268-272, 1979.