import numpy as np
import numpy.linalg as alg


def ps (n,m):
    return 1+n*m+2**(n+m)+n*m*2**(n+m-2)

A = np.array([[ps(i,j) for j in range(3)] for i in range(3)])
B = np.array([ps(i,3) for i in range(3)])
X = alg.solve(A,B)
print('matrice=',A,'\n','second membre=',B,'\n','coeff du po=',X)