Generalized Hermite formula for the two-sided Lagrange-Sylvester interpolation

The classical Hermite formula of the scalar polynomial interpolation is generalized to the setting of two-sided Lagrange-Sylvester interpolation for matrix polynomials. The formula then is applied to a description of the set of all matrix polynomial solutions to that two-sided interpolation.