Puissance n d’une matrice carrée : matpow

matpow élève une matrice carrée à la puissance n en la jordanisant On tape :

matpow([[1,2],[2,1]],n)

On obtient :

[[((-1)^n+3^n)/2,(-(-1)^n+3^n)/2],[(-(-1)^n+3^n)/2,((-1)^n+3^n)/2]]

On a en effet :
jordan([[1,2],[2,1]]) renvoie :
[[1,-1],[1,1]],[[3,0],[0,-1]]