[Maxima] proposal to cut out romberg and quanc8 functions

Robert Dodier robert.dodier at gmail.com
Wed Jan 3 06:29:55 CST 2007


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


More information about the Maxima mailing list