Dear sage-devel,

Applying full_simplify() to gamma functions sometimes causes an
error.  For example, this is correct:

sage: gamma(4/3).full_simplify()
1/3*gamma(1/3)

but this fails:

sage: gamma(1/3).full_simplify()
boom!

The error occurs because full_simplify() here runs the following
Maxima command:

minfactorial(factcomb(makefact(gamma(1/3))));

This returns (-2/3)!  [in Maxima] which then gets converted to Sage as
factorial(-2/3).  Since the Sage function factorial() is only defined
on non-negative integers, this causes an error.

I have created a Trac ticket for this (Trac #9240):

http://trac.sagemath.org/sage_trac/ticket/9240

and uploaded a patch to fix it.

Best,

Tom

---
Tom Coates
Royal Society University Research Fellow
Reader in Pure Mathematics
Imperial College London

-- 
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

Reply via email to