#範例12-1：高斯喬丹消去法最後化簡的梯形矩陣
#指令：M.rref()
from sympy import *
M = Matrix([
    [1,1,2,9],
    [2,4,-3,1],
    [3,6,-5,0]
    ])
M_echelon = M.echelon_form()
print('簡化的梯形矩陣的echelon form =\n', M_echelon)
M_simple = M.rref()
print('簡化的梯形矩陣的echelon form =\n', M_simple)
print(M_simple[0].col(0))
print(M_simple[0].col(1))
print(M_simple[0].col(2))