A = np.zeros((d,d))
 Pi = np.zeros(d)