I have three vertices of a triangle $ A = (x_a, y_a) $, $ B = (x_b, y_b) $, and $ C = (x_c, y_c) $, and I need to find the coordinates of a point $ P $ such that the incircles of the triangles $ ACP $, $ BCP $, and $ ABP $ are congruent.

In barycentric coordinates, any point $ P $ can be expressed as:
$ P = \alpha A + \beta B + \gamma C $
where $ \alpha $, $ \beta $, and $ \gamma $ are proportional to the areas of triangles $ PBC $, $ PAC $, and $ PAB $, respectively, and satisfy $ \alpha + \beta + \gamma = 1 $.
Now, we know that the areas of these triangles should be proportional to the semiperimeter times the inradius. If the incircles of the triangles are congruent, this gives us the relations:
$ \text{Area of } PBC \propto (a + PB + PC), \quad \text{Area of } PAC \propto (b + PA + PC), \quad \text{Area of } PAB \propto (c + PA + PB) $
Here, $ a =BC$, $ b =AC$, and $ c =AB$ represent the side lengths of the triangle $ ABC $.
Then I define the following function for any point $ X $:
$ f(X) = w_a A + w_b B + w_c C $
where:
$ w_a = \frac{a + XB + XC}{a + b + c + 2(XA + XB + XC)}, \quad w_b = \frac{b + XA + XC}{a + b + c + 2(XA + XB + XC)}, \quad w_c = \frac{c + XA + XB}{a + b + c + 2(XA + XB + XC)} $
$P$ is the solution to $f(P)=P$ (it is a fixed point of the function).
However, I am stuck at this point due to the nonlinearity of the problem. The direct use of lengths complicates the calculations, leading to square roots and multiple variables in the equations. I would prefer a solution that involves other geometric properties, such as relevant angles or the inradius, which might lead to more elegant equations—even if these equations are only solvable numerically, perhaps involving trigonometric or higher-degree equations with just one variable. Could I approach this problem using such parameters to solve for $ \alpha $, $ \beta $, and $ \gamma $ instead of relying on the lengths?