I would like to see a nice wrapper function for quadpack integration
methods. I have written a nintegrate function:
(%i1) load(nintegrate)$
(%i2) nintegrate(x, x, 0, 1);
(%o2) 0.5
(%i3) nintegrate(exp(-x^2), x, 0, inf);
(%o3) .8862269254527579
(%i4) nintegrate(x, x, 0, inf);
Numerical integration failed with error code 5
-- an error. To debug this try debugmode(true);
(%i5) nintegrate(x, x, 0, a);
(%o5) nintegrate(x,x,0,a)
(%i6) subst(a=1, %);
(%o6) 0.5
(%i7) nintegrate(nintegrate(x+y, x, 0, y), y, 0, 1);
(%o7) .5000000000000001
(%i8) plot2d(nintegrate(sin(t)/t, t, 0, x), [x,-10,10])$
It is the first time I tried to write a simplifying function so maybe
I missed something. I would like to see a function like nintegrate
defined in quadpack, it would make numerical integration much easier
to use in maxima (in my opinion anyway).
--
Andrej
-------------- next part --------------
A non-text attachment was scrubbed...
Name: nintegrate.lisp
Type: application/x-extension-lisp
Size: 1187 bytes
Desc: not available
Url : http://www.math.utexas.edu/pipermail/maxima/attachments/20070523/d92e2a49/attachment.bin