Yes, I have simply tweaked log_integral, to get Li and it is the linear term that is my headache. I can't see Li in any of the integrators, but I can get a valid integration if I return Li as li(x)-li(2) in _eval(). But to me that seems untidy. It would be nice if _eval returned Li(x) but pass to integrate() li(x)-li(2). I'll post code to trac when I get an account set up. The #11143 patches I downloaded, applied ok on 5.2.rc02, I'll see if there are any updates but given I am a newbie with Mercurial I may go horribly wrong.
On Friday, July 27, 2012 9:24:35 PM UTC+1, kcrisman wrote: > > > > I don't know that this is something we have - integrate just goes > through Maxima (or Sympy or whatever if you choose a different > algorithm). Pynac doesn't support integrate methods. I suppose if > one were sufficiently clever one could do this (_integrate_), but that > shouldn't hold up #3401, which should be much easier now that #11143 > is ready to go in - that is, if Volker or someone can verify which > patches should actually be applied there ;-) > > Also, how would the integral of Li be different from the integral of > li - by a linear term, presumably? If we send it to Maxima > appropriately, this would probably solve the issue (not because Maxima > has both Li and li, but because we could perhaps get around that). > Why don't you put up your code at #3401 and we'll see if there is a > way to massage this into shape! > -- -- 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
