| 1 | initial version |
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
Copyright Sage, 2010. Some rights reserved under creative commons license. Content on this site is licensed under a Creative Commons Attribution Share Alike 3.0 license.