Hi Andrej,
On Fri, Jan 15, 2010 at 5:43 PM, andrejv <andrej.vodopi...@gmail.com> wrote:
<SNIP>
> d2 is defined in the testsuite for the Zeilberger algorithm. It is not
> necessary to load the tests, in share/contrib/Zeilberger/
> zeilberger.mac remove the last line which loads them.
Thank you very much for this pointer. I have created an updated Maxima
spkg [1]. Testing this updated spkg with Sage 4.3.1.alpha2, I see that
the problem originally reported in this thread is solved:
[mv...@mod sage-4.3.1.alpha2-maxima]$ ./sage
----------------------------------------------------------------------
| Sage Version 4.3.1.alpha2, Release Date: 2010-01-13 |
| Type notebook() for the GUI, and license() for information. |
----------------------------------------------------------------------
**********************************************************************
* *
* Warning: this is a prerelease version, and it may be unstable. *
* *
**********************************************************************
sage: var("a, b, c, d, e, f, g, h");
sage: S = matrix([[a,b], [c,d]])
sage: D = matrix([[e,f], [g,h]])
sage: K = S*D
sage: y = expand(K.det()); y
a*d*e*h - a*d*f*g - b*c*e*h + b*c*f*g
sage: simplify(y)
-a*d*f*g + a*d*h*e + b*c*f*g - b*c*h*e
sage:
sage:
sage: reset()
sage: var("a1, b1, c1, d1, a2, b2, c2, d2");
sage: A = matrix([[a1,b1], [c1,d1]])
sage: B = matrix([[a2,b2], [c2,d2]])
sage: M = A*B
sage: y = expand(M.det()); y
a1*a2*d1*d2 - a1*b2*c2*d1 - a2*b1*c1*d2 + b1*b2*c1*c2
sage: simplify(y)
a1*a2*d1*d2 - a1*b2*c2*d1 - a2*b1*c1*d2 + b1*b2*c1*c2
sage: bool(simplify(y) == y)
True
sage:
sage: var("a_1, b_1, c_1, d_1, a_2, b_2, c_2, d_2");
sage: A = matrix([[a_1, b_1], [c_1, d_1]])
sage: B = matrix([[a_2, b_2], [c_2, d_2]])
sage: M = A*B
sage: y = expand(M.det()); y
a_1*a_2*d_1*d_2 - a_1*b_2*c_2*d_1 - a_2*b_1*c_1*d_2 + b_1*b_2*c_1*c_2
sage: simplify(y)
a_1*a_2*d_1*d_2 - a_1*b_2*c_2*d_1 - a_2*b_1*c_1*d_2 + b_1*b_2*c_1*c_2
sage: bool(simplify(y) == y)
True
sage: d2 = var('d2')
sage: d2.simplify()
d2
However, even this updated Maxima spkg doesn't fix the problem
reported by William:
sage: factorial = var('factorial')
sage: factorial^3 + factorial/7 - 1
factorial^3 + 1/7*factorial - 1
sage: (factorial^3 + factorial/7 - 1).simplify()
BOOM!
<SNIP>
TypeError: no canonical coercion from <class
'sage.functions.other.Function_factorial'> to Symbolic Ring
[1]
http://boxen.math.washington.edu/home/mvngu/spkg/standard/maxima/maxima-5.20.1.p0.spkg
--
Regards
Minh Van Nguyen
--
To post to this group, send an email to sage-devel@googlegroups.com
To unsubscribe from this group, send an email to
sage-devel+unsubscr...@googlegroups.com
For more options, visit this group at http://groups.google.com/group/sage-devel
URL: http://www.sagemath.org
