Ask Your Question
1

calculating residue with maxima_methods

asked 2012-07-07 01:34:55 -0500

Stefán Ingi gravatar image

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?

edit retag flag offensive close delete

2 answers

Sort by » oldest newest most voted
2

answered 2012-07-07 03:42:41 -0500

achrzesz gravatar image
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))
edit flag offensive delete publish link more

Comments

I didn't know Sympy had this... see also http://trac.sagemath.org/sage_trac/ticket/11210

kcrisman ( 2012-07-07 16:17:33 -0500 )edit
0

answered 2012-07-07 16:24:09 -0500

kcrisman gravatar image

updated 2012-07-07 16:24:24 -0500

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.

edit flag offensive delete publish link more

Comments

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.

Stefán Ingi ( 2012-07-24 00:18:18 -0500 )edit

Hmm, that's interesting. I'll add that to the bug report.

kcrisman ( 2012-07-24 04:52:33 -0500 )edit

Your Answer

Please start posting anonymously - your entry will be published after you log in or create a new account.

Add Answer

[hide preview]

Stats

Asked: 2012-07-07 01:34:55 -0500

Seen: 272 times

Last updated: Jul 07 '12