> I think this is a nice convenience function. I like how it uses
> quad_qags for bounded intervals, and quad_qagi for unbounded. I
> rarely ever use any other integration routine.
>
> Some comments:
>
> o What desired precision is used?
> o Should the user be allowed to choose that?
> o What about the max number of intervals allowed for the integrator?
> o Should that be user selectable?
Currently they are the default from quadpack (epsrel=1e-8 and
lim=200). My preference for optional arguments is nintegrate(expr,
var, lo, hi, epsrel=1e-8, lim=200) and not that the meaning of an
optional argument is defined by position (epsrel the fifth argument
and lim the sixth like in quadpack). I can add this.
--
Andrej