(49G) Function of a Matrix

[Gilles59]

Hi Gerald

I think you can do this with the integrated DIAGMAP function. With your example :

Code:

[[ -9 0 -1 -4 ]
[ 4 9 3 3 ]
[ 0 4 -9 9 ]
[ 5 -1 -6 -4 ]]

'LN(X)' DIAGMAP XNUM

takes ~11,6s in approx mode on a 50g (i dont test with the 49g)

Then :

Code:

 'EXP(X)' DIAGMAP XNUM

and your retrieve exactly the real parts and x.xxxxE-12 for the imaginary parts.