Valery Pipin <pip@iszf.irk.ru> writes:
> On Thursday 20 November 2003 13:14, Raymond Toy wrote:
> >I've played around a bit with quadpack and have a very rough interface
> >to quadpack for maxima. Quadpack is an adaptive numerical quadrature
> >package.
> The maxima's qromb is slow and give wrong answers on some integrand
There is also share/numeric/qq.lisp, which is described in the manual.
The CVS version doesn't work, but it has been fixed
https://sourceforge.net/tracker/?func=detail&atid=104933&aid=842863&group_id=4933
and it is quite a bit faster than romberg on certain problems (try
qq.dem). I will commit it soon.
It does not handle Ray's integral so gracefully (although the
singularities could easily be avoided by an obvious change of
variables, so I dare say that his example is a bit contrived ;-).
Wolfgang