Re: Hankel transform question

• To: mathgroup at smc.vnet.net
• Subject: [mg82437] Re: Hankel transform question
• From: danl at wolfram.com
• Date: Sat, 20 Oct 2007 05:52:16 -0400 (EDT)
• References: <ff9t5n\$6cu\$1@smc.vnet.net>

```On Oct 19, 4:25 am, Jim Rockford <jim.rockfo... at gmail.com> wrote:
> I'm getting some disparate (or at least different looking) results for
> a particular Hankel transform related integral in Mathematica 5.2
> versus Mathematica 6.01.
>
> The integral I'm dealing with is the order-V Hankel transform of a
> constant ( f(r) = 1, say )
>
> Define
> g[s_] = s BesselJ[V,w s]
>
> and I need  the output for
>
> Integrate[g[s],{s,0,Infinity}]
>
> Both versions of Mathematica complain about not being able to verify
> convergence.  I believe this integral should give back a Dirac delta
> function.  Here are the outputs:
>
> *****************
>    version 5.2
> *****************
> V/w^2      ( with V and w real and positive)
>
> *****************
>   version 6.01
> *****************
> Integrate::idiv:  Integral of s BesselJ[V,w s] does not converge on
> (0,Infinity) >>
>
> (1) First of all, what accounts for the differences in output between
> the two versions of Mathematica?
> (2) Second, I know that   Integrate[s BesselJ[0,s],{s,0,Infinity}]
> should give Delta(s), but Mathematica 5.2  gives back the answer 0,
> and Mathematica 6.01 again balks and says that the integral doesn't
> converge.
>
> Is there any way for me to get Mathematica to "look for" Delta
> functions, instead of trying to grind out a numerical integration?
>
> Thanks,
> Jim

What I get in the development Mathematica kernel is simply an
unevaluated Integrate.

In[12]:= Integrate[s*BesselJ[v,w*s], {s,0,Infinity},
Assumptions->{v>0,w>0}]

Out[12]= Integrate[s BesselJ[v, s w], {s, 0, Infinity},
Assumptions -> {v > 0, w > 0}]

But the integral is in fact divergent. To see this, check the series
expansion at infinity.

In[13]:= InputForm[Series[s*BesselJ[v,w*s], {s,Infinity,2},
Assumptions->{v>0,w>0}]]

Out[13]//InputForm=
Cos[Pi/4 + (Pi*v)/2 - s*w]*
SeriesData[s, Infinity, {Sqrt[2/Pi]/Sqrt[w]}, -1, 3, 2] +
SeriesData[s, Infinity,
{(-1 + 4*v^2)/(4*Sqrt[2*Pi]*w^(3/2))}, 1, 3, 2]*
Sin[Pi/4 + (Pi*v)/2 - s*w]

Notice that the "main" term is, up to scaling and phase change, of the
form Cos[s]*Sqrt[s]. This diverges in classical integration.

To get a regularized result, you can do as below.

In[14]:= InputForm[Integrate[s*BesselJ[v,w*s], {s,0,Infinity},
Assumptions->{v>0,w>0}, GenerateConditions->False]]

Out[14]//InputForm= v/w^2

Offhand I do not see any reason to expect a delta function result.
Certainly DiracDelta[s] would be incorrect, because s is the integral
variable of (definite) integration, hence cannot appear in the result.

Daniel Lichtblau
Wolfram Research

```

• Prev by Date: Re: A riddle: Functions that return unevaluated when they cannot
• Next by Date: Re: Ver 6.0 Subsets and Packed Array Question
• Previous by thread: Hankel transform question
• Next by thread: Re: Hankel transform question