Hola.
Hmm... no se si me estoy perdiendo de algo, pero en tu mensaje inicial no mencionas que

y

son simétricas, de hecho no mencionas que

sea un espacio con producto interno. Pero obviando ésto, si ya se puede diagonalizar simultaneamente a los operadores, como te dije en mi anterior mensaje y lo haz expuesto muy bien en tu anterior respuesta, todo se reduce al caso de dos matrices diagonales.
Ahora, justamente para ésto nos valemos de los polinomios de interpolación de Lagrange. Puedes comprobar que dados los

puntos del plano

, su polinomio de Lagrange de grado

(si los puntos son diferentes) asociado,

, es justamente un polinomio que verifica las igualdades

, para todo

. Éstas últimas igualdades las puedes comprobar a mano, justamente por la forma que tiene el polinomio asociado a la lista de puntos que tomamos. (Siendo coloquiales, digamos que éstas igualdades son justamente la "gracia" que tienen los polinomios de interpolación de Lagrange).
Entonces, debido a las anteriores igualdades, podemos verificar que
![[S]_{1}=L\big([T]_{1}\big)](/foros/Sources/latexrender/pictures/e6b1164d3e42fbc10d546dd06989890a.png)
. Y ésto terminaría la prueba, cualquier duda, pregunta.
Saludos.