The problem becomes more managable if you transform to a coordinate system where $\mathbf B$ is at rest. Then $\mathbf A$ has velocity
$$
\pmatrix{\dot x\\\dot y}=-\pmatrix{\frac xr+1\\\frac yr}\;.
$$
Now transform to polar coordinates $r=\sqrt{x^2+y^2}$ and $\phi=\arctan\left(\frac yx\right)$, with
$$
\dot r=\frac{x\dot x}r+\frac{y\dot y}r=-1-\frac xr=-1-\cos\phi
$$
and
$$
\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\;.
$$
So we have
$$
\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}\;.
$$
Dividing by $r$ and integrating both sides yields
$$
\log r=C-2\log\sin\frac\phi2\;.
$$
The initial value is $r\left(\frac\pi2\right)=1$, which yields $C=-\log 2$, so
$$
r=\frac1{2\sin^2\frac\phi2}\;.
$$
In the limit $t\to\infty$ we have $\phi\to\pi$, and thus $r\to\frac12$.