Re: How to deal with this integral
- To: mathgroup at smc.vnet.net
- Subject: [mg97837] Re: How to deal with this integral
- From: m.r at inbox.ru
- Date: Mon, 23 Mar 2009 04:03:59 -0500 (EST)
- References: <gq2eva$edb$1@smc.vnet.net>
paperkite rainyday wrote:
> hi everybody,
>
> I want to deal with this integral with Mathematica:
>
> Integrate[a*b*Exp[-(a+b)]*DiracDelta[1/(4*a)+1/(4*b)+c-t],{a,0,Infinity},{b,0,Infinity},{c,0,Infinity},Assumptions->t>0],
>
> Mathematica 6.0 gives 0. But it is wrong, the correct result is
> (1/2)*x^2*Exp[-x]*((1+x)*BesselK[0,x]+(2+x+2/x)*BesselK[1,x]),
> where x=1/(2*t).
>
> What's wrong?
>
> Thanks,
>
>
> paperkite
It's not as fun when you know the answer, but you can evaluate the
integral by making the change of variables u == a + b, v == a b. The
(u, v) region is given by
In[1]:= Reduce[t > 0 && b > a > 0 && 1/a + 1/b < 4 t && a + b == u &&
a b == v,
{u, v}, {a, b}, Reals]
Out[1]= t > 0 && u > 1/t && u/(4 t) < v < u^2/4
(we're integrating over b > a, where the substitution function is one-
to-one).
In[2]:= Assuming[t > 0, 2 Integrate[v E^-u/Sqrt[u^2 - 4 v],
{u, 1/t, Infinity}, {v, u/(4 t), u^2/4}]] /.
HoldPattern@ MeijerG[{{}, {a_}}, {{b_, c_}, {}}, z_] :>
z^b E^-z HypergeometricU[a - c, b - c + 1, z] // FullSimplify
Out[2]= (E^(-(1/(2 t))) (BesselK[1, 1/(2 t)] + (1 + 2 t) BesselK[2, 1/
(2 t)]))/(16 t^3)
Maxim Rytin
m.r at inbox.ru