(DM42) Matrix exponential

[Gil]

Thanks, Gjermund.

I changed the program, using "your" DIAG—> command instead of an Identity Matrix.

Depending on the size of the Matrix, the execution time now might be much shorter than the previous version.

\<< "
1 Arg
:

MAT
[n x n]" DROP EGV DUP SIZE OBJ\-> DROP \-> ev d
\<< DUP INV 1 d
FOR i ev i GET EXP
NEXT d \->ARRY d DIAG\-> SWAP * * \->NUM
\>>
\>>

Or
« "
1 Arg
:

MAT
[n x n]" DROP EGV DUP SIZE OBJ—> DROP —> ev d
« DUP INV 1 d
FOR i ev i GET EXP
NEXT d —>ARRY d DIAG—> SWAP * * —>NUM
»
»

Gil