import numpy as np
np.set_printoptions(precision=3)
np.set_printoptions(suppress=True)

A=np.array([[1,3,-2,-3],[-1,0,1,0],[0,3,-1,-3],[-1,-2,1,2]])

x=np.array([0,1,0,0])

x=np.matmul(A,x)
n=np.linalg.norm(x,np.infty)
print("Norma: %.3f" % n)
x=x/n
print(x)

np.linalg.eig(A)