Aloha :)
Die gegebene Dichtefunktion hat Kugelsymmetrie, weil sie nur vom Betrag des Orsvektors abhängt. Daher drängen sich hier Kugelkoordinaten geradezu auf:$$\vec r=\begin{pmatrix}r\cos\varphi\sin\vartheta\\r\sin\varphi\sin\vartheta\\r\cos\vartheta\end{pmatrix}\quad;\quad r\in[0;\infty)\quad;\quad\varphi\in[0;2\pi]\quad;\quad\vartheta\in[0;\pi]$$Die Masse \(M\) der Gaswolke ist damit:
$$M=\int\limits_V\rho(\vec r)\,dV=\int\limits_{r=0}^\infty\;\int\limits_{\varphi=0}^{2\pi}\;\int\limits_{\vartheta=0}^\pi\frac{2}{r^2(r+1)^2}\,r^2\sin\vartheta\,dr\,d\varphi\,d\vartheta$$$$\phantom{M}=\int\limits_{r=0}^\infty\frac{2}{(r+1)^2}\,dr\cdot\int\limits_{\varphi=0}^{2\pi}d\varphi\cdot\int\limits_{\vartheta=0}^\pi\sin\vartheta\,d\vartheta=\left[-\frac{2}{r+1}\right]_0^\infty\cdot\left[\varphi\right]_0^{2\pi}\cdot\left[-\cos\varphi\right]_0^\pi$$$$\phantom{M}=\left(0+2\right)\cdot(2\pi-0)\cdot(1-(-1))=8\pi$$