The following used to work around six months ago on sagenb.org:
var('x')
var('t')
( x/((4*x - 1)*(t^2 - t + x)) ).maxima_methods().residue(t,(1-sqrt(1-4*x))/2)
but now it returns zero which is incorrect. Can anybody help me get the correct answer again?
sage: from sympy import Symbol, residue
sage: x = Symbol("x")
sage: t = Symbol("t")
sage: residue(x/((4*x - 1)*(t^2 - t + x)),t,(1-sqrt(1-4*x))/2)
-x/((-4*x + 1)**(1/2)*(4*x - 1))
I didn't know Sympy had this... see also http://trac.sagemath.org/sage_trac/ticket/11210
I get this in a very old version of Sage's Maxima, the current Sage Maxima, and the most recent one I have.
Maxima 5.27.0 http://maxima.sourceforge.net
using Lisp SBCL 1.0.24
Distributed under the GNU Public License. See the file COPYING.
Dedicated to the memory of William Schelter.
The function bug_report() provides bug reporting information.
(%i1) residue(( x/((4*x - 1)*(t^2 - t + x)) ),t,(1-sqrt(1-4*x))/2);
(%o1) 0
What/when was the answer you got before? I can't get this to do anything else in any version of Sage, either, not just the Maxima in it.
This is now Maxima artifact 3541292; I've also put it at the relevant Sage ticket.
I was reexamining my old code and found the culprit. When it worked before I had run "assume(x<1/8)". That makes maxima give the same answer as sympy. Sorry for the trouble, but at least I learned something about sympy.
Asked: 2012-07-07 08:34:55 +0100
Seen: 1,309 times
Last updated: Jul 07 '12
Copyright Sage, 2010. Some rights reserved under creative commons license. Content on this site is licensed under a Creative Commons Attribution Share Alike 3.0 license.
