Как упростить матричные выражения в SymPy?

Рассмотрим следующий пример

import sympy as sy
n = sy.symbols('n')
A = sy.MatrixSymbol("A",n,n)
B = sy.MatrixSymbol("B",n,n)
C = sy.MatrixSymbol("C",n,n)
M = A.inverse()*B.inverse() - A.inverse()*C*B.inverse()
B.inverse()*M.inverse()*A.inverse()

В этом примере B^-1*(A^-1*B^-1 - A^-1*C*B^-1)^-1*A^-1.

Может ли SymPy упростить выражение для (IC)^-1? Если нет, как насчет каких-либо промежуточных результатов, например, сбора общих факторов в M?

from sympy import * n = symbols('n') A = MatrixSymbol("A",n,n) B = MatrixSymbol("B",n,n) C = MatrixSymbol("C",n,n) M = A.inverse()*B.inverse() - A.inverse()*C*B.inverse() expression = B.inverse()*M.inverse()*A.inverse() # convert expression to string then simplify simplify_expression = simplify(str(expression)) pprint(simplify_expression)

Ответ 1