v=vector(RR,[0,0,25]) A=matrix(RR,[[0.8,0.3,0.2],[0.1,0.5,0.6],[0.1,0.2,0.2]]) #v #A show(A) #for i in range(20): # v*A^i