I will first numerically/graphically verify the statement to be proved. Given $\triangle ABC$, it can be transformed into an isosceles right triangle, with vertices $L(0,0), M(1,0), N(0, 1)$ using the affine transformation $\mathbf{y} = \mathbf{A x + b}$, for some matrix $\mathbf{A}$ and some vector $\mathbf{b}$.
Clearly, an inscribed ellipse of the original triangle is transformed into an inscribed ellipse of the transformed triangle, and the line segments connecting the vertices to the touching points of the inscribed ellipse in the original triangle will be transformed into line segments in the transformed triangle. So if this set of the three line segments is concurrent, then so will be the set of the corresponding three line segments in the transformed triangle and ellipse, and vice versa.
From this, it is sufficient to show that the line segments connecting the vertices of the transformed triangle to the touching points of the transformed ellipse are concurrent, to show that they are concurrent in the original triangle.
Now suppose we are given the relative location of the touching points of the ellipse on sides $AB$ and $AC$, then these are transformed into known points $T_1$ and $T_2$ on $LM$ and $LN$ respectively, because the relative location of points on line segments does not change after an affine transformation, as defined above. The location of the touching point $T_3$ on side $MN$ is not arbitrary and can be found analytically.
Now we have three touching points, and three vertices, so we just connect them in pairs and see if they intersect at the same point. This is what I have done, and the figure below illustrates this. Thus I have verified that the line segments connecting vertices to the touching points of an in-ellipse on the opposite sides do indeed intersect in a single point.

To turn the above verification into a proof, let the touching point of the inscribed ellipse with the horizontal side be $T_1 = (t_1, 0)$ and the touching point with the vertical side be $T_2 = (0, t_2)$, and let the equation of the inscribed ellipse be
$ (r - r_0)^T Q (r - r_0) = 1$
where $r = [x,y]^T $,
Now the gradient of the ellipse at $T_1$ is given by
$ g_1 = 2 Q (T_1 - r_0) $
And we know that this vector will be pointing along $(-e_2) = - (0, 1)$, therefore,
$ Q (T_1 - r_0) = - \alpha e_2 , \hspace{10pt} \alpha \gt 0 $
From this
$ T_1 - r_0 = - \alpha Q^{-1} e_2 $
Since $T_1$ is on the ellipse, it follows that
$\alpha = \dfrac{1}{\sqrt{ e_2^T Q^{-1} e_2} } $
So that
$ T_1- r_0 = - \dfrac{ Q^{-1} e_2 }{ \sqrt{ e_2^T Q^{-1} e_2 } } $
Premultiplying this equation with $e_1^T $ and $e_2^T$ respectively, results in
$ e_1^T (T_1 - r_0) = t_1 - r_{0x} = - \dfrac{ e_1^T Q^{-1} e_2}{ \sqrt{ e_2^T Q^{-1} e_2 } } $
and
$ e_2^T (T_1 - r_0 ) = 0 - r_{0y} = - \sqrt{e_2^T Q^{-1} e_2 } $
Similarly, the gradient at $T_2$ is pointing in the negative $e_1$ direction, and applying the same method as above, we get
$ - r_{0x} = - \sqrt{ e_1^T Q^{-1} e_1 } $
and
$ t_2 - r_{0y} = - \dfrac{ e_2^T Q^{-1} e_1 }{ e_1^T Q^{-1} e_1 } $
Now let matrix $P = Q^{-1}$ , then
$ e_1^T Q^{-1} e_1 = P_{11} , e_1^T Q^{-1} e_2 = e_2^T Q^{-1} e_1 = P_{12} = P_{21} , e_2^T Q^{-1} e_2 = P_{22} $
Therefore, the equations we have now are
$t_1 - r_{0x} = - \dfrac{ P_{12} }{\sqrt{ P_{22} } } $
$ - r_{0y} = - \sqrt{ P_{22}} $
$ - r_{0x} = - \sqrt{P_{11}} $
$ t_2 - r_{0y} = - \dfrac{ P_{12} } { \sqrt{ P_{11} }} $
Eliminating $r_{0x}$ using the first and third equations, and eliminating $r_{0y}$ using the second and fourth equations, gives us,
$ t_1 = - \dfrac{ P_{12} }{\sqrt{ P_{22} } } + \sqrt{P_{11}} $
$ t_2 = - \dfrac{ P_{12} } { \sqrt{ P_{11} } } + \sqrt{P_{22}} $
And these two equations yield
$ t_1 \sqrt{P_{22} } = - P_{12} + \sqrt{P_{11}} \sqrt{P_{22}} $
$ t_2 \sqrt{P_{11}} = - P_{12} + \sqrt{P_{11}} \sqrt{P_{22}} $
Subtracting,
$ t_1 \sqrt{P_{22}} = t_2 \sqrt{P_{11} } $
Let $X = \sqrt{P_{11}} $, then
$P_{11} = X^2 $
$ P_{22} = \left( \dfrac{t_2}{t_1} \right)^2 X^2 $
$ P_{12} = - t_2 X + \left( \dfrac{t_2}{t_1} \right) X^2 $
and we also have
$ r_{0x} = \sqrt{ P_{11} } = X $
$ r_{0y} = \sqrt{ P_{22} } = \left(\dfrac{t_2}{t_1} \right) X $
What is left is to determine the value of $X$. For that, consider the unknown point $T_3$ which lies on side $MN$, the outward normal is
$ u = \dfrac{1}{\sqrt{2}} [ 1, 1]^T $
Applying the same method given above, we deduce that
$ T_3 - r_0 = \dfrac{ P u }{\sqrt{ u^T P u } } \hspace{25pt} (*) $
Premultiplying by $u^T $
$ u^T (T_3 - r_0) = u^T (T_3 - M + M - r_0) = 0 + \dfrac{1}{\sqrt{2}} \left( 1 - r_{0x} - r_{0y} \right) = \sqrt{ u^T P u } $
Using the above equations, this equation becomes
(where $\alpha = \dfrac{t_2}{t_1} $ )
$ \dfrac{1}{\sqrt{2}}( 1 - X - \alpha X ) = \sqrt{ \frac{1}{2}\left( X^2 + \alpha^2 X^2 + 2 ( - t_1 X + \alpha X^2 ) \right) } \hspace{25pt}(**) $
Squaring both sides of this equation gives
$ ( 1 + (1 + \alpha)^2 X^2 - 2 (1 + \alpha) X ) = X^2 (1 + \alpha)^2 - 2 t_1 X $
from this it follows that
$ X = \dfrac{1}{ 2 ( 1 + \alpha - t_1 ) } = \dfrac{ t_1}{ 2 (t_1 + t_2 - t_1 t_2 ) }$
Let $ d= 2 (t_1 + t_2 - t_1 t_2 ) $, then $ X = \dfrac{t_1}{d} $
thus we have shown that the center of the inscribed ellipse is
$ r_0 = [ \dfrac{t_1}{d} , \dfrac{t_2}{d} ] $
And
$ Q = P^{-1} = \bigg(\dfrac{1}{d^2} \begin{bmatrix} t_1^2 && t_1 t_2 (1 - d) \\ t_1 t_2 (1 - d) && t_2^2 \end{bmatrix} \bigg) ^{-1} = \dfrac{d}{ t_1^2 t_2^2 (2 - d) } \begin{bmatrix} t_2^2 && t_1 t_2 (d-1) \\ t_1 t_2 (d -1) && t_1^2 \end{bmatrix}$
To find the third tangency point, we use equation $(*)$, and we have
$ \dfrac{ P u}{\sqrt{u^T P u} } = \dfrac{[ t_1^2 + t_1 t_2 (1 - d) , t_2^2 + t_1 t_2 (1 - d) ]^T}{d^2 (1 - (1 + \alpha) X)}$
With $\alpha = \dfrac{t_2}{t_1} , X = \dfrac{t_1}{d} $ this becomes
$ \dfrac{ P u}{\sqrt{u^T P u} } = \dfrac{[ t_1^2 + t_1 t_2 (1 - d) , t_2^2 + t_1 t_2 (1 - d) ]^T}{ d (d - (t_1 + t_2)) }$
Hence,
$T_3 = r_0 + \dfrac{ P u}{\sqrt{u^T P u} } = [\dfrac{t_1}{d} , \dfrac{t_2}{d} ]^T + \dfrac{[ t_1^2 + t_1 t_2 (1 - d) , t_2^2 + t_1 t_2 (1 - d) ]^T}{ d (d - (t_1 + t_2)) }$
And the right hand side simplifies to
$T_3 = \dfrac{1}{d (d- t_1 - t_2)} [ t_1 (d - t_1 - t_2) + t_1^2 + t_1 t_2 (1 - d) , t_2 (d - t_1 - t_2) + t_2^2 + t_1 t_2 (1- d) ]^T $
And this simplifies further to the final form
$ T_3 = \dfrac{1}{d - t_1 - t_2} ( t_1 (1 - t_2), t_2 (1 - t_1) ) $
Now the intersection of $NT_1$ and $MT_2$ is easily shown to be
$ G = \dfrac{1}{1 - t_1 t_2} (t_1 ( 1- t_2), t_2 (1 - t_1) ) $
Clearly, point $G$ lies on $LT_3$. And this concludes the proof.