Hello!
I'm new in Sagemath and I want to calculate the limit of the power of a given 3x3 diagonal matrix:
M=matrix([[0,5,10],[1/10,0,0],[0,1/2,0]]); #Given matrix
J=M.eigenmatrix_right()[0]; #Diagonalization
n=var('n');
limit(J^n,n=infinity) #The limit command fails =(
Is there any way to do this?
Trank you so much!