FUNMV - computes trigonometric and hyperbolic matrix functions times vectors
[C,S,s,m,mv,mvd] = FUNMV(t,A,B,FLAG,[],PREC) computes f(tA)B or f(tSQRTM(A))B, where f is a trigonometric or hyperbolic matrix function and SQRTM denotes any matrix square root, without explicitly forming f(tA) or SQRTM(A). PREC is the required accuracy, 'double', 'single' or 'half', and defaults to CLASS(A).
Details on the underlying algorithms can be found in
A. H. Al-Mohy, "A Truncated Taylor Series Algorithm for Computing the Action of Trigonometric and Hyperbolic Matrix Functions
Read More: https://epubs.siam.org/doi/abs/10.1137/17M1145227",