On Friday, August 29, 2014 7:15:27 PM UTC+2, Robert Dodier wrote:
>
> QUADPACK ...
>
I've tried this in mpmath's quad. I think it works there, but maybe I've
overlooked the actual problem.
sage: import mpmath as mp
sage: f1 = lambda _ : 1. / mp.sqrt(_^3 + 2)
sage: f2 = lambda _ : 1. / mp.sqrt(_^4 + 2)
sage: mp.quad(f1, [1, 10^9])
mpf('1.8280370240753521')
sage: mp.quad(f1, [1, 10^10])
mpf('1.8280802692664053')
sage: mp.quad(f1, [1, 4.5 * 10^10])
mpf('1.8280908397677311')
sage: mp.quad(f1, [1, 10^11])
mpf('1.8280939410898087')
sage: mp.quad(f1, [1, 10^12])
mpf('1.8280982294340853')
sage: mp.quad(f2, [1, 10^9])
mpf('0.88156905940188868')
sage: mp.quad(f2, [1, 10^10])
mpf('0.88156905993967871')
sage: mp.quad(f2, [1, 4.5 * 10^10])
mpf('0.88156905860896162')
sage: mp.quad(f2, [1, 10^11])
mpf('0.88156905640783767')
-- H
--
You received this message because you are subscribed to the Google Groups
"sage-support" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to sage-support+unsubscr...@googlegroups.com.
To post to this group, send email to sage-support@googlegroups.com.
Visit this group at http://groups.google.com/group/sage-support.
For more options, visit https://groups.google.com/d/optout.
