Status: New Owner: ---- Labels: Type-Defect Priority-Medium
New issue 2184 by [email protected]: Wrong rref() http://code.google.com/p/sympy/issues/detail?id=2184 The code: ############################### CODE ############################### import sympy as sy qc = sy.Symbol('qc') qs = sy.Symbol('qs') M=sy.Matrix(([2*qc-qs,-qc,qs,-qc],[-qc,2*qc-qs,-qc,qs],[qs,-qc,2*qc-qs,-qc],[-qc,qs,-qc,2*qc-qs])) print M print "Determinant of M = ",M.det() print M.rref() ############################### returns this output: ############################### OUT ############################### [-qs + 2*qc, -qc, qs, -qc] [ -qc, -qs + 2*qc, -qc, qs] [ qs, -qc, -qs + 2*qc, -qc] [ -qc, qs, -qc, -qs + 2*qc] Determinant of M = 0 ([1, 0, 0, 0] [0, 1, 0, 0] [0, 0, 1, 0] [0, 0, 0, 1], [0, 1, 2, 3]) ############################### The rref seems to be wrong (I also checked with matlab and sage). -- 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.
