>Split up the integral (or integrand) to compute the part near x=0
>using the asymptotic form of j(x,n) for small x to avoid underflow.
>Or disable the underflow error if it doesn't affect the final results.
>
>--
>Brian Gough
Hi. Thanks for your reply. I implemented a solution equivalent to yo
Jorge Talamantes writes:
> Dear all,
>
> I am trying to integrate the following function:
>
> I = \int_0^{x1} M (D, alpha, x, n) dx,
>
> where D, alpha and n are parameters to be passed to M, and
>
> M = x^(D-alpha-1) * [ j(x,n+0.5) ]^2.
>
> Here, j is the Bessel function of or
Dear all,
I am trying to integrate the following function:
I = \int_0^{x1} M (D, alpha, x, n) dx,
where D, alpha and n are parameters to be passed to M, and
M = x^(D-alpha-1) * [ j(x,n+0.5) ]^2.
Here, j is the Bessel function of order (n + 0.5).
For some combinations of D and alpha, the