What is the radius of the circle $|z - z_0| = \rho |z - z_1|$? My conjecture is $$\text{radius} = |z_1 - z_0| \cdot \frac {\rho ^2}{(\rho -1)^2(\rho + 1)^2}$$ which I request help verifying, proving, or correcting.
Note: Formulas exist; I'm asking for help with this approach, as described below.
Let $d = |z_1 - z_0|$. Using the approach taught in Street Fighting Mathematics, I deduced the following:
- The radius is obviously a function of $d$ and $\rho$.
- It must be undefined when $\rho = 1$.
- This also implies that $\rho$ is unitless.
- It must be linear in $d$, since the radius is in units of linear distance, and $d$ is as well (whereas $\rho$ is unit less).
- It must be symmetric when switching $z_0, z_1$, implying $f(\rho) = f(\frac 1 \rho)$.
- It must approach $0$ when $\rho \to 0$ or $\rho \to \infty$.
- It must approach $\infty$ when $\rho$ approaches $1$ (illustrated by Circle of Appolonius diagrams).
Asking What function $f(x) = f(1/x)$ has limits $\lim_{x \to 0} f(x) = 0, \lim_{x \to \infty} f(x) = 0$, and $\lim_{x \to 1} f(x) = \infty$? , I got some candidate functions, took the simplest, realized I needed to square $\rho$ to makes things work (which is intuitive given that distance requires taking a square root, and so e.g. the standard equation of circles is in $r^2$), and got this candidate: $$\text{radius} = \frac d {\rho^2 + \frac 1 {\rho^2} - 2}\\ = d \cdot \frac {\rho ^2}{(\rho -1)^2(\rho + 1)^2}.$$
Can you verify or correct this? Is there a way I could develop this approach into a proof?

You have made an amazing list of 7 observations on the radius of the circle $|z - z_0| = \rho |z - z_1|$.
Unfortunately, they are not enough to single out the right formula for the radius.
Here is another observation.
Fix $d=|z_0- z_1|$. When $\rho$ is very small, $z$ is very near to $z_0$. Hence, $|z-z_1|$ is roughly the same as $|z_0-z_1|=d$. So, $|z-z_0|=\rho|z-z_1|\approx\rho d$. That is, the radius $r\approx\rho d$. More precisely, when $\rho\to0$, $\frac{r}{\rho d}\to1$.
Since $r$ is linear in $d$, let $r(\rho,d)=g(\rho)d$ for some function $g(\rho)$, $\rho>0$. Then our observations mean
The simplest function that satisfies the observations above is $$g(\rho)=\frac1{|\rho-\frac1\rho|}.$$
Hence, we could hope the answer might be
$$r(\rho,d)=\frac d{|\rho-\frac1\rho|}.$$
Fortunately and magically, it is correct!