This is the mail archive of the
gsl-discuss@sources.redhat.com
mailing list for the GSL project.
Re: some problem with QAWF integration
- From: Brian Gough <bjg at network-theory dot co dot uk>
- To: xie qiong <xieq at yahoo dot com>
- Cc: gsl-discuss at sources dot redhat dot com
- Date: Sun, 7 Apr 2002 21:48:32 +0100 (BST)
- Subject: Re: some problem with QAWF integration
- References: <20020406210558.52869.qmail@web10005.mail.yahoo.com>
xie qiong writes:
> Hello there,
>
> I have an urgent question about QAWF function.
>
> I need to computer a semi-infinit integration
> \int_{0} ^{\infty} f(x)*sin(w*x)/x dx
> The integrand f(x) * sin(w*x) / x has a finite value
> limit(equal to w) at x = 0.0 but f(x) / x doesn't a
> finite value limit at x = 0.0.
>
> In this case, what can I do to make QAWF work?
Break the integration range into two parts (0,pi/w) (pi/w,inf) and
compute the first with a standard integrator.
Brian