3
$\begingroup$

Here $A>0$, $w$-real,$\mathtt{i}$-complex. Mathematica gives the answer: $$\frac{1}{2}e^{2\pi w}(\mathtt{i}\pi+2\Gamma(0,2\pi w)-2\Gamma(0,2(1+\mathtt{i}A)\pi w)+2\ln(-\mathtt{i}+A)+2\ln(w)-2\ln(w+iA)) $$ My questions: 1. How to obtain this results without mathematica? 2. Why does this integral grow exponentially as a function of $w$?

Mostly I don't understand why this integral explodes.The real part is:$$\int_{0}^{A}\frac{x\cos(2\pi wx)+\sin(2\pi wx)}{1+x^{2}}\ {d}x $$ and geometrically I don't understand why this integral grows so fast

  • 0
    What do you mean by Gamma[0,2πw]?2011-05-05
  • 0
    incomplete gamma function (upper)2011-05-05
  • 0
    @AD and others: For clarification the Mathematica Manual say $\Gamma(a,z_0,z_1)$ is the generalized incomplete gamma function $\Gamma(a,z_0)-\Gamma(a,z_1)$.2011-05-05
  • 0
    If $A$ were 1, this is surely easy to do... let me take a crack at it later. But wait, you really have $i=\sqrt{-1}$ in both the numerator and denominator?2011-05-05
  • 0
    yes, it is a root of my polynomial2011-05-05
  • 1
    For now, I'll only note that your results can be recast in terms of sine and cosine integrals; I'm peering at your problem now, and will see what comes up.2011-05-05
  • 1
    In particular, have a look at [this](http://functions.wolfram.com/ElementaryFunctions/Sin/21/01/02/02/02/) and [this](http://functions.wolfram.com/ElementaryFunctions/Sin/21/01/02/02/02/)...2011-05-05

1 Answers 1

2

You have a typo in your post. The first exponential in the answer given by Mathematica should read $e^{-2\pi w}$. In total the solution does NOT grow exponentially. The exponential growth of the $\Gamma$-function is "cured" by the factor $e^{-2\pi w}$.

To show that the integral does not grow exponential, it is possible to obtain an asymptotic expression for $w\to\infty$ without resorting to the explicit expression given in your post. To this end, we use integration by parts $$\int_{0}^{A}\frac{e^{-2\pi iwx}}{x-i}dx = \frac{ie^{-2\pi iwx}}{2\pi w(x-i)} \Biggr|_{x=0}^A + \frac{i}{2\pi w} \int_0^A \frac{e^{-2\pi iwx}}{(x-i)^2} dx.$$ Repetitive integration by parts yields an asymptotic expansion for $w\to\infty$. The boundary term is the leading term, thus we have $$\int_{0}^{A}\frac{e^{-2\pi iwx}}{x-i}dx \sim \frac{1}{2\pi w} \left( 1- \frac{e^{-2\pi i w A}}{1+i A}\right).$$

For fun, I also give the next order term (if one is only interested in large $w$ this asymptotic expansion may prove more useful than the exact expression in terms of not so elementary functions). The next integration by parts yields $$\int_{0}^{A}\frac{e^{-2\pi iwx}}{x-i}dx=\frac{ie^{-2\pi iwx}}{2\pi w(x-i)} \Biggr|_{x=0}^A - \frac{e^{-2\pi iwx}}{[2\pi w(x-i)]^2} \Biggr|_{x=0}^A- \frac{1}{(2\pi w)^2} \int_0^A \frac{e^{-2\pi iwx}}{(x-i)^3} dx,$$ which gives the next term in the asymptotic expansion $$\int_{0}^{A}\frac{e^{-2\pi iwx}}{x-i}dx \sim \frac{1}{2\pi w} \left( 1- \frac{e^{-2\pi i w A}}{1+i A}\right) - \frac{1}{(2\pi w)^2} \left( 1+ \frac{e^{-2\pi i w A}}{ (A-i)^2}\right). $$

  • 0
    hm, now I am confused. I have checked,Mathematica gives me $e^{2\pi w}$. I would prefer $e^{-2\pi w}$ too2011-05-05
  • 0
    @Katja: $e^{2\pi w}$ is wrong. My post shows how the integral should behave for large $w$.2011-05-05
  • 0
    @Katja: the two terms which I have given in my post approximate the integral reasonable well for almost all values of $w\gtrsim 1$ and $A$.2011-05-05
  • 0
    Tsk, I already upvoted this (and thus cannot upvote it again), but I agree; if $w$ is big, this is a better choice for the numerics.2011-05-05