sympy gradient of function of n-dim vector

I would simply like to define a symbolic function of this signature

Sympy can’t notice elementary algebraic identity: I – A*(A+I).inverse() == (A+I).inverse()

The following sympy code outputs False instead of True

Sympy can’t notice elementary algebraic identity: I – A*(A+I).inverse() == (A+I).inverse()

The following sympy code outputs False instead of True