It seems that sagemath is unable to calculate Fermic-Dirac type integrals, e. g. integrate(x^2/(1+ e^x),x,0,oo) => limit(1/3x^3 - x^2log(e^x + 1) - 2xdilog(-e^x) + 2polylog(3, -e^x), x, +Infinity, minus) + 3/2zeta(3)
integrate(x^3/(1+ e^x),x,0,oo)==> -7/120pi^4 + limit(1/4x^4 - x^3log(e^x + 1) - 3x^2dilog(-e^x) + 6xpolylog(3, -e^x) - 6polylog(4, -e^x), x, +Infinity, minus)
Normally, the frormer valuates to (3Zeta[3])/2 and the latter to (7Pi^4)/120, using Mathematica.