El problema se vuelve más manejable si usted transformar a un sistema de coordenadas donde $\mathbf B$ está en reposo. A continuación, $\mathbf A$ tiene velocidad
$$
\pmatrix{\dot x\\\dot y}=-\pmatrix{\frac xr+1\\\frac yr}\;.
$$
Ahora transformación a coordenadas polares $r=\sqrt{x^2+y^2}$$\phi=\arctan\left(\frac yx\right)$, con
$$
\dot r=\frac{x\dot x}r+\frac{y\dot y}r=1-\frac xr=-1-\cos\phi
$$
y
$$
\dot\phi=\frac{\frac{\dot y}x-\frac{\dot xy}{y^2}}{1+\left(\frac yx\right)^2}=\frac{x\dot y-y\dot x}{x^2+y^2}=\frac y{r^2}=\frac{\sin\phi}r\;.
$$
Así tenemos
$$
\frac{\mathrm dr}{\mathrm d\phi}=\frac{\dot r}{\dot\phi}=-r\cdot\frac{1+\cos\phi}{\sin\phi}=-r\frac{\cos\frac\phi2}{\sin\frac\phi2}\;.
$$
Dividiendo por $r$ y la integración de ambos lados de los rendimientos
$$
\log r=C-2\log\sin\frac\phi2\;.
$$
El valor inicial es $r\left(\frac\pi2\right)=1$, que los rendimientos de $C=-\log 2$, por lo que
$$
r=\frac1{2\sin^2\frac\phi2}\;.
$$
En el límite de $t\to\infty$ tenemos $\phi\to\pi$, y por lo tanto $r\to\frac12$.