Two algorithms for computing the matrix cosine function, Jorge Sastre, Javier Ibáñez, Pedro Alonso, Jesús Peinado, Emilio Defez. Applied Mathematics and Computation, Volume 312, pp. 66-77, November 2017, Preprint.  Matlab code cosmtay.m.