Interpolation of cosine operator functions