Comment #1 on issue 2231 by [email protected]: Matrix determinant properties are not satisfied.
http://code.google.com/p/sympy/issues/detail?id=2231

    >>> B = Matrix([[a,b],[c,d]])


Try to use `expand()` method after multiplication. For me:

In [14]: (B*B).det()
Out[14]: -2*a*b*c*d + a**2*d**2 + b**2*c**2

In [15]: B.det()*B.det()
Out[15]: (a*d - b*c)**2

In [16]: lhs =  (B*B).det()

In [17]: rhs =  B.det()*B.det()

In [18]: lhs
Out[18]: -2*a*b*c*d + a**2*d**2 + b**2*c**2

In [19]: rhs
Out[19]: (a*d - b*c)**2

In [20]: rhs.expand()
Out[20]: -2*a*b*c*d + a**2*d**2 + b**2*c**2

In [21]: rhs.expand() == lhs
Out[21]: True


--
You received this message because you are subscribed to the Google Groups 
"sympy-issues" group.
To post to this group, send email to [email protected].
To unsubscribe from this group, send email to 
[email protected].
For more options, visit this group at 
http://groups.google.com/group/sympy-issues?hl=en.

Reply via email to