(38G) 40G: Matrix Power Programme

[Gerald H]

Edit: Amended to include negative powers.

Edit: Programme slightly more elegant.

Edit: Amended to avoid errors for some powers > 5*10^11.

A lack in the 38G & 40G (not the 40gs) is the inability to find integer powers ≠ 2 or -1 of matrices.

This programme supplies the lack.

Store the matrix to power in M1, enter the integer power on the stack in HOME view & activate the programme.

The powered matrix is recalled as Ans in HOME view & is also stored in M2.

M1 remains unchanged.

Ans►N:
IFTE(N<0,INVERSE(M1),M1)►M3:
ABS(N)►N:
SIZE(M1):
IDENMAT(Ans(1))►M2:
WHILE N
REPEAT
IF
N MOD 2
THEN
M3*M2►M2:
N-1►N:
END:
M3^2►M3:
N/2►N:
END:
M2:

[Francois Lanciault]

Thank you Gerald for these programs. Programs for these calculators are few and far between.

François

[Gerald H]

(03-15-2015 04:45 AM)Francois Lanciault Wrote:  Thank you Gerald for these programs. Programs for these calculators are few and far between.

François

Thanks for the recognition.

It would be nice if others published their programmes.

[Gerald H]

Just to advise of the latest edit to include negative powers.

[Eddie W. Shore]

Much gratitude, Gerald.

[Gerald H]

Just to advise of latest edit, avoiding errors for some powers > 5*10^11.