We consider the problem of interpolation of matrix functions in case of ordinary, Jordan, Hadamard and Frobenius multiplication rule. We give interpolation formulas of Lagrange and Hermite type, obtain sufficient conditions for convergence of interpolation processes and estimate interpolation errors in the class of analytic functions.