Mathematics Stack Exchange is a question and answer site for people studying math at any level and professionals in related fields. Join them; it only takes a minute:

Sign up
Here's how it works:
  1. Anybody can ask a question
  2. Anybody can answer
  3. The best answers are voted up and rise to the top

Is there an asymptotic expansion for the function: \begin{equation} g(x)=\frac{1}{2\zeta(3)}\int_x^\infty \frac{u^2}{e^u - 1} du, \end{equation} over the domain $x\in [0,\infty)$ in terms of elementary functions? Here $\zeta$ is the Riemann zeta function and $1/2\zeta(3)$ is a normalization factor included to ensure $g(0)=1$.

share|cite|improve this question
This is sometimes referred to as the incomplete Bose-Einstein integral. – J. M. Apr 4 '13 at 15:13
up vote 2 down vote accepted

We have $$g(x) = \dfrac1{2 \zeta(3)}\underbrace{\int_x^{\infty} \dfrac{u^2}{e^u-1} du}_{I(x)}$$ We will now obtain a series expansion for $I(x)$. We have $$I(x) = \int_x^{\infty} \dfrac{u^2}{e^u-1} du = \int_x^{\infty} \dfrac{u^2 e^{-u}}{1-e^{-u}} du = \int_x^{\infty} \sum_{k=0}^{\infty} u^2 e^{-(k+1)u} du = \sum_{k=1}^{\infty} \int_x^{\infty} u^2 e^{-ku} du$$ We now have that $$\int_x^{\infty} u^2 e^{-ku} du = \dfrac{e^{-kx} \left(k^2 x^2 + 2kx + 2\right)}{k^3}$$ We have have $$I(x) = \sum_{k=1}^{\infty} \dfrac{e^{-kx} \left(k^2 x^2 + 2kx + 2\right)}{k^3}\tag{$\star$}$$ For a given $x$, truncating $(\star)$ will give you an exponentially converging approximation. Note that $$I(0) = \displaystyle \sum_{k=1}^{\infty} \dfrac2{k^3} = 2 \zeta(3)$$

share|cite|improve this answer
Great! But surely you should substitute $k+1\rightarrow k$ to simplify $(\star)$? Also, the sum is missing in $l(0)$. – Douglas B. Staple Apr 4 '13 at 15:58
@DouglasB.Staple Oh yes:)!. Thanks. Done. – user17762 Apr 4 '13 at 16:00
Beautiful. Also, $(\star)$ gives $ l(x)=2\operatorname{polylog}(3,e^{-x}) + 2x\operatorname{polylog}(2,e^{-x}) - x^2\ln(1-e^{-x})$ using the definition of the $\operatorname{polylog}$ function. – Douglas B. Staple Apr 4 '13 at 16:15
@DouglasB.Staple Yes. I always forget the polyLog expansions, though I am encountering it frequently in recent times. – user17762 Apr 4 '13 at 16:18

yes is equal to the incomplete gamma function

$$ \Gamma (3,x) $$

share|cite|improve this answer
No, it's not.$\phantom{}$ – J. M. Apr 4 '13 at 15:39

Your Answer


By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.