Produto infinito para pi
Comparação da convergência do produto Wallis (asteriscos roxos) e várias séries históricas infinitas para
π .
S n é a aproximação após tomar
n termos. Cada subtrama subseqüente amplia a área sombreada horizontalmente em 10 vezes.
(clique para detalhes)
Em matemática , o produto Wallis para π , publicado em 1656 por John Wallis , afirma que
π
2
=
∏
n
=
1
∞
4
n
2
4
n
2
-
1
=
∏
n
=
1
∞
(
2
n
2
n
-
1
⋅
2
n
2
n
+
1
)
=
(
2
1
⋅
2
3
)
⋅
(
4
3
⋅
4
5
)
⋅
(
6
5
⋅
6
7
)
⋅
(
8
7
⋅
8
9
)
⋅
⋯
{\ displaystyle {\ begin {align} {\ frac {\ pi} {2}} & = \ prod _ {n = 1} ^ {\ infty} {\ frac {4n ^ {2}} {4n ^ {2 } -1}} = \ prod _ {n = 1} ^ {\ infty} \ left ({\ frac {2n} {2n-1}} \ cdot {\ frac {2n} {2n + 1}} \ right ) \\ [6pt] & = {\ Grande (} {\ frac {2} {1}} \ cdot {\ frac {2} {3}} {\ Grande)} \ cdot {\ Grande (} {\ frac {4} {3}} \ cdot {\ frac {4} {5}} {\ Grande)} \ cdot {\ Grande (} {\ frac {6} {5}} \ cdot {\ frac {6} { 7}} {\ Big)} \ cdot {\ Big (} {\ frac {8} {7}} \ cdot {\ frac {8} {9}} {\ Big)} \ cdot \; \ cdots \\ \ end {alinhado}}}
Prova usando integração
Wallis derivou esse produto infinito como é feito nos livros de cálculo hoje, examinando os valores pares e ímpares de , e observando que, para grande , o aumento de 1 resulta em uma mudança que se torna cada vez menor à medida que aumenta. Deixar
∫
0
π
pecado
n
x
d
x
{\ displaystyle \ int _ {0} ^ {\ pi} \ sin ^ {n} x \, dx}
n
{\ displaystyle n}
n
{\ displaystyle n}
n
{\ displaystyle n}
n
{\ displaystyle n}
eu
(
n
)
=
∫
0
π
pecado
n
x
d
x
.
{\ displaystyle I (n) = \ int _ {0} ^ {\ pi} \ sin ^ {n} x \, dx.}
(Esta é uma forma de integrais de Wallis .) Integrar por partes :
você
=
pecado
n
-
1
x
⇒
d
você
=
(
n
-
1
)
pecado
n
-
2
x
cos
x
d
x
d
v
=
pecado
x
d
x
⇒
v
=
-
cos
x
{\ displaystyle {\ begin {alinhado} u & = \ sin ^ {n-1} x \\\ Rightarrow du & = (n-1) \ sin ^ {n-2} x \ cos x \, dx \\ dv & = \ sin x \, dx \\\ Rightarrow v & = - \ cos x \ end {alinhado}}}
⇒
eu
(
n
)
=
∫
0
π
pecado
n
x
d
x
=
-
pecado
n
-
1
x
cos
x
|
0
π
-
∫
0
π
(
-
cos
x
)
(
n
-
1
)
pecado
n
-
2
x
cos
x
d
x
=
0
+
(
n
-
1
)
∫
0
π
cos
2
x
pecado
n
-
2
x
d
x
,
n
>
1
=
(
n
-
1
)
∫
0
π
(
1
-
pecado
2
x
)
pecado
n
-
2
x
d
x
=
(
n
-
1
)
∫
0
π
pecado
n
-
2
x
d
x
-
(
n
-
1
)
∫
0
π
pecado
n
x
d
x
=
(
n
-
1
)
eu
(
n
-
2
)
-
(
n
-
1
)
eu
(
n
)
=
n
-
1
n
eu
(
n
-
2
)
⇒
eu
(
n
)
eu
(
n
-
2
)
=
n
-
1
n
{\ displaystyle {\ begin {alinhado} \ Rightarrow I (n) & = \ int _ {0} ^ {\ pi} \ sin ^ {n} x \, dx \\ [6pt] {} & = - \ sin ^ {n-1} x \ cos x {\ Biggl |} _ {0} ^ {\ pi} - \ int _ {0} ^ {\ pi} (- \ cos x) (n-1) \ sin ^ {n-2} x \ cos x \, dx \\ [6pt] {} & = 0+ (n-1) \ int _ {0} ^ {\ pi} \ cos ^ {2} x \ sin ^ { n-2} x \, dx, \ qquad n> 1 \\ [6pt] {} & = (n-1) \ int _ {0} ^ {\ pi} (1- \ sin ^ {2} x) \ sin ^ {n-2} x \, dx \\ [6pt] {} & = (n-1) \ int _ {0} ^ {\ pi} \ sin ^ {n-2} x \, dx- (n-1) \ int _ {0} ^ {\ pi} \ sin ^ {n} x \, dx \\ [6pt] {} & = (n-1) I (n-2) - (n- 1) I (n) \\ [6pt] {} & = {\ frac {n-1} {n}} I (n-2) \\ [6pt] \ Rightarrow {\ frac {I (n)} { I (n-2)}} & = {\ frac {n-1} {n}} \\ [6pt] \ end {alinhado}}}
Agora, fazemos duas substituições de variáveis por conveniência para obter:
eu
(
2
n
)
=
2
n
-
1
2
n
eu
(
2
n
-
2
)
{\ displaystyle I (2n) = {\ frac {2n-1} {2n}} I (2n-2)}
eu
(
2
n
+
1
)
=
2
n
2
n
+
1
eu
(
2
n
-
1
)
{\ displaystyle I (2n + 1) = {\ frac {2n} {2n + 1}} I (2n-1)}
Obtemos valores para e para uso posterior.
eu
(
0
)
{\ displaystyle I (0)}
eu
(
1
)
{\ displaystyle I (1)}
eu
(
0
)
=
∫
0
π
d
x
=
x
|
0
π
=
π
eu
(
1
)
=
∫
0
π
pecado
x
d
x
=
-
cos
x
|
0
π
=
(
-
cos
π
)
-
(
-
cos
0
)
=
-
(
-
1
)
-
(
-
1
)
=
2
{\ displaystyle {\ begin {alinhados} I (0) & = \ int _ {0} ^ {\ pi} dx = x {\ Biggl |} _ {0} ^ {\ pi} = \ pi \\ [6pt ] I (1) & = \ int _ {0} ^ {\ pi} \ sin x \, dx = - \ cos x {\ Biggl |} _ {0} ^ {\ pi} = (- \ cos \ pi ) - (- \ cos 0) = - (- 1) - (- 1) = 2 \\ [6pt] \ end {alinhado}}}
Agora, calculamos para valores pares aplicando repetidamente o resultado da relação de recorrência da integração por partes. Eventualmente, acabamos chegando , o que calculamos.
eu
(
2
n
)
{\ displaystyle I (2n)}
eu
(
0
)
{\ displaystyle I (0)}
eu
(
2
n
)
=
∫
0
π
pecado
2
n
x
d
x
=
2
n
-
1
2
n
eu
(
2
n
-
2
)
=
2
n
-
1
2
n
⋅
2
n
-
3
2
n
-
2
eu
(
2
n
-
4
)
{\ displaystyle I (2n) = \ int _ {0} ^ {\ pi} \ sin ^ {2n} x \, dx = {\ frac {2n-1} {2n}} I (2n-2) = { \ frac {2n-1} {2n}} \ cdot {\ frac {2n-3} {2n-2}} I (2n-4)}
=
2
n
-
1
2
n
⋅
2
n
-
3
2
n
-
2
⋅
2
n
-
5
2
n
-
4
⋅
⋯
⋅
5
6
⋅
3
4
⋅
1
2
eu
(
0
)
=
π
∏
k
=
1
n
2
k
-
1
2
k
{\ displaystyle = {\ frac {2n-1} {2n}} \ cdot {\ frac {2n-3} {2n-2}} \ cdot {\ frac {2n-5} {2n-4}} \ cdot \ cdots \ cdot {\ frac {5} {6}} \ cdot {\ frac {3} {4}} \ cdot {\ frac {1} {2}} I (0) = \ pi \ prod _ {k = 1} ^ {n} {\ frac {2k-1} {2k}}}
Repetindo o processo para valores ímpares ,
eu
(
2
n
+
1
)
{\ displaystyle I (2n + 1)}
eu
(
2
n
+
1
)
=
∫
0
π
pecado
2
n
+
1
x
d
x
=
2
n
2
n
+
1
eu
(
2
n
-
1
)
=
2
n
2
n
+
1
⋅
2
n
-
2
2
n
-
1
eu
(
2
n
-
3
)
{\ displaystyle I (2n + 1) = \ int _ {0} ^ {\ pi} \ sin ^ {2n + 1} x \, dx = {\ frac {2n} {2n + 1}} I (2n- 1) = {\ frac {2n} {2n + 1}} \ cdot {\ frac {2n-2} {2n-1}} I (2n-3)}
=
2
n
2
n
+
1
⋅
2
n
-
2
2
n
-
1
⋅
2
n
-
4
2
n
-
3
⋅
⋯
⋅
6
7
⋅
4
5
⋅
2
3
eu
(
1
)
=
2
∏
k
=
1
n
2
k
2
k
+
1
{\ displaystyle = {\ frac {2n} {2n + 1}} \ cdot {\ frac {2n-2} {2n-1}} \ cdot {\ frac {2n-4} {2n-3}} \ cdot \ cdots \ cdot {\ frac {6} {7}} \ cdot {\ frac {4} {5}} \ cdot {\ frac {2} {3}} I (1) = 2 \ prod _ {k = 1} ^ {n} {\ frac {2k} {2k + 1}}}
Fazemos a seguinte observação, com base no fato de que
pecado
x
≤
1
{\ displaystyle \ sin {x} \ leq 1}
pecado
2
n
+
1
x
≤
pecado
2
n
x
≤
pecado
2
n
-
1
x
,
0
≤
x
≤
π
{\ displaystyle \ sin ^ {2n + 1} x \ leq \ sin ^ {2n} x \ leq \ sin ^ {2n-1} x, 0 \ leq x \ leq \ pi}
⇒
eu
(
2
n
+
1
)
≤
eu
(
2
n
)
≤
eu
(
2
n
-
1
)
{\ displaystyle \ Rightarrow I (2n + 1) \ leq I (2n) \ leq I (2n-1)}
Dividindo por :
eu
(
2
n
+
1
)
{\ displaystyle I (2n + 1)}
⇒
1
≤
eu
(
2
n
)
eu
(
2
n
+
1
)
≤
eu
(
2
n
-
1
)
eu
(
2
n
+
1
)
=
2
n
+
1
2
n
{\ displaystyle \ Rightarrow 1 \ leq {\ frac {I (2n)} {I (2n + 1)}} \ leq {\ frac {I (2n-1)} {I (2n + 1)}} = { \ frac {2n + 1} {2n}}}
, onde a igualdade vem de nossa relação de recorrência.
Pelo teorema de compressão ,
⇒
lim
n
→
∞
eu
(
2
n
)
eu
(
2
n
+
1
)
=
1
{\ displaystyle \ Rightarrow \ lim _ {n \ rightarrow \ infty} {\ frac {I (2n)} {I (2n + 1)}} = 1}
lim
n
→
∞
eu
(
2
n
)
eu
(
2
n
+
1
)
=
π
2
lim
n
→
∞
∏
k
=
1
n
(
2
k
-
1
2
k
⋅
2
k
+
1
2
k
)
=
1
{\ displaystyle \ lim _ {n \ rightarrow \ infty} {\ frac {I (2n)} {I (2n + 1)}} = {\ frac {\ pi} {2}} \ lim _ {n \ rightarrow \ infty} \ prod _ {k = 1} ^ {n} \ left ({\ frac {2k-1} {2k}} \ cdot {\ frac {2k + 1} {2k}} \ right) = 1}
⇒
π
2
=
∏
k
=
1
∞
(
2
k
2
k
-
1
⋅
2
k
2
k
+
1
)
=
2
1
⋅
2
3
⋅
4
3
⋅
4
5
⋅
6
5
⋅
6
7
⋅
⋯
{\ displaystyle \ Rightarrow {\ frac {\ pi} {2}} = \ prod _ {k = 1} ^ {\ infty} \ left ({\ frac {2k} {2k-1}} \ cdot {\ frac {2k} {2k + 1}} \ right) = {\ frac {2} {1}} \ cdot {\ frac {2} {3}} \ cdot {\ frac {4} {3}} \ cdot { \ frac {4} {5}} \ cdot {\ frac {6} {5}} \ cdot {\ frac {6} {7}} \ cdot \ cdots}
Prova usando o produto infinito de Euler para a função seno
Embora a prova acima seja tipicamente apresentada em livros modernos de cálculo, o produto de Wallis é, em retrospecto, um corolário fácil do produto infinito de Euler posterior para a função seno .
pecado
x
x
=
∏
n
=
1
∞
(
1
-
x
2
n
2
π
2
)
{\ displaystyle {\ frac {\ sin x} {x}} = \ prod _ {n = 1} ^ {\ infty} \ left (1 - {\ frac {x ^ {2}} {n ^ {2} \ pi ^ {2}}} \ right)}
Deixe :
x
=
π
2
{\ displaystyle x = {\ frac {\ pi} {2}}}
⇒
2
π
=
∏
n
=
1
∞
(
1
-
1
4
n
2
)
⇒
π
2
=
∏
n
=
1
∞
(
4
n
2
4
n
2
-
1
)
=
∏
n
=
1
∞
(
2
n
2
n
-
1
⋅
2
n
2
n
+
1
)
=
2
1
⋅
2
3
⋅
4
3
⋅
4
5
⋅
6
5
⋅
6
7
⋯
{\ displaystyle {\ begin {align} \ Rightarrow {\ frac {2} {\ pi}} & = \ prod _ {n = 1} ^ {\ infty} \ left (1 - {\ frac {1} {4n ^ {2}}} \ right) \\ [6pt] \ Rightarrow {\ frac {\ pi} {2}} & = \ prod _ {n = 1} ^ {\ infty} \ left ({\ frac {4n ^ {2}} {4n ^ {2} -1}} \ right) \\ [6pt] & = \ prod _ {n = 1} ^ {\ infty} \ left ({\ frac {2n} {2n- 1}} \ cdot {\ frac {2n} {2n + 1}} \ right) = {\ frac {2} {1}} \ cdot {\ frac {2} {3}} \ cdot {\ frac {4 } {3}} \ cdot {\ frac {4} {5}} \ cdot {\ frac {6} {5}} \ cdot {\ frac {6} {7}} \ cdots \ end {alinhados}}}
Relação com a aproximação de Stirling
A aproximação de Stirling para a função fatorial afirma que
n
!
{\ displaystyle n!}
n
!
=
2
π
n
(
n
e
)
n
[
1
+
O
(
1
n
)
]
.
{\ displaystyle n! = {\ sqrt {2 \ pi n}} {\ left ({\ frac {n} {e}} \ right)} ^ {n} \ left [1 + O \ left ({\ frac {1} {n}} \ direita) \ direita].}
Considere agora as aproximações finitas do produto Wallis, obtidas tomando os primeiros termos do produto
k
{\ displaystyle k}
p
k
=
∏
n
=
1
k
2
n
2
n
-
1
2
n
2
n
+
1
,
{\ displaystyle p_ {k} = \ prod _ {n = 1} ^ {k} {\ frac {2n} {2n-1}} {\ frac {2n} {2n + 1}},}
onde pode ser escrito como
p
k
{\ displaystyle p_ {k}}
p
k
=
1
2
k
+
1
∏
n
=
1
k
(
2
n
)
4
[
(
2
n
)
(
2
n
-
1
)
]
2
=
1
2
k
+
1
⋅
2
4
k
(
k
!
)
4
[
(
2
k
)
!
]
2
.
{\ displaystyle {\ begin {alinhados} p_ {k} & = {1 \ over {2k + 1}} \ prod _ {n = 1} ^ {k} {\ frac {(2n) ^ {4}} { [(2n) (2n-1)] ^ {2}}} \\ [6pt] & = {1 \ sobre {2k + 1}} \ cdot {{2 ^ {4k} \, (k!) ^ { 4}} \ over {[(2k)!] ^ {2}}}. \ End {alinhado}}}
Substituindo a aproximação de Stirling nesta expressão (para e ) pode-se deduzir (após um pequeno cálculo) que converge para as .
k
!
{\ displaystyle k!}
(
2
k
)
!
{\ displaystyle (2k)!}
p
k
{\ displaystyle p_ {k}}
π
2
{\ displaystyle {\ frac {\ pi} {2}}}
k
→
∞
{\ displaystyle k \ rightarrow \ infty}
Derivada da função zeta de Riemann em zero
A função zeta de Riemann e a função eta de Dirichlet podem ser definidas:
ζ
(
s
)
=
∑
n
=
1
∞
1
n
s
,
ℜ
(
s
)
>
1
η
(
s
)
=
(
1
-
2
1
-
s
)
ζ
(
s
)
=
∑
n
=
1
∞
(
-
1
)
n
-
1
n
s
,
ℜ
(
s
)
>
0
{\ displaystyle {\ begin {alinhados} \ zeta (s) & = \ sum _ {n = 1} ^ {\ infty} {\ frac {1} {n ^ {s}}}, \ Re (s)> 1 \\ [6pt] \ eta (s) & = (1-2 ^ {1-s}) \ zeta (s) \\ [6pt] & = \ sum _ {n = 1} ^ {\ infty} { \ frac {(-1) ^ {n-1}} {n ^ {s}}}, \ Re (s)> 0 \ end {alinhado}}}
Aplicando uma transformada de Euler à última série, obtém-se o seguinte:
η
(
s
)
=
1
2
+
1
2
∑
n
=
1
∞
(
-
1
)
n
-
1
[
1
n
s
-
1
(
n
+
1
)
s
]
,
ℜ
(
s
)
>
-
1
⇒
η
′
(
s
)
=
(
1
-
2
1
-
s
)
ζ
′
(
s
)
+
2
1
-
s
(
em
2
)
ζ
(
s
)
=
-
1
2
∑
n
=
1
∞
(
-
1
)
n
-
1
[
em
n
n
s
-
em
(
n
+
1
)
(
n
+
1
)
s
]
,
ℜ
(
s
)
>
-
1
{\ displaystyle {\ begin {alinhados} \ eta (s) & = {\ frac {1} {2}} + {\ frac {1} {2}} \ sum _ {n = 1} ^ {\ infty} (-1) ^ {n-1} \ left [{\ frac {1} {n ^ {s}}} - {\ frac {1} {(n + 1) ^ {s}}} \ right], \ Re (s)> - 1 \\ [6pt] \ Rightarrow \ eta '(s) & = (1-2 ^ {1-s}) \ zeta' (s) + 2 ^ {1-s} (\ ln 2) \ zeta (s) \\ [6pt] & = - {\ frac {1} {2}} \ sum _ {n = 1} ^ {\ infty} (- 1) ^ {n-1} \ esquerda [{\ frac {\ ln n} {n ^ {s}}} - {\ frac {\ ln (n + 1)} {(n + 1) ^ {s}}} \ direita], \ Re ( s)> - 1 \ end {alinhado}}}
⇒
η
′
(
0
)
=
-
ζ
′
(
0
)
-
em
2
=
-
1
2
∑
n
=
1
∞
(
-
1
)
n
-
1
[
em
n
-
em
(
n
+
1
)
]
=
-
1
2
∑
n
=
1
∞
(
-
1
)
n
-
1
em
n
n
+
1
=
-
1
2
(
em
1
2
-
em
2
3
+
em
3
4
-
em
4
5
+
em
5
6
-
⋯
)
=
1
2
(
em
2
1
+
em
2
3
+
em
4
3
+
em
4
5
+
em
6
5
+
⋯
)
=
1
2
em
(
2
1
⋅
2
3
⋅
4
3
⋅
4
5
⋅
⋯
)
=
1
2
em
π
2
⇒
ζ
′
(
0
)
=
-
1
2
em
(
2
π
)
{\ displaystyle {\ begin {alinhados} \ Rightarrow \ eta '(0) & = - \ zeta' (0) - \ ln 2 = - {\ frac {1} {2}} \ sum _ {n = 1} ^ {\ infty} (- 1) ^ {n-1} \ left [\ ln n- \ ln (n + 1) \ right] \\ [6pt] & = - {\ frac {1} {2}} \ sum _ {n = 1} ^ {\ infty} (- 1) ^ {n-1} \ ln {\ frac {n} {n + 1}} \\ [6pt] & = - {\ frac {1 } {2}} \ left (\ ln {\ frac {1} {2}} - \ ln {\ frac {2} {3}} + \ ln {\ frac {3} {4}} - \ ln { \ frac {4} {5}} + \ ln {\ frac {5} {6}} - \ cdots \ right) \\ [6pt] & = {\ frac {1} {2}} \ left (\ ln {\ frac {2} {1}} + \ ln {\ frac {2} {3}} + \ ln {\ frac {4} {3}} + \ ln {\ frac {4} {5}} + \ ln {\ frac {6} {5}} + \ cdots \ right) \\ [6pt] & = {\ frac {1} {2}} \ ln \ left ({\ frac {2} {1}} \ cdot {\ frac {2} {3}} \ cdot {\ frac {4} {3}} \ cdot {\ frac {4} {5}} \ cdot \ cdots \ right) = {\ frac {1} {2}} \ ln {\ frac {\ pi} {2}} \\\ Rightarrow \ zeta '(0) & = - {\ frac {1} {2}} \ ln \ left (2 \ pi \ right ) \ end {alinhado}}}
Veja também
Notas
links externos
<img src="https://en.wikipedia.org/wiki/Special:CentralAutoLogin/start?type=1x1" alt="" title="" width="1" height="1" style="border: none; position: absolute;">