>>> np.linalg.matrix_power(A,3)
array([[-54,  63, 117],
       [-36,  36,  72],
       [-45,  63, 126]])