Aloha :)
Zur Berechnung der Divergenz des Vektorfeldes \(\vec B(\vec r)\) bringen wir den konstanen Vorfaktor auf die linke Seite und definieren uns ein neues Vektorfeld:$$\vec f(\vec r)\coloneqq\frac{4\pi}{\mu_0}\,\vec B(\vec r)=\frac{3(\vec m\,\vec r)\,\vec r}{r^5}-\frac{\vec m}{r^3}$$
Die Komponenten von \(\vec f(\vec r)\) schreiben wir einzeln auf:$$f_i(x_1;x_2;x_3)=\frac{\overbrace{3(m_1x_1+m_2x_2+m_3x_3)\cdot x_i}^{=u}}{\underbrace{\left(x_1^2+x_2^2+x_3^2\right)^{5/2}}_{=v}}-\frac{m_i}{\left(x_1^2+x_2^2+x_3^2\right)^{3/2}}$$und bestimmen die partielle Ableitung von \(f_i\) nach \(x_i\):$$\frac{\partial f_i}{\partial x_i}=\frac{\overbrace{\left(3m_i\cdot x_i+3(m_1x_1+m_2x_2+m_3x_3)\cdot1\right)}^{=u'}\cdot\overbrace{\left(x_1^2+x_2^2+x_3^2\right)^{5/2}}^{=v}}{\underbrace{\left(x_1^2+x_2^2+x_3^2\right)^5}_{=v^2}}$$$$\phantom{\frac{\partial f_i}{\partial x_i}}-\frac{\overbrace{3(m_1x_1+m_2x_2+m_3x_3)\cdot x_i}^{=u}\cdot\overbrace{\frac52\left(x_1^2+x_2^2+x_3^2\right)^{3/2}\cdot2x_i}^{=v'}}{\underbrace{\left(x_1^2+x_2^2+x_3^2\right)^5}_{=v^2}}$$$$\phantom{\frac{\partial f_i}{\partial x_i}}+\frac{m_i\cdot\frac32\left(x_1^2+x_2^2+x_3^2\right)^{1/2}\cdot2x_i}{\left(x_1^2+x_2^2+x_3^2\right)^3}$$$$\phantom{\frac{\partial f_i}{\partial x_i}}=\frac{\left(3m_ix_i+3\vec m\,\vec r\right)\cdot r^5}{r^{10}}-\frac{15\vec m\,\vec r\,r^3\,\cdot x_i^2}{r^{10}}+\frac{3m_ix_i}{r^5}$$$$\phantom{\frac{\partial f_i}{\partial x_i}}=\frac{6m_ix_i+3\vec m\,\vec r}{r^5}-\frac{15\vec m\,\vec r\cdot x_i^2}{r^7}$$
Die Divergenz ist die Summe der drei partiellen Ableitungen \(i=1,2,3\):$$\operatorname{div}(\vec f)=\frac{6m_\pink1x_\pink1+3\vec m\,\vec r}{r^5}+\frac{6m_\green2x_\green2+3\vec m\,\vec r}{r^5}+\frac{6m_\blue3x_\blue3+3\vec m\,\vec r}{r^5}$$$$\phantom{\operatorname{div}(\vec f)}-\frac{15\vec m\,\vec r\cdot x_{\pink1}^2}{r^7}-\frac{15\vec m\,\vec r\cdot x_{\green2}^2}{r^7}-\frac{15\vec m\,\vec r\cdot x_{\blue3}^2}{r^7}$$$$\phantom{\operatorname{div}(\vec f)}=\frac{6(m_1x_1+m_2x_2+m_3x_3)+9\vec m\,\vec r}{r^5}-\frac{15\vec m\,\vec r\,(x_1^2+x_2^2+x_3^2)}{r^7}$$$$\phantom{\operatorname{div}(\vec f)}=\frac{6\vec m\,\vec r+9\vec m\,\vec r}{r^5}-\frac{15\vec m\,\vec r\,r^2}{r^7}=\frac{15\vec m\,\vec r}{r^5}-\frac{15\vec m\,\vec r}{r^5}=0$$
Damit verschwindet auch die Divergenz des Vektorfeldes \(\vec B(\vec r)\).