On Tue, 2004-04-27 at 12:35, Stavros Macrakis wrote:
> Detailed discussion of Gosse's numerical integration question below.
I respectfully claim that we have two totally different problems than
the ones discussed in detail by Stavros. Ray Toy added quadpack routines
to Maxima a while ago. quadpack is a great deal more sophisticated than
romberg.
--------------------------------------------------------
|abacus3>rmaxima
Maxima 5.9.0.1cvs http://maxima.sourceforge.net
Using Lisp CMU Common Lisp 18e fedora release 0.fdr.2.e.1
Distributed under the GNU Public License. See the file COPYING.
Dedicated to the memory of William Schelter.
This is a development version of Maxima. The function bug_report()
provides bug reporting information.
(C1) quad_qag(t**(10.84)*2.71828**(-t*t),t,0,1000,1);
(D1) [52.3756706393706, 1.099172453199315e-7, 375, 0]
(C2)
--------------------------------------------------------
So I think we have two problems:
1) The quadpack routines are entirely undocumented. One has to read
numerical/slatec/quadpack.lisp to find them.
2) There should probably be one higher-level interface to quadpack. I
don't want to get rid of the current low-level interfaces, but I think a
simpler starting point would be helpful.
--Jim