"""
import numpy as np
# Input matrix A
A = np.array([[0.3, 0.3, 0.2], [0.2, 0.1, 0.2], [0, 0.2, 0.2]])
d = np.array([1, 2, 3])
# Check if A is a square matrix
if A.shape[0] != A.shape[1]:
print("A is not a square matrix!")
else:
n = A.shape[0]
#Calculate Leontief inverse matrix
x = np.matmul(np.linalg.inv(np.identity(n) - A), d)
print(x)
"""
leontief inverse行列を計算する。
・経済学的な話題で、ある商品達を生産するのに必要な資源（商品）の量をずらりと並べたら、行列が出来る。
・そして、社会全体でのそれぞれの商品の需要の量を並べると、ベクトルになる。
・生産量と生産するための消費、そして、需要の三者の関係を方程式で繋げると、平衡状態における生産量が見えてくる。
・この時、Leontief inverse matrix なる逆行列が重要になる。
・入力出力解析と言われる分野らしい。
・実際は、貯蔵量とか季節変動性とか色々ありそうだ。
詳しくはこちら。