Approximation of functions of large matrices with Kronecker structure