[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [Help-gsl] Oscillating integral: Product of Bessel functions

From: Inigo Aldazabal Mensa
Subject: Re: [Help-gsl] Oscillating integral: Product of Bessel functions
Date: Thu, 9 Mar 2006 10:57:07 +0100
User-agent: KMail/1.8


As x->\infty your integrand diverges oscillating something like exponentially, 
as gsl_sf_bessel_In does, and nothing kills this behaviour. Are you sure this 
integral should converge?



El Jueves, 9 de Marzo de 2006 09:54, Mehmet Suzen escribió:
> Hello All,
> I have tried to use QAWO or QAWF for integrating function which is
> a product of modified bessel functions. Integral looks like
>  f = int_{0}^{infty} x * cos(-x*Z) * gsl_sf_bessel_In(m,x*R)  dx
> X and R are parameters.
> but I am getting  error message ;
> gsl: qawo.c:436: ERROR: cannot reach tolerance because of roundoff error
> I have tried different epsabs or TOL, it doesn't help. How do I choose
> right tolerance and bin sizes? And I have seen on the net that product of
> bessel functions are problematic in integration. Any comment would be
> helpful. Cheers,
> --Mehmet
> _______________________________________________
> Help-gsl mailing list
> address@hidden

reply via email to

[Prev in Thread] Current Thread [Next in Thread]