aheight said:
Depends on what you mean by "solve". Just because an integral cannot be expressed in terms of simple functions like polynomials or radicals or trig functions doesn't mean it can't be "solved". There are many special functions that represent such integrals. For example:
##\int \frac{\sin(x)}{x}dx=\text{Si}(x)##
And the sine function Si(x) is just a function like ##\sin(x)##. Why then cannot we just say:
##\int \log(e^{x^2}-1)dx=\text{myfunction}(x)##
Also, how about a (convergent) power series solution? Is that not a legitimate solution just like sine or cosine? For example:
##\log(e^{x^2}+1)=\log(2)+\frac{x^2}{2}+\frac{x^4}{8}+\frac{x^8}{192}+\frac{x^{12}}{2280}+\cdots##
then can't we write:
## \int \log(e^{x^2}+1)dx=\int \big[ \log)(2)+\frac{x^2}{2}+\frac{x^4}{8}+\frac{x^8}{192}+\frac{x^{12}}{2280}+\cdots \big]dx##
Can you integrate that power series solution and come out with a closed-form solution like:
##\int \log(e^{x^2}+1)dx=\log(2)x+\sum_{n=1}^{\infty} (\text{some expression here})##?
Let
$$F(x) = \int_0^x \ln \left( 1 + e^{\phi^2} \right) \, d\phi .$$
The series solution proposed above will converge only for small values of ##|x|##. A much more effective series solution (involving non-elementary but well-studied and readily-available functions) is obtained by writing
$$ \ln\left(1+e^{\phi^2}\right) = \ln\left( e^{\phi^2} \right) + \ln \left( 1 + e^{-\phi^2} \right) = \phi^2 + \ln \left( 1 + e^{-\phi^2} \right).$$
Thus
$$F(x) = \frac{x^3}{3} + \int_0^x \ln \left( 1 + e^{-\phi^2} \right) \, d\phi.$$
Now expand the second logarithm, using ##\ln(1+u) = u - u^2/2 + u^3/3 - \cdots## with ##u = e^{-\phi^2} < 1## (so: convergent). Note that ##u^n = e^{-n \phi^2}##, so we can express the ##n##th term of the series for ##F(x)## as
$$t_n(x) = \frac{(-1)^{n-1}}{n} \int_0^x e^{-n \phi^2} \, d\phi = \frac{(-1)^{n-1} \sqrt{\pi}}{n^{3/2}} \text{erf}( x \sqrt{n}) .$$
Here, "erf" is the non-elementary function
$$\text{erf}(w) = \frac{2}{\sqrt{\pi}} \int_0^w e^{-t^2} \, dt$$.
The infinite series above converges for all ##x > 0.##
In sophisticated computer packages the function "erf" is built-in, so numerical evaluation of the terms of ##F(x)## is not very difficult. (However, to be fair, I must say that for sizable values of ##x>0## we may need to include a large number of terms in order to obtain good accuracy.)