ay([[5, 40,50],[50, 10,130],[15, 15,20]]) ys = np.array([100, 200,120]) print(a) inp=np.identity(3) for i in range(3): for j in range(3): inp[i][j]=a[i][j]/ys[j] res=np.identity(3) res=res-inp rev=np.linalg.inv(res) print(rev)