Ytukyg

I've found a number of closely related functions, each of which takes three complex numbers $z_1,,z_2$, and $z_3$ (which we can consider as the vertices of a triangle) as its arguments, and outputs a complex number $w$. A few of these functions represent well-known triangle centers $w$, such as the centroid (that is, $w=(z_1+z_2+z_3)/3)$), or the zeros and critical points of the polynomial $R(z) := (z-z_1)(z-z_2)(z-z_3)$. Other functions among them are less obvious, such as the following point $w=p$:

$$p := frac{(z_2+z_3) z_1^2+left(z_2^2-6z_2z_3+z_3^2right) z_1+z_2z_3(z_2+z_3)sqrt{-3(z_1-z_2)^2(z_1-z_3)^2(z_2-z_3)^2}}{2 left(z_1^2+z_2^2+z_3^2-z_1z_2-z_1z_3-z_2z_3right)}$$

As far as I can tell, by comparing the point $p$ to different triangle centers in GeoGebra, it seems to be the first isodynamic point to a high degree of precision for the triangles that I've thrown at it so far. My problem is twofold:

$(1)$ How do I prove that $p$ is (or isn't) the first isodynamic point of the triangle ($z_1,,z_2,,z_3$)?

$(2)$ How do I (numerically) match $w$ with likely triangle centers based on the Encyclopedia of Triangle Centers? GeoGebra uses a small portion of the material in the ETC, but unfortunately, it fails to implement some important triangle centers with low indices. Suggestions for other, similar resources are welcome, too.

As for $(1)$, I think a viable starting point might be to convert $p$ into barycentric or normalized trilinear coordinates and directly compare the result to the coordinates given in the ETC, but I'm not sure how to do this. It seems to be slightly easier to use Cartesian coordinates as a starting point, but Mathematica (for example) generally doesn't seem to be very inclined to return the real and imaginary parts of $p$ (for general $z_k$, e.g. when you let $z_k = a_k + i b_k$ for real numbers $a_k,,b_k,,k=1,2,3,$ and simplify the expressions).

Regarding $(2)$, the following website can supposedly be used to compare points to almost all entries in the ETC: https://faculty.evansville.edu/ck6/encyclopedia/Search_6_9_13.html

Unfortunately, I seem to get inconsistent results from my attempts to implement the algorithm on the website above in Mathematica. More specifically, certain triangle centers that I've tested (such as the centroid) yield the correct coordinates in the list, while other obvious ones that I've tested are either not listed, or have the wrong index. I suspect that I've misunderstood the information on the website in some elementary way, so feel free to correct my algorithm (for $p$ above, in the example that follows):

$(i)$ Choose $z_1,,z_2,,z_3$ as the vertices of a triangle with side lengths $6,,9,$ and $13$, e.g. $z_1 = 0,,z_2 = 9,$ and $z_3 = (-26/9) + (8/9)sqrt{35}i$.

$(ii)$ Solve the following linear system of equations for $u,v,w$:

$$begin{cases}
u,text{Re}(z_1) + v,text{Re}(z_2) + w,text{Re}(z_3) = text{Re}(p)\
u,text{Im}(z_1) + v,text{Im}(z_2) + w,text{Im}(z_3) = text{Im}(p)\
u + v + w = 1,
end{cases}$$

where $u:v:w$ are the barycentric coordinates for $p$. (Here, $p$ can be approximated to a high degree of precision, e.g. to avoid problems with Mathematica...)

$(iii)$ Let $a=6,,b=9$, and $c=13$ (or, more ideally, just use the Pythagorean theorem...), and define $x = u/a,,y = v/b,,z = w/c$. Furthermore, calculate the area of the triangle as $A = 4sqrt{35}$.

$(iv)$ Calculate $kx = 2Ax/(ax+by+cz) = 2Ax/(u+v+w)$, which should be the sought "coordinate" in the table.

Carrying out the algorithm above yields $kx approx 0.14368543660$, whereas the coordinate for the first isodynamic point (with index $15$ in the table) is $sim 3.10244402065$. Perhaps $p$ isn't actually the first isodynamic point, but I like to think that I've done something wrong, like interpreting $a,,b,$ and $c$ as a triangle's side lengths...