       Re: integral

• To: mathgroup at smc.vnet.net
• Subject: [mg46603] Re: integral
• From: Paul Abbott <paul at physics.uwa.edu.au>
• Date: Wed, 25 Feb 2004 13:06:58 -0500 (EST)
• Organization: The University of Western Australia
• References: <c1ehlu\$kfu\$1@smc.vnet.net>
• Sender: owner-wri-mathgroup at wolfram.com

```In article <c1ehlu\$kfu\$1 at smc.vnet.net>,

> I need to evaluate the the following integral:
>
> Integrate[Exp[-a^2*(x-x0)^2] BesselJ[1,bx] Sin[xY]}]
> x0, b and y are fixed
>
> Does anyone know, if there is a primitive to it ?

First note that you need to include a space between the b and x and x
and Y. Secondly, capital letter variables are generally not a good idea
(especially C, D, E, I).

Mathematica can compute the following integral:

int[a_,x0_][y_] = Integrate[x Exp[-a^2 (x-x0)^2] Sin[x y], x]

Since BesselJ[1,z] is equivalent to

z/2 Hypergeometric0F1[2, -z^2/4] // FunctionExpand

the general term in the series expansion of BesselJ[1,z] is

z/2 (-z^2/4)^n/n!^2/(n + 1),

and expansion in odd powers of z. Furthermore, since

D[Sin[x y], {y, 2}]

is

(-x^2) Sin[x y]

you can obtain the integral you want as an infinite sum of (parametric)
derivatives with respect to y of int[a,x0][y].

However, even the basic integral is rather complicated, and computing
all the required derivatives is unlikely to simplify things.

Instead, if you want to compute this integral over a specified domain
then using NDSolve (instead of NIntegrate) is probably the best approach:

With[{a=1, x0=2, b=3, y=1},
NDSolve[{f==0,f'[x]==Exp[-a^2 (x-x0)^2] BesselJ[1, b x] Sin[x y]},
f,{x,0,16}]]

This returns the antiderivative as an InterpolatingFunction.

Cheers,
Paul

--
Paul Abbott                                   Phone: +61 8 9380 2734
School of Physics, M013                         Fax: +61 8 9380 1014
The University of Western Australia      (CRICOS Provider No 00126G)
35 Stirling Highway
Crawley WA 6009                      mailto:paul at physics.uwa.edu.au
AUSTRALIA                            http://physics.uwa.edu.au/~paul

```

• Prev by Date: Re: Inversion using Cholesky Decomposition
• Next by Date: RE: RE: Re: corrected RE: Re: Computing sets of equivalences