(DM42) Matrix exponential

[Gil]

Thanks for your explanations.

As suggested, I used the EGV command and it works nice.

\<< EGV DUP SIZE OBJ\-> DROP \-> ev d
\<< DUP INV d IDN 1 d
FOR i { i i } ev i GET EXP PUT
NEXT SWAP * * \->NUM
\>>
\>>

Or

« EGV DUP SIZE OBJ—> DROP —> ev d
« DUP INV d IDN 1 d
FOR i { i i } ev i GET EXP PUT
NEXT SWAP * * —>NUM
»
»



With
[[ 1. 2. 3. ]
[ 4. 5. 6. ]
[ 7. 8. 9. ]], we get:
[[ 1118906.69941 1374815.06293 1630724.42645 ]
[ 2533881.04188 3113415.03136 3692947.02084 ]
[ 3948856.38436 4852012.99979 5755170.61523 ]].

By the way, if I get 3 values 10, 20 and 30, are there commands to build a Matrix
[10 0 0
0 20 0
0 0 30], without using the above loop n (n=3)?

Regards,
Gil