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ática] G (z-1) [/ matemática], mientras que la primera integral converge solo si [matemática] z <1 [/ matemática].
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 !