On Sat, 20 Feb 2010 13:13:09 -0800 Robert Bradshaw <rober...@math.washington.edu> wrote:
> On Feb 20, 2010, at 12:40 PM, John H Palmieri wrote: > > > > I was curious about this, so I tried specifying the number of > > digits: > > > > sage: h = integral(sin(x)/x^2, (x, 1, pi/2)); h > > integrate(sin(x)/x^2, x, 1, 1/2*pi) > > sage: h.n() > > 0.33944794097891573 > > sage: h.n(digits=14) > > 0.33944794097891573 > > sage: h.n(digits=600) > > 0.33944794097891573 > > sage: h.n(digits=600) == h.n(digits=14) > > True > > sage: h.n(prec=50) == h.n(prec=1000) > > True > > > > Is there an inherit limit in Sage on the accuracy of numerical > > integrals? > > > I don't know if this above uses ginac, maxima, or gsl, but it seems > limited to double precision (which wouldn't surprise me for any of > those three systems). It's also very new. It should probably be > raising a NotImplemented error, but at least it's not casting the > result to a higher precision ring. This is now #8321: http://trac.sagemath.org/sage_trac/ticket/8321 Thanks. Burcin -- 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
