Integrate sqrt(n)*frac{1}{(1+x^2)^n}

[davyyao]

Hi, I don't know how to calculate the following integral of an even function

limit{n goes to ∞} int_{-∞}^{∞}frac{sqrt(n)}{(x^2+1)^n}=?

Can this calculus be done by complex integral?

I was thinking for a few days and no answer so far : (

Thanks!

 

[mfb]

Complex integrals look like a good approach for natural n.
With the right integration parts, and a residual at the right point...

As this looks like homework, I moved your thread to the homework section.

 

[jackmell]

Hi, I don't know how to calculate the following integral of an even function

limit{n goes to ∞} int_{-∞}^{∞}frac{sqrt(n)}{(x^2+1)^n}=?

Can this calculus be done by complex integral?

I was thinking for a few days and no answer so far : (

Thanks!

Hi davy,

You can calculate it by brute force using the Residue Theorem to first compute the integral as a function of $n$, then take the limit. You'll need some Gamma function relationships and know about the rising factorial function and it gets messy. Write it first as:

$$\int_{-\infty}^{\infty} \frac{1}{(1+z^2)^n}dz=\int_{-\infty}^{\infty} \left(\frac{1}{(z+i)^n}\right) \left(\frac{1}{(z-i)^n}\right)dz$$

And consider a contour integral over the upper half-disc. And just for now assume the integral over the half-circle arc goes to zero as it's size goes to infinity. Can prove it latter. Then the integral resolves to just computing the residue. However, if you know about the Residue Theorem, then you'll realize you have to compute the general expression for the n-1 derivative of the first expression:

$$\frac{d^{(n-1)}}{dz^{(n-1)}}(z+i)^{-n}$$

at $z=i$. Ok, that's a little tough but you can do it. First try it for say n=5, n=6, the derivatives are easy to compute and resolve to basically the rising factorial function and some other factors. Then use deduction to derive the general expression for any n. Let's see if you can get that far.

 

[davyyao]

Hi , Thanks for mfb and jackmell, I use the residue theorem can obtain final result
lim_{n goes to infinity}sqrt(n)*pi*2^(-2n+2)
As I see, I guest it goes to zero as n goes to to infinity.
However, I numerically plot this integrand and it looks delta-like function.
Maybe this integral is finite...I am not sure~~

 

[Ray Vickson]

Hi, I don't know how to calculate the following integral of an even function

limit{n goes to ∞} int_{-∞}^{∞}frac{sqrt(n)}{(x^2+1)^n}=?

Can this calculus be done by complex integral?

I was thinking for a few days and no answer so far : (

Thanks!

If you change variables, first to y = x^2+1, then to u = 1/y, your integral from 0 to ∞ (which is 1/2 of what you want) reduces to a Beta function B(a,b), for some a and b defined in terms of n. The Beta can, in turn be re-expressed in terms of the Gamma function, using the identity
$$B(a,b) = \frac{\Gamma(a) \Gamma(b)}{\Gamma(a+b)}$$
You end up with a ratio of Gamma functions involving n, multiplied by some computable constants. Now you can take n → ∞ and get a fairly simple limit (using, for example, Stirling's formula to get limits of the Gamma functions for large n).

 

[jackmell]

Hi , Thanks for mfb and jackmell, I use the residue theorem can obtain final result
lim_{n goes to infinity}sqrt(n)*pi*2^(-2n+2)
As I see, I guest it goes to zero as n goes to to infinity.
However, I numerically plot this integrand and it looks delta-like function.
Maybe this integral is finite...I am not sure~~

No and no. First it's not hard to just get a ball-park idea of where the limit "might" be going by just numerically integrating say the first 500 terms and plotting them:

mytable = Table[Sqrt[n]*NIntegrate[1/(1 + z^2)^n, 
      {z, -Infinity, Infinity}], {n, 1, 500}]; 
mylimit = Graphics[{Red, Line[{{0, Sqrt[Pi]}, 
       {500, Sqrt[Pi]}}]}]; 
Show[{ListPlot[mytable, Joined -> True], mylimit}, 
  PlotRange -> All]

That's Mathematica code and if you don't know it then you have got to go out of your way to find a machine running it and learn it and practice coding it if you want to do math in the 21st century. Anyway, the plot does not seem to be going to zero but that's no guarantee it's not. Remember, just a rough idea.

Secondly, why don't we just then flat-out compute it in Mathematica then figure out how to arrive at that answer. That's not cheating.

In[38]:=
myintegrand = Integrate[1/(1 + z^2)^n, 
   {z, -Infinity, Infinity}]
Limit[Sqrt[n]*myintegrand[[1]], n -> Infinity]

Out[38]=
ConditionalExpression[(Sqrt[Pi]*Gamma[-(1/2) + n])/
   Gamma[n], Re[n] > 1/2]

Out[39]=
Sqrt[Pi]

Alright, fine, we got the answer. Now my third point depends if you just want the answer or if you want to learn how to get good in math. If it's the later, then you have got to try really hard to either muscle though my way and arrive at it or Ray's way even if it takes you a week to do so.

 

[Ray Vickson]

No and no. First it's not hard to just get a ball-park idea of where the limit "might" be going by just numerically integrating say the first 500 terms and plotting them:

mytable = Table[Sqrt[n]*NIntegrate[1/(1 + z^2)^n, 
      {z, -Infinity, Infinity}], {n, 1, 500}]; 
mylimit = Graphics[{Red, Line[{{0, Sqrt[Pi]}, 
       {500, Sqrt[Pi]}}]}]; 
Show[{ListPlot[mytable, Joined -> True], mylimit}, 
  PlotRange -> All]

That's Mathematica code and if you don't know it then you have got to go out of your way to find a machine running it and learn it and practice coding it if you want to do math in the 21st century. Anyway, the plot does not seem to be going to zero but that's no guarantee it's not. Remember, just a rough idea.

Secondly, why don't we just then flat-out compute it in Mathematica then figure out how to arrive at that answer. That's not cheating.

In[38]:=
myintegrand = Integrate[1/(1 + z^2)^n, 
   {z, -Infinity, Infinity}]
Limit[Sqrt[n]*myintegrand[[1]], n -> Infinity]

Out[38]=
ConditionalExpression[(Sqrt[Pi]*Gamma[-(1/2) + n])/
   Gamma[n], Re[n] > 1/2]

Out[39]=
Sqrt[Pi]

Alright, fine, we got the answer. Now my third point depends if you just want the answer or if you want to learn how to get good in math. If it's the later, then you have got to try really hard to either muscle though my way and arrive at it or Ray's way even if it takes you a week to do so.

You don't need to learn Mathematica to do math in the 21st century; you can use Maple instead, which is how I got the answer for this problem (while hiding the solution details from the OP).

 

[mfb]

A general remark: In homework questions (and only there), you will rarely see a factor like sqrt(n) with no purpose. If this does not change the convergence, the question would likely have a 1 there. This is a very strong hint that the the limit is probably a non-zero real number, with the sqrt-term designed to counter the running value of the integral.