(DM42) Matrix exponential

[Werner]

(08-21-2023 11:49 AM)Albert Chan Wrote:  1118905.69945      1374815.06297      1630724.42651
2533881.04197      3113414.03144      3692947.02095
3948856.38452      4852012.99997      5755169.61548
(..)
HP71B, SCALE BY 2^16 or 2^17

-1.73575875814      0.551819099660
-1.47151759909      0.103638240710

I can reproduce these exactly on the 42S using the program below.
(well making sure I have P=13 in the first example..)
If you want MEXP, perform an extra M+I

00 { 99-Byte Prgm }
01▸LBL "MXP-1"
02 ENTER
03 FNRM
04 LN
05 2
06 LN
07 ÷
08 10
09 +
10 IP
11 X<0?
12 CLX
13 X<>Y
14 2
15 RCL ST Z
16 +/-
17 Y^X
18 ×
19 RCL ST X
20 4
21 ÷
22 XEQ 01
23 3
24 ÷
25 RCL× ST Y
26 XEQ 01
27 2
28 ÷
29 RCL× ST Y
30 XEQ 01
31 X<>Y
32 ×
33▸LBL 04
34 ENTER
35 +
36 LASTX
37 STOx ST X
38 +
39 DSE ST Y
40 GTO 04
41 RTN
42▸LBL "M+I"
43▸LBL 01 @ add I
44 EDIT
45 ^
46▸LBL 02
47 ↓
48 ENTER @ no, you can't remove this
49 1
50 +
51 →
52 FC? 77
53 GTO 02
54 EXITALL
55 END

Cheers, Werner