Cómo resolver [math] \ displaystyle \ int_ {0} ^ {\ infty} \! \ Frac {x ^ 3} {e ^ x-1} \ dx [/ math]

Una definición muy conocida de la función Riemann Zeta …

[matemáticas] \ en caja {\ color {azul} {\ implica \ displaystyle \ int_ {0} ^ {\ infty} \ dfrac {u ^ {x-1}} {e ^ {u} -1} \, du = \ zeta (x). \ Gamma (x)}} [/ math]

aquí [matemáticas] (x-1) = 3 \ implica x = 4 [/ matemáticas]

entonces la integral será

[matemáticas] \ large \ boxed {I = \ zeta (4) \ Gamma (4) = 6. \ dfrac {\ pi ^ {4}} {90} = \ dfrac {\ pi ^ {4}} {15} }[/matemáticas]


Prueba de la propiedad anterior de Riemann Zeta :

Primero, la función gamma se define como ..

[matemáticas] \ displaystyle \ int_ {0} ^ {+ \ infty} z ^ {x-1} e ^ {- z} \, dz = \ Gamma (x) [/ math]

ahora tomemos una serie …

[matemáticas] \ dfrac {1} {e ^ u-1} = e ^ {- u} + e ^ {- 2u} + e ^ {- 3u} + \ ldots [/ matemáticas]

entonces la integral se convierte

[matemáticas] \ displaystyle \ int_ {0} ^ {\ infty} \ dfrac {u ^ {x-1}} {e ^ {u} -1} \, du = \ displaystyle \ sum_ {n \ ge 1} \ displaystyle \ int_ {0} ^ {\ infty} u ^ {x-1} e ^ {- nu} \, du [/ math]

ahora reemplaza [math] u = \ frac {z} {n} [/ math], tendremos

[matemáticas] I = \ displaystyle \ sum_ {n> 1} \ left (\ dfrac {1} {n} \ right) ^ {x} \ displaystyle \ int_ {0} ^ {\ infty} z ^ {x-1 } e ^ {- z} \, dz = \ zeta (x) \ Gamma (x) [/ math]

-> Para otra prueba ver sección de comentarios

[matemáticas] \ displaystyle \ int_0 ^ {\ infty} \ frac {x ^ 3e ^ {- x}} {1 – e ^ {- x}} dx = \ int_0 ^ {\ infty} x ^ 3e ^ {- x } \ sum_ {i = 0} ^ {\ infty} e ^ {- ix} dx = \ sum_ {i = 0} ^ {\ infty} \ int_0 ^ {\ infty} x ^ 3e ^ {- (i + 1 ) x} dx [/ matemáticas]

Sustituya [matemáticas] (i + 1) x = y [/ matemáticas] en cada integral para obtener,

[matemáticas] \ displaystyle \ int_0 ^ {\ infty} \ frac {x ^ 3e ^ {- x}} {1 – e ^ {- x}} dx = \ sum_ {i = 0} ^ {\ infty} \ int_0 ^ {\ infty} \ frac {y ^ 3} {(i + 1) ^ 4} e ^ {- y} dy [/ math]

[matemáticas] \ displaystyle = \ left (\ int_0 ^ {\ infty} y ^ 3e ^ {- y} dy \ right) \ left (\ sum_ {i = 0} ^ {\ infty} \ frac {1} {( i + 1) ^ 4} \ right) = \ Gamma (4) \ zeta (4) [/ math]


En general,

[matemáticas] \ displaystyle \ int_ {0} ^ {\ infty} \ frac {x ^ {s-1} e ^ {- x}} {1 – e ^ {- x}} dx = \ Gamma (s) \ zeta (s) \; \; \; \ forall s> 1 [/ matemáticas]

[Editar] La pregunta ha sido modificada por un bot de Quora. Por lo tanto, lo estoy actualizando en consecuencia. Con suerte, no volverá a cambiar. [\ Editar]

Para comenzar, he respondido una pregunta similar aquí [1].

Marque esta respuesta si desea más detalles (técnicos).

Bien, el truco aquí es mostrar que la cantidad [matemática] I_ {s}: = \ int_ {0} ^ {\ infty} \ dfrac {x ^ {s}} {e ^ {x} -1} dx [ / math] para [math] s \ geq 0 [/ math] puede escribirse como:

[math] I_ {s} = \ Gamma (s + 1) \ zeta (s + 1) [/ math] donde [math] \ zeta [/ math] es la función zeta de Riemann [2] y [math] \ Gamma [/ math] es la función gamma [3].

De manera similar a lo que hice en la respuesta anterior, expanda la fracción en una serie geométrica (después de factorizar [math] e ^ {x} [/ math]), intercambie signos de suma e integración, calcule la integral usando la integración por sustitución y reconociendo la función [math] \ Gamma [/ math].

Eso es. Para obtener la respuesta, inserte [math] s = 3 [/ math] en la identidad y obtendrá:

[matemáticas] I: = I_ {3} = \ Gamma (4) \ zeta (4) = 6 \ dfrac {\ pi ^ {4}} {90} = \ dfrac {\ pi ^ {4}} {15} [/matemáticas]

En términos más generales, se puede demostrar que, para cualquier número entero estrictamente positivo [matemática] n [/ matemática]:

[matemáticas] \ zeta (2n) = (-1) ^ {n + 1} \ dfrac {B_ {2n} (2 \ pi) ^ {2n}} {2 (2n)!} [/ matemáticas] donde [matemáticas ] B_ {2n} [/ math] son ​​los números de Bernoulli [4].

Así [matemáticas] I_ {2n-1} = \ Gamma (2n) \ zeta (2n) = (2n-1)! (- 1) ^ {n + 1} \ dfrac {B_ {2n} (2 \ pi) ^ {2n}} {2 (2n)!} = (- 1) ^ {n + 1} \ dfrac {B_ {2n} (2 \ pi) ^ {2n}} {4n} [/ matemáticas]

[matemática] n = 1 [/ matemática] da [matemática] \ dfrac {B_ {2} (2 \ pi) ^ {2}} {4} = B_ {2} \ pi ^ {2} = \ dfrac {\ pi ^ {2}} {6} [/ math] ya que [math] B_ {2} = \ dfrac {1} {6} [/ math]

[matemáticas] n = 2 [/ matemáticas] da [matemáticas] – \ dfrac {B_ {4} (2 \ pi) ^ {4}} {8} = – 2B_ {4} \ pi ^ {4} = \ dfrac {\ pi ^ {4}} {15} [/ math] ya que [math] B_ {4} = – \ dfrac {1} {30} [/ math]

Y así.

Consulte el siguiente enlace si tiene curiosidad acerca de esta representación especial de la función zeta de Riemann [5].

Notas al pie

[1] ¿Cómo resuelvo: [matemáticas] \ displaystyle \ int_0 ^ {\ infty} \ frac {x} {e ^ {x} -1} dx [/ math]?

[2] Función zeta de Riemann – Wikipedia

[3] Función gamma – Wikipedia

[4] Número de Bernoulli – Wikipedia

[5] Valores particulares de la función zeta de Riemann – Wikipedia

Algunas grandes respuestas ya han aparecido en este hilo. Aquí hay otra forma posible de obtener toda la familia de integrales.

[matemáticas] \ displaystyle I_n: = \ int_0 ^ {\ infty} \! \! \! \ mathrm {d} x ~ \ frac {x ^ {n + 1}} {e ^ x-1}, \ qquad n = 0,1,2… [/ matemáticas]

Comencemos definiendo la serie de potencia formal (llamada función generadora)

[matemáticas] \ displaystyle G (z): = \ sum_ {n = 0} ^ {\ infty} \ frac {I_n} {n!} z ^ n, \ qquad (A) [/ math]

donde los valores permitidos de [math] z [/ math] que garantizan la convergencia de esta serie se especificarán en un segundo. Al conectar la definición de [matemáticas] I_n [/ matemáticas] anterior, obtenemos

[matemáticas] \ displaystyle G (z) = \ sum_ {n = 0} ^ {\ infty} \ frac {z ^ n} {n!} \ int_0 ^ {\ infty} \! \! \! \ mathrm {d } x ~ \ frac {x ^ {n + 1}} {e ^ x-1} [/ math]

[matemáticas] \ displaystyle = \ int_0 ^ {\ infty} \! \! \! \ mathrm {d} x ~ \ frac {x} {e ^ x-1} \ sum_ {n = 0} ^ {\ infty} \ frac {x ^ nz ^ n} {n!} [/ math]

[matemáticas] \ displaystyle = \ int_0 ^ {\ infty} \! \! \! \ mathrm {d} x ~ \ frac {x} {e ^ x-1} \, e ^ {zx}. \ qquad (B )[/matemáticas]

Promover, adicional,

[matemáticas] \ displaystyle G (z) = \ int_0 ^ {\ infty} \! \! \! \ mathrm {d} x ~ \ frac {x} {e ^ x-1} \, e ^ {zx + xx }[/matemáticas]

[matemáticas] \ displaystyle = \ int_0 ^ {\ infty} \! \! \! \ mathrm {d} x ~ \ frac {x} {e ^ x-1} e ^ {(z-1) x} (e ^ x-1 + 1) [/ matemáticas]

[matemáticas] \ displaystyle = \ int_0 ^ {\ infty} \! \! \! \ mathrm {d} x ~ x \, e ^ {(z-1) x} + \ int_0 ^ {\ infty} \! \ ! \! \ mathrm {d} x ~ \ frac {x} {e ^ x-1} \, e ^ {(z-1) x} [/ math]

Tenga en cuenta que la segunda integral es simplemente [matemáticas] G (z-1) [/ matemáticas], mientras que la primera integral converge solo si [matemáticas] z <1 [/ matemáticas].

Dejando [math] y = (1-z) x [/ math] en la primera integral, tenemos

[matemáticas] \ displaystyle \ int_0 ^ {\ infty} \! \! \! \ mathrm {d} x ~ x \, e ^ {(z-1) x} = \ int_0 ^ {\ infty} \ frac {\ mathrm {d} y} {(1-z) ^ 2} \, y \, e ^ {- y} = \ frac {1} {(1-z) ^ 2}, [/ math]

que produce la ecuación funcional

[matemática] \ displaystyle G (z) -G (z-1) = \ frac {1} {(1-z) ^ 2} [/ matemática]

satisfecho por [matemáticas] G (z) [/ matemáticas]. Para resolver esta ecuación, sumamos las ecuaciones

[matemática] \ displaystyle G (z) -G (z-1) = \ frac {1} {(1-z) ^ 2} [/ matemática]

[matemáticas] \ displaystyle G (z-1) -G (z-2) = \ frac {1} {(2-z) ^ 2} [/ matemáticas]

[matemáticas] \ displaystyle G (z-2) -G (z-3) = \ frac {1} {(3-z) ^ 2} [/ matemáticas]

ad infinitum, obteniendo

[matemáticas] \ displaystyle G (z) -G (- \ infty) = \ sum_ {k = 1} ^ {\ infty} \ frac {1} {(kz) ^ 2} [/ math]

donde lhs resulta de la cancelación de términos adyacentes. Además, vemos en la ecuación [matemáticas] (B) [/ matemáticas], que [matemáticas] G (- \ infty) = 0 [/ matemáticas]. Por lo tanto,

[matemáticas] \ displaystyle G (z) = \ sum_ {k = 1} ^ {\ infty} \ frac {1} {(kz) ^ 2} [/ matemáticas]

Ahora, podemos cosechar los frutos de [matemáticas] G (z) [/ matemáticas].

De la ecuación [matemáticas] (A) [/ matemáticas], tenemos

[matemáticas] \ displaystyle I_n = G ^ {(n)} (0), [/ matemáticas]

donde [math] G ^ {(n)} (z) [/ math] es la [math] n ^ {\ text {th}} [/ math] derivada de [math] G [/ math] wrt [math] z [/ matemáticas]. Usando la identidad

[matemáticas] \ displaystyle \ frac {\ mathrm {d} ^ n} {\ mathrm {d} z ^ n} \ frac {1} {(kz) ^ 2} = \ frac {(n + 1)!} { (kz) ^ {n + 2}} [/ matemáticas]

encontramos

[matemáticas] \ displaystyle G ^ {(n)} (z) = \ sum_ {k = 1} ^ {\ infty} \ frac {(n + 1)!} {(kz) ^ {n + 2}}. [/matemáticas]

Finalmente tenemos

[matemáticas] \ displaystyle I_n = \ sum_ {k = 1} ^ {\ infty} \ frac {(n + 1)!} {k ^ {n + 2}} = (n + 1)! \ zeta (n + 2) [/ matemáticas]

Este resultado también sigue la conocida representación integral de la función zeta de Riemann utilizada por Siddhartha Ganguly anteriormente.

En particular, tenemos

[matemáticas] \ displaystyle I_0 = \ int_0 ^ {\ infty} \! \! \! \ mathrm {d} x ~ \ frac {x} {e ^ x-1} = \ zeta (2) = \ frac {\ pi ^ 2} {6} [/ matemáticas]

[matemáticas] \ displaystyle I_2 = \ int_0 ^ {\ infty} \! \! \! \ mathrm {d} x ~ \ frac {x ^ 3} {e ^ x-1} = 3! \ zeta (4) = \ frac {\ pi ^ 4} {15} [/ matemáticas]

[matemáticas] \ displaystyle I_4 = \ int_0 ^ {\ infty} \! \! \! \ mathrm {d} x ~ \ frac {x ^ 5} {e ^ x-1} = 5! \ zeta (6) = \ frac {8 \ pi ^ 6} {63} [/ matemáticas]

PD: Un físico identificaría [matemática] I_n [/ matemática] como el área bajo la curva de distribución de Bose-Einstein para una partícula spin [matemática] n + 1 [/ matemática] con potencial químico [matemática] \ mu [/ matemática] establecido en [matemáticas] 0 [/ matemáticas].

Salud !

A la luz de otras respuestas proporcionadas aquí, señalaría que la integral dada es un caso especial para la transformación Mellin de la función [matemáticas] f (x) = \ dfrac {1} {e ^ x – 1}. [/ matemáticas]

Definimos la transformación de Mellin de una función como:

[matemáticas] \ displaystyle \ large \ boxed {\ mathcal {M} \ left \ {f (x) \ right \} = \ int_ {0} ^ {\ infty} x ^ {p-1} f (x) \ , dx} \\ [/ math]

[matemáticas] \ displaystyle \ mathcal {M} \ left \ {\ dfrac {1} {e ^ x-1} \ right \} = \ int_ {0} ^ {\ infty} x ^ {p-1} \ dfrac {1} {e ^ x-1} \, dx [/ math]

Para evaluar esto, utilizamos el siguiente resultado Infinite GP:

[matemáticas] \ displaystyle \ sum_ {n = 0} ^ {\ infty} e ^ {- nx} = \ dfrac {1} {1-e ^ {- x}} = \ dfrac {e ^ x} {e ^ x – 1} [/ matemáticas]

Dividiendo ambos lados por [matemáticas] e ^ x [/ matemáticas] observamos,

[matemáticas] \ displaystyle \ dfrac {1} {e ^ x – 1} = \ sum_ {n = 0} ^ {\ infty} \ dfrac {e ^ {- nx}} {e ^ x} = \ sum_ {n = 0} ^ {\ infty} e ^ {- (n + 1) x} = \ sum_ {n = 1} ^ {\ infty} e ^ {- nx} [/ math]

Por lo tanto, nuestra transformación Mellin se simplifica como,

[matemáticas] \ displaystyle \ mathcal {M} \ left \ {\ dfrac {1} {e ^ x-1} \ right \} = \ sum_ {n = 1} ^ {\ infty} \ int_ {0} ^ { \ infty} x ^ {p-1} e ^ {- nx} \, dx = \ sum_ {n = 1} ^ {\ infty} \ dfrac {\ Gamma (p)} {n ^ p} = \ Gamma ( p) \ zeta (p) \ tag * {(1)} [/ math]

En la pregunta dada, [matemáticas] p = 4 [/ matemáticas], y por lo tanto la respuesta es

[matemáticas] \ displaystyle \ Gamma (4) \ zeta (4) = \ dfrac {\ pi ^ 4} {15} [/ matemáticas]


Explicando un poco el paso [matemáticas] (1) [/ matemáticas]:

Sabemos que [math] \ displaystyle \ Gamma (p) = \ int_ {0} ^ {\ infty} e ^ {- x} x ^ {p-1} \, dx. [/ Math] Entonces,

[matemáticas] \ displaystyle \ Gamma (p) = \ int_ {0} ^ {\ infty} e ^ {- nx} (nx) ^ {p-1} \, d (nx) [/ math]

[matemáticas] \ displaystyle \ implica \ Gamma (p) = n ^ p \ int_ {0} ^ {\ infty} e ^ {- nx} x ^ {p-1} \, dx \ implica \ int_ {0} ^ {\ infty} e ^ {- nx} x ^ {p-1} \, dx = \ dfrac {\ Gamma (p)} {n ^ p} [/ math]

Y por lo tanto,

[matemáticas] \ displaystyle \ sum_ {n = 1} ^ {\ infty} \ int_ {0} ^ {\ infty} x ^ {p-1} e ^ {- nx} \, dx = \ sum_ {n = 1 } ^ {\ infty} \ dfrac {\ Gamma (p)} {n ^ p} = \ Gamma (p) \ zeta (p) [/ math]

donde [matemáticas] \ displaystyle \ zeta (p) = \ sum_ {n = 1} ^ {\ infty} \ dfrac {1} {n ^ p} \ left (\ Re (p)> 1 \ right) [/ math ] es la función Riemann Zeta.


Solo quiero agregar que si la pregunta se modificó a [matemáticas] \ displaystyle \ int_ {0} ^ {\ infty} \ dfrac {x ^ 3} {e ^ x + 1} \, dx [/ matemáticas] O [matemáticas ] \ displaystyle \ int_ {0} ^ {\ infty} \ dfrac {x ^ 3} {(e ^ x-1) ^ 2} \, dx [/ math] los siguientes resultados pueden ser útiles:

  • [matemáticas] \ matemáticas {M} \ left \ {\ dfrac {1} {e ^ x + 1} \ right \} = \ left (1-2 ^ {1-p} \ right) \ Gamma (p) \ zeta (p) \\ [/ matemáticas]
  • [matemáticas] \ matemáticas {M} \ left \ {\ dfrac {1} {\ left (e ^ x-1 \ right) ^ 2} \ right \} = \ Gamma (p) \ left [\ zeta (p- 1) – \ zeta (p) \ right] [/ math]

La generalización se deja como un ejercicio para el lector 🙂


Gracias a mi hermano Sid por presentarme la Transformación Mellin. Como aspirante a ingeniero eléctrico, diría que es una gran herramienta.

Expande usando la serie geométrica: [matemáticas] \ frac {1} {e ^ x-1} = \ sum_ {n = 0} ^ \ infty e ^ {- (n + 1) x} [/ matemáticas]. Intercambio de suma e integración, y su integral se puede reescribir como: [matemáticas] \ sum_ {n = 0} ^ \ infty \ frac {1} {(n + 1) ^ 4} \ cdot \ int_0 ^ \ infty u ^ 3 e ^ {- u} du [/ math], donde hice la sustitución [math] u = (n + 1) x [/ math].


Esta es una identidad bien conocida que involucra la función gamma y la función zeta. Uno puede probarlo escribiendo el denominador en la integral como una suma infinita (después de las manipulaciones apropiadas).

[matemáticas] \ int_ {0} ^ {\ infty} \ frac {x ^ 3} {e ^ x-1} dx [/ matemáticas]

Una definición muy conocida de la función Riemann Zeta, es decir

[matemáticas] \ zeta (x) = \ frac {1} {\ Gamma (x)} \ int {0} ^ {\ infty} \ frac {u ^ {x-1}} {e ^ u-1} du [/matemáticas]

[matemáticas] \ zeta (x) \ Gamma (x) = \ int {0} ^ {\ infty} \ frac {u ^ {x-1}} {e ^ u-1} du [/ math]

entonces tenemos,

[matemáticas] x-1 = 3 [/ matemáticas]

[matemáticas] \ implica x = 4 [/ matemáticas]

Así, por la función Riemann Zeta tenemos

[matemáticas] \ int_ {0} ^ {\ infty} \ frac {x ^ {4-1}} {e ^ x-1} dx = \ zeta (4) \ Gamma (4) [/ matemáticas]

usa series geométricas infinitas.

Espero que esto responda a su pregunta, en caso de que desee la prueba de la suma de la serie armónica. Házmelo saber.