Comment #13 on issue 2130 by asmeurer: wrong limit http://code.google.com/p/sympy/issues/detail?id=2130
I don't like that. It is basically oo - oo (if you call .doit() for example, it will be converted to nan). I think we should either special case the code to return EulerGamma, or else just raise NotImplementedError as before.
A possibility btw would be to make the limit code work with polygamma(0, x) and then use the identity summation(1/k, (k, 1, n)) == polygamma(0, n + 1) + EulerGamma. Just a though; I don't actually know enough about the Gruntz algorithm to know how to best proceed with this.
-- You received this message because you are subscribed to the Google Groups "sympy-patches" group. To post to this group, send email to sympy-patches@googlegroups.com. To unsubscribe from this group, send email to sympy-patches+unsubscr...@googlegroups.com. For more options, visit this group at http://groups.google.com/group/sympy-patches?hl=en.
