わらｓぎ９９９＋ - 経済統計練習帳

import numpy as np
y=[20,40]
a=np.identity(2)
a[0][0]=4
a[1][0]=12
a[0][1]=8
a[1][1]=18
inp=np.identity(2)
for i in range(2):
    for j in range(2):
        inp[i][j]=a[i][j]/y[j]
 
res=np.identity(2)
res=res-inp
rev=np.linalg.inv(res)
print(rev)