### Limit of a matrix power

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]]); ~~J=matrix([[1,0,0],[0,-1/2-1/2*I,0],[0,0,-1/2+1/2*I]]); #Given diagonal 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!