Re: [Help-gsl] Fwd: Modifying adaptive integration function

[Francesco Florian]

> Hello,
>
> What exactly are you trying to do that QAG cannot? Perhaps there is 
> an alternate way without rewriting the function. 

I have to compute (for many t ∈ [0,1]) exp(∫_0^t f(x) dx), where lim_{x→1} 
f(x)=-∞.
Any adaptive quadrature computing the integral will use many intervals near the 
singularity, since it is the greatest source of errors; however, the 
exponential actually makes most errors int that region negligible.

What about my other question? I.e. I noticed that `retrieve (workspace, &a_i, 
&b_i, &r_i, &e_i);' is called at both the beginning and the end of the `do.., 
while' iteration (lines 181 and 231), and the second call seems useless. Is 
that correct?

Thank you
-- 
Francesco Florian