Pelo jeitão deve ter a ver com Zeta[4], concorda?
(x^3)/(exp(x) - 1) = (x^3)*exp(-x) /(1-exp(-x) ) =>
(x^3)*exp(-x) /(1-exp(-x) ) = SOMA{k=1..oo}(x^3*exp(-k*x)) =>
Agora integrando a série:
int{0..oo}( SOMA{k=1..oo}(x^3*exp(-k*x))) =>
SOMA{k=1..oo}( int{0..oo}(x^3*exp(-k*x))
Bem, agora tabela de integrais ou maple... daí temos:
int(x^3*exp(-k*x)) = -(6+6*k*x+3*x^2*k^2+x^3*k^3)*exp(-k*x)/k^4 =>
int{0..oo}(x^3*exp(-k*x)) = 6/k^4
Voltando à série:
SOMA{k=1..oo}( int{0..oo}(x^3*exp(-k*x))=SOMA{k=1..oo}(6/k^4)
Somas de inversos de potências pares de inteiros é algo bem conhecido, então:
SOMA{k=1..oo}(6/k^4) = 6*Zeta(4)=6*Pi^4/90 = Pi^4/15
Porém esta série não é muito prática para calcular a integral entre outros limites.
[]’s Demetrio
//-----------------------------------------------------------------------