USAGE y = thfm ( x, MODE ) trigonometric/hyperbolic functions of square matrix x MODE cos sin tan sec csc cot cosh sinh tanh sech csch coth acos asin atan asec acsc acot acosh asinh atanh asech acsch acoth sqrt log exp NOTE --- IMPORTANT --- This algorithm does NOT USE an eigensystem similarity transformation. It maps the MODE functions to functions of expm, logm and sqrtm, which are known to be robust with respect to non-diagonalizable ('defective') x EXA thfm( x ,'cos' ) calculates matrix cosinus EVEN IF input matrix x IS NOT DIAGONALIZABLE ASSOC expm, sqrtm, logm, funm Copyright (C) 2001 Rolf Fabian <fabian@tu-cottbus.de> 010213 published under current GNU GENERAL PUBLIC LICENSE