I am trying to find three independent infinitesimal transformations that leave my lagrangian invariant.
$$\mathcal{L}=-\frac{1}{2}\partial_{\mu}\phi_I\partial^{\mu}\phi^I-\frac{1}{2}m^2\phi_I\phi^I $$ The indecis $I=1,2,3$ represent three different scalar fields. I have assumed that the transformation is $$\phi^I\rightarrow (\mathbb{1}-i\epsilon_AT^{IJ}_A)\phi^J$$ where A is a label $A=1,2,3$ for the three different transformations, $\epsilon$ is infinitesimal and $T^{IJ}$ is a $3\times 3$ matrix.
I have assumed that I can drop higher order $\epsilon$ so after insertion I get that the variation of the lagrangian is $$\delta\mathcal{L}=\frac{1}{2}\partial_{\mu}\phi_I\partial^{\mu}i\epsilon T_A^{IJ}\phi^J+\frac{1}{2}\partial_{\mu}i\epsilon T_A^{IJ}\phi^I\partial^{\mu}\phi_J $$ which I want to be equal to zero. Now I have matrix equation but it is here that I run into trouble. I think I am confusing myself in trying to solve this matrix equation. I would appreciate some help in solving it and some conformation that my analysis/method is correct.
EDIT
To clarify my problem: we have a matrix equation with the matrix sizes $1\times2\cdot2\times2\cdot2\times1+2\times2\cdot2\times1\cdot1\times2=1\times1+2\times2=0$ and I do not see how to solve it.