Determination of eigenvectors with Lagrange multipliers

We present a method to determine the eigenvectors of an n× n Hermitian matrix by introducing Lagrange undetermined multipliers. In contrast to a usual Lagrange multiplier that is a number, we introduce matrix-valued multipliers with a constraint equation, which make the eigenvalue equation directly solvable. Then, there exists a unique solution for each eigenvalue equation and the eigenvectors are obtained by imposing the constraint limit. This method is in clear contrast to the conventional approach of Gaussian elimination and it will be a good pedagogical example for conceptual understanding for the gauge symmetry just with the knowledge of quantum physics and linear algebra at the undergraduate level.