On July 29, 2009, Robert Dodier wrote:
> The quadpack functions don't need to be protected against
> unknown limits of integration (not since a few revisions ago).
> While use of lambda is good form in this context,
> the example can be simplified a little by omitting lambda's.
> (%i3) quad_qags (quad_qags (exp (x - y^2), y, 1, 2 + x)[1], x, 0, 1);
> (%o3) [.2384683615524653, 2.647530656401834e-15, 21, 0]
thanks for the heads-up on this flexibility of quad_qags, which I have not
been aware of.
> I can't think of a problem for which I'd recommend romberg over quadpack.
I will rewrite ch. 8 Sections to emphasize the attributes of quad_qags
compared
to romberg.
Ted