import numpy as np

A = np.array([[2,1,1],[1,2,1],[1,1,2]])
v = np.random.randint(100, size=3)
print(v)

for t in range(20):
	v = np.matmul(A,v) / np.linalg.norm(np.matmul(A,v))
	print(v, np.linalg.norm(v))

print(1/np.sqrt(3))