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: xie qiong <xieq at yahoo dot com>
- To: Brian Gough <bjg at network-theory dot co dot uk>
- Cc: gsl-discuss at sources dot redhat dot com
- Date: Mon, 8 Apr 2002 11:08:26 -0700 (PDT)
- Subject: Re: some problem with QAWF integration
Thanks first.
This method seems great but it breaks down when w
increases beyond 1e4.
Have you any suggestion?
Qiong
--- Brian Gough <bjg@network-theory.co.uk> wrote:
> 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
>
>
>
__________________________________________________
Do You Yahoo!?
Yahoo! Tax Center - online filing with TurboTax
http://taxes.yahoo.com/