On 17 dic, 11:48, "ma...@mendelu.cz" <ma...@mendelu.cz> wrote: > And another observation: > > maxima returns answer immediatelly (with a lag necessary to start > maxima) > m is the original matrix from x.py > > sage: m._maxima_().determinant().expand().sage() > x0^2*x2^2*x3^2*x7^2 - 2*x0*x1*x2*x3*x4*x5*x6*x7 + x1^2*x4^2*x5^2*x6^2 > > Anyway, the answer is different from expected one. I do not konw which > one is correct. > > Robert
If you perform the computations over the polynomial ring, sage gives an answer sage: R=FractionField(PolynomialRing(GF(2),",".join(map(genVar,range (0,10))))) sage: n=matrix(R.base(),m) sage: n-m [0 0 0 0 0 0 0 0] [0 0 0 0 0 0 0 0] [0 0 0 0 0 0 0 0] [0 0 0 0 0 0 0 0] [0 0 0 0 0 0 0 0] [0 0 0 0 0 0 0 0] [0 0 0 0 0 0 0 0] [0 0 0 0 0 0 0 0] sage: det(n) x1^2*x4^2*x5^2*x6^2 + x0^2*x2^2*x3^2*x7^2 However, the result differs from the one in maxima and maple About the result, the matrix is mod 2, so the expected anwer in maxima equals the one computed above sage: factor(det(n)) (x1*x4*x5*x6 + x0*x2*x3*x7)^2 -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support+unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sage-support URL: http://www.sagemath.org
