Why doesn't Taylor expansion work for integrating to infinity?

  • Thread starter Thread starter pivoxa15
  • Start date Start date
  • Tags Tags
    Integral Weird
pivoxa15
Messages
2,250
Reaction score
1
Consider integrating sin(t)/t from 0 to infinity via the Taylor expansion method. If you do that then you will get infinite terms in the series which clearly is wrong as that integral equals pi/2. Can anyone tell me why using the Taylor expansion method dosen't work? It always dosen't work when integrating to infinity. But other methods like complex integral does work. Why?

Taylor expansion use to be a sure thing for me every time but now my innocence has been broken.
 
Physics news on Phys.org
It's probably a uniform convergence thing - or you're using an invalid Taylor expansion - or possibly both. In general you can't interchange limits and integrals.

What does
"It always dosen't work when integrating to infinity"
mean?
 
Just renormalize...

If you're brave, you'll just work out how to sum divergent series. I believe some little known physicist once did that for a special and got some award out of it or something. :-p
 
matt grime said:
It's probably a uniform convergence thing - or you're using an invalid Taylor expansion - or possibly both. In general you can't interchange limits and integrals.

What does
"It always dosen't work when integrating to infinity"
mean?

Taylor expansion expands a function in terms of polynomials. Integrate any polynomial from 0 to infinity and you get infinity hence "It always dosen't work when integrating to infinity". So a series with infinity in every term is not pleasant to say the least. In the case of sint/t, the terms are oscillating between negative and positive so the integration of the taylor polynomial on face value is 0 as you might like to subtract infinities after grouping the terms together, two by two. However infinities shouldn't be subtracted to the integral is not 0 but pi/2 when evaluated using another method.
 
You don't evaluate integrals over the range 0 to infinity just by putting in infinity and 0.
 
To the OP: What's the radius of convergence for the MacLaurin series expansion of \frac{\sin x}{x} ?

Daniel.
 
genneth said:
If you're brave, you'll just work out how to sum divergent series. I believe some little known physicist once did that for a special and got some award out of it or something. :-p

what are you talking about? dimensional regularization?
 
dextercioby said:
To the OP: What's the radius of convergence for the MacLaurin series expansion of \frac{\sin x}{x} ?

Daniel.

infinity. Or for all values of x. That is after we remove the singularity at x=0.
 
matt grime said:
You don't evaluate integrals over the range 0 to infinity just by putting in infinity and 0.

Its got to do with putting the limit as x->infinity. You are suggesting the limit doesn't exist?
 
  • #10
ice109 said:
what are you talking about? dimensional regularization?

I was thinking of Feynman and QED renormalisation...
 
  • #11
I think

<br /> \lim_{R\to\infty} \lim_{N\to\infty}\int\limits_{-R}^{R} \Big(\sum_{k=0}^{N} a_k x^k\Big) dx<br />

should be giving correct answer. But it is probably not possible to calculate it reasonably.
 
  • #12
Indeed -- after taking the N->inf limit the function is what Mathematica called SinIntegral, i.e. the integral of sinc. It's a special function, so no expression in terms of "normal" functions. It's pretty much for this purpose of calculating these random integrals that complex contour integration was pursued.
 
  • #13
there is another way to do this integral (that does not involve contour integration). define

<br /> f(x, \alpha) = \frac{\sin (x)}{x} e^{-\alpha x}<br />

where you recover your function when \alpha = 0. then differentiate past the integral sign with respect to \alpha. i'll leave the rest to you.
 
  • #14
quetzalcoatl9 said:
there is another way to do this integral (that does not involve contour integration). define

<br /> f(x, \alpha) = \frac{\sin (x)}{x} e^{-\alpha x}<br />

where you recover your function when \alpha = 0. then differentiate past the integral sign with respect to \alpha. i'll leave the rest to you.

Good idea, but it still doesn't give the right value for the integral.



\frac{d}{d\alpha }\int_{0}^{\infty }\frac{\sin x}{x}e^{-\alpha x}dx=\int_{0}^{\infty }\left( -\left( \sin x\right) e^{-\alpha x}\right) \,dx= \frac{1}{\alpha ^{2}+1}\lim_{x\rightarrow \infty }\left( e^{-\alpha x}\cos x+\alpha e^{-\alpha x}\sin x-1\right) =-\frac{1}{\alpha ^{2}+1}

\int_{0}^{\infty }\frac{\sin x}{x}e^{-\alpha x}dx=-\int \frac{d\alpha }{\alpha ^{2}+1}=\allowbreak -\arctan \alpha +C
 
  • #15
It does. The constant turns out to be pi/2, but the way I do that is not very rigorous...
 
  • #16
Ok Here it goes: Ignoring x between -1 and 1 :( maybe someone can fill that in later for me, but anyway: The original integral is less than the negative of its derivative, since its all values in the integrand are being divided by x.

So basically it becomes \frac{1}{\alpha^2 +1} &gt; -\arctan \alpha + C

Take the limit as alpha approaches infinity, the LHS approaches infinity, hence
\lim_{\alpha \to \infty} -\arctan \alpha + C = 0.

Knowing that as x approaches pi/2 on the left, it goes to infinity, the limit of x as x goes to infinity for arctan x is pi/2.

Ie C = pi/2.

I know there are many holes in the argument, but it works.
 
  • #17
that is correct. good job.
 
  • #18
to be clear:

f(x, \alpha) = -\tan^{-1} (\alpha) + C
\lim_{\alpha \rightarrow \infty} f(x, \alpha) = \lim_{\alpha \rightarrow \infty} \int_0^{\infty} \frac{\sin(x)}{x} e^{-\alpha x} dx = 0

but also

\lim_{\alpha \rightarrow \infty} f(x, \alpha) = -\tan^{-1} (\alpha) + C = -\frac{\pi}{2} + C
\Rightarrow C = \frac{\pi}{2}
\Rightarrow f(x, \alpha) = \frac{\pi}{2} - \tan^{-1} (\alpha)

and

f(x, 0) = \frac{\pi}{2}
 
  • #19
I see it now, nice job, quetzalcoatl9 !
 
Back
Top