>>>>> "Edwin" == Edwin Woollett <woollett at charter.net> writes:
Edwin> Sorry for the delay - I didn't see your original reply
Edwin> in the Mailing List until today - been rather
Edwin> over-committed...
No problem.
Edwin> I am getting a Lisp error message when I use your version:
Oops. I've uploaded a new version that fixes this.
Unfortunately, it still takes a very long time. My original test used
the default parameters and it converged very quickly. With your new
parameters, I gave up after a minute or so.
But I think the problem is with the integrand, 1/(1+sqrt(x)). I'm
guessing the infinite slope at 0 is going to give lots of trouble.
bromberg(g, 1/2, 1) returns instantly. bromberg(g,1/100,1) takes a
little bit longer.
Don't know if this will help or not, but it might be worthwhile to
implement tanh-sinh quadrature
(http://en.wikipedia.org/wiki/Tanh-sinh_quadrature). (I think Richard
has mentioned this before.)
Ray