(DM42) Matrix exponential

[John Keith]

(08-26-2023 10:20 PM)Albert Chan Wrote:  If n×n elements sized about the same, say k, FNORM(X) ≈ k*n

Code:

< 40 P=MAX(0,IROUND(LOG2(MAXAB(X)))+10) @ MAT X=(1/2^P)*X

> 40 P=MAX(0,IROUND(LOG2(FNORM(X)))+9) @ MAT X=(1/2^P)*X

I like this, especially since the HP-28/48 don't have MAXAB but do have ABS (same as FNORM ).