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
 
//-----------------------------------------------------------------------