Hi,
I am trying to numerically evaluate the integral of
1/(x^2-1)
from x=minf to x=inf. I tried romberg integration, eg.
g(x):=(mode_declare(x,float),1/(x^2-1))$
translate(g)$
rombergit:15$
rombergtol:1.e-4$
romberg(g,-5,5);
but I can't get it to converge. I played with romergit, rombergtol, and
rombergabs but there doesn't seem to be a parameter combination that solves
the integral in a reasonable period of time on my computer. Am I doing
something wrong or is romberg just inappropriate for the above type of
integral?
I also found a function quanc8 for adaptive integration in the
documentation. But loading the corresponding package fails:
(C13) load("qq");
Load failed for
/home/felix/local/share/maxima/5.9.0/share/numeric/qq.lisp
-- an error. Quitting. To debug this try DEBUGMODE(TRUE);)
Any ideas how to get this to work? Or won't it be of any help for the above
problem?
BTW, of course the above integral can be solved analytically. The result is
0.
Felix