Hi,
When I try to have Maxima do the Planck integral in dimensionless form, integrate(x^3/(exp(x)-1),x,0,inf), I get an answer that involves limits of polylog functions plus the correct answer(pi^4/15) prefaced with a negative sign. Maxima couldn't evaluate the limits of the polylog functions, and so I was disappointed.
An old version of Maple that I have computes the correct symbolic value of this integral quickly and cleanly. I'm not a sophisticated user of Maxima and so I thought that there might be something wrong that I was doing here. Any thoughts would be greatly appreciated and very helpful. Thanks - Jim