import numpy as np
import scipy
from scipy import linalg

# 1. 求解Ax = b
A = np.array([[3, 1], 
　　　　　　　　[1, 2]])
b = np.array([9, 8])
x = np.linalg.solve(A, b)
print(x)

# 2. 常用方法
A = np.array([[1, 2, 3], 
              [4, 5, 6], 
              [7, 8, 9]])
I = np.eye(3) # 單位矩陣
A.T # A的轉置
3 * A # 普通乘
np.dot(A, I) # 矩陣乘法
np.linalg.inv(A) # A的逆

# 3. A的LU分解（有時並不標準L）
L, U = scipy.linalg.lu(A, True)
print(L)
print(U)