Use Python to calculate P**n while n approach +oo and P is a matrix

The code is

代码是

import sympy as sm
import numpy as np

k=sm.Symbol('k')
p=np.matrix([[1./2,1./4,1./4],[1./2,0,1./2],[1./4,0,3./4]])
sm.limit(p**k,k,sm.oo)
import sympy as sm
import