# Solving and Plotting an integral function

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

I want to integrate and plot the following integral function:

$$F(x)=\displaystyle\int_x^\infty\dfrac{t\sqrt{t^2-x^2}}{e^t-1}\text{d}t$$:

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.

• Possible duplicate of What are the most common pitfalls awaiting new users? Apr 15, 2019 at 18:59
• @MarcoB. Please indicate which answer(s) are relevant. Apr 15, 2019 at 21:34
• Please provide the code you used. Apr 28, 2019 at 14:20

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}]


• 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?
– Lele
Apr 15, 2019 at 14:21
• @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. Apr 15, 2019 at 16:16
• 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?
– Lele
Apr 16, 2019 at 10:22
• The sum(!) doesn't depend on x and therefor doesn't define F[x]! Apr 16, 2019 at 10:31
• 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$
– Lele
Apr 16, 2019 at 10:45