The eigenvectors ($\psi _i$) of the following matrix are rather complicated $$ \rho =\left( \begin{array}{cccc} G_+ & k & k & G_- \\ k & L_+ & L_- & k \\ k & L_- & L_+ & k \\ G_- & k & k & G_+ \\ \end{array} \right) $$ So, I want to suppose that $ E_i $ is the eigenvalue of $\rho$ that corresponding to the eigencfunction $\psi _i$, then to evaluate ($\psi _i$) as a function of $E_i $.
Is there a way I can do this with Mathematica?