# Créé par shemon, le 26/01/2026 en Python 3.7
from math import*
import numpy as np
I3=np.eye(3)
B=np.array([[-2,-4,1],[1,3,-1],[1,0,2]])
B2=np.dot(B,B)
B3=np.dot(B2,B)
PB=B3-3*B2-B+3*I3