(DM42) Matrix exponential

[Gil]

I tried for the HP50G a very simple approach :

exp(M) = {(I+M/9000/9000)^9000} ^ 9000
(twice 9000, as powers for n>= 10000 are not allowed in HP50G).

However the results with the above matrixes are not so accurate.

Is there a simple approach to get better results?

Thanks for your help.

Regards,
Gil