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