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