Subject: proposal to cut out romberg and quanc8 functions
From: Robert Dodier
Date: Wed, 3 Jan 2007 05:29:55 -0700
On 1/3/07, Andrej Vodopivec <andrej.vodopivec at gmail.com> wrote:
> I don't have any strong opinions about this. But romberg and quanc8
> have the advantage that they can be nested:
>
> (%i1) romberg(romberg(x^2+y^2, x, 0, y), y, 0, 1);
> (%o1) 0.33333333333333
> (%i2) integrate(integrate(x^2+y^2, x, 0, y), y, 0, 1);
> (%o2) 1/3
By itself it's not enough to keep romberg and quanc8. I don't think
we want to tell people, "You can simply nest the numerical integration
calls if you're OK with using methods which don't work so well."
Quoting the inner call yields the same result with quad_qags --
quad_qags ('(quad_qags (x^2 + y^2, x, 0, y)[1]), y, 0, 1);
=> [.3333333333333333, 3.70074341541719E-15, 21, 0]
best,
Robert