Find the following sum:
$$\sum_{k=1}^\infty \frac{1}{k^2}\frac{x^k}{k!}$$
where $x$ is a real number. This is a power series in $x$. In particular, I'm interested in the case $x>0$.
Disclaimer: This is not a homework exercise, I do not know if a closed form solution exists. If it doesn't, exist, then an approximation in terms of well-known functions (not the all-mighty general hypergeometric $_pF_q$, something simpler please) would be desired.
We know that
$$e^x=1+\sum_{k=1}^\infty\frac{x^k}{k!}$$
Divide both sides to get
$$\frac{e^x}x=\frac1x+\sum_{k=1}^\infty\frac{x^{k-1}}{k!}$$
Integrate both sides to get
$$\int_1^x\frac{e^y}ydy+c_1=\ln(x)+\sum_{k=1}^\infty\frac1k\frac{x^k}{k!}$$
Divide both sides again.
$$\frac1x\left(\int_1^x\frac{e^y}ydy+c_1-\ln(x)\right)=\sum_{k=1}^\infty\frac1k\frac{x^{k-1}}{k!}$$
Integrate both sides to get
which is quite messy?