Dear sage-devel, Calculating the Taylor expansion of gamma functions sometimes causes an error. For example, this is correct:
sage: taylor(gamma(1/2+x),x,0,2) 1/4*(pi^2 + 2*euler_gamma^2 + 8*euler_gamma*log(2) + 8*log(2)^2)*sqrt(pi)*x^2 - (euler_gamma + 2*log(2))*sqrt(pi)*x + sqrt(pi) but this fails: sage: taylor(gamma(1/3+x),x,0,2) boom! The error occurs because when parsing the output from Maxima we should replace all occurrences of the polygamma function psi[n](foo) by psi(n,foo), but currently we don't. I created a Trac ticket for this (Trac #9217): http://trac.sagemath.org/sage_trac/ticket/9217 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
