summary:We propose a simple method to obtain sharp upper bounds for the interpolation error constants over the given triangular elements. These constants are important for analysis of interpolation error and especially for the error analysis in the Finite Element Method. In our method, interpolation constants are bounded by the product of the solution of corresponding finite dimensional eigenvalue problems and constant which is slightly larger than one. Guaranteed upper bounds for these constants are obtained via the numerical verification method. Furthermore, we introduce remarkable formulas for the upper bounds of these constants