I am new in the use of Mathematica, so please forgive me..

I want to integrate and plot the following integral function:


where $x>0$ is the argument of my integral function.

Mathematica doesn't manage to do the calculation, not even if I substitute the x inside the square root with a number, e.g. 1.

Can anyone please help me? Thanks


1 Answer 1


Your problem can be solved numerically:

F[x_?NumericQ] :=NIntegrate[(t Sqrt[t^2 - x^2])/(Exp[t] - 1), {t, x, Infinity}]
Plot[F[x], {x, 0, 10}]

enter image description here

  • $\begingroup$ Thank you very much, it was silly I guess... Just a brief question: what is the meaning of using "?NumericQ" in the argument of the function? $\endgroup$
    – Lele
    Apr 15, 2019 at 14:21
  • $\begingroup$ @EmanueleCopello - Look at the documentation for PatternTest. Any function which uses numeric techniques such as NIntegrate can only evaluate if its argument(s) is/are numeric, so the argument(s) should be restricted to being numeric. $\endgroup$
    – Bob Hanlon
    Apr 15, 2019 at 16:16
  • $\begingroup$ Ok thank you. What if I have to evaluate it in this way: $F(x)=\sum_{j=1}^n\displaystyle\int_{x_j}^\infty\dfrac{t\sqrt{t^2-x_j^2}}{e^t-1}\text{d}t$ where $x_j$ should be a list of values I have to insert somehow... What is the best computational way to do it? $\endgroup$
    – Lele
    Apr 16, 2019 at 10:22
  • $\begingroup$ The sum(!) doesn't depend on x and therefor doesn't define F[x]! $\endgroup$ Apr 16, 2019 at 10:31
  • $\begingroup$ Ok sorry, I have to specify that better.. My function is in fact this one: $F(T)=\displaystyle\sum_{j=1}^n\displaystyle\int_{m_j/T}^\infty\dfrac{t\sqrt{t^2-(m_j/T)^2}}{e^t-1}\text{d}t$ $\endgroup$
    – Lele
    Apr 16, 2019 at 10:45

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

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