answered 4 years ago
Looks like a bug indeed. You may use instead
sage: u=integrate(x**3/(exp(x)-1),(x,0,oo),algorithm='giac');u 1/15*pi^4