Detail from Integral Evaluation in Szabo (Appendix)
As commented, do note that the integral (A.20) is taken over all $\mathbf{k}$-space, whereas the integral (A.21) is one-dimensional. The explicit calculation of this integral follows as 

$$\begin{align} \int\text{d}\mathbf{k}\,\Psi(k)e^{-i\mathbf{k\cdot r}} &amp;=\int_0^\infty\text{d}k\,k^2\int_0^\pi\text{d}\theta\int_0^{2\pi}\text{d}\varphi\,\sin\theta\Psi(k)e^{-ikr\cos\theta}\\
&amp;= 2\pi\int_0^\infty\text{d}k\,k^2\Psi(k)\int_0^\pi\text{d}\theta\,\sin\theta e^{-ikr\cos\theta}\\
&amp;= 2\pi\int_0^\infty\text{d}k\,k^2\Psi(k)\left[\frac{e^{ikru}}{ikr}\right]_{u=-1}^1\\
&amp;= \frac{4\pi}{r}\int_0^\infty\text{d}k\,k\Psi(k)\sin(kr).\end{align}$$

In the first equality we convert this integral from Cartesian to spherical coordinates, in the second equality we evaluate the integral over $\varphi$, in the third equality we substitute $u = -\cos\theta$, and in the fourth equality we make use of Euler's formula to get the exponential representation of $\sin(kr)$.