saltydog
Science Advisor
Homework Helper
- 1,590
- 3
I wish to finally evaluate:
[tex] y(t)=\frac{e^{-\pi/t}}{2 i t^{3/2}}\mathop\lim\limits_{R\to\infty}\left\{\text{Erfi}\left(v-\sqrt{\pi/t}\right)_{\sqrt{Rt}e^{-\pi i/4}}^{\sqrt{Rt}e^{\pi i/4}}\right\}\tag{1}[/tex]
First note:
[tex] \text{Erfi}(z)=\frac{\text{Erf(iz)}}{i}=\frac{1}{i}\frac{2}{\sqrt{\pi}}\int_0^{iz} e^{-w^2}dw[/tex]
Using this relation, I can re-write the above limit as:
[tex] \frac{2}{i\sqrt{\pi}}\left\{\mathop\lim\limits_{R\to\infty}\int_0^{i(\sqrt{Rt}e^{\pi i/4})} e^{-w^2}dw-\mathop\lim\limits_{R\to\infty}\int_0^{i(\sqrt{Rt}e^{-\pi i/4})} e^{-w^2}dw\right\}[/tex]
Using Euler's relation, I convert the upper limits to real and imaginary parts and multiply by i:
[tex] \begin{align*}<br /> i\sqrt{Rt}e^{\pi i/4}&=i\left[\left(\sqrt{Rt/2}-\sqrt{\pi/2}\right)+i\sqrt{Rt/2}\right] \\<br /> &=-\sqrt{Rt/2}+i\left(\sqrt{Rt/2}-\sqrt{\pi/t}\right)<br /> \end{align}[/tex]
Now, as R goes to infinity, this quantity goes to [itex]-\infty+i\infty[/itex].
And the lower limit:
[tex] \begin{align*}<br /> i\sqrt{Rt}e^{-\pi i/4}&=i\left[\sqrt{Rt/2}+\left(\sqrt{Rt/2}-\sqrt{\pi/2}\right)\right] \\<br /> &=\sqrt{Rt/2}+i\left(\sqrt{Rt/2}-\sqrt{\pi/2}\right)<br /> \end{align}[/tex]
as R goes to infinity, this limit goes to [itex]\infty+i\infty[/itex].
Thus the integrals can be rewritten as:
[tex] \frac{1}{i}\frac{2}{\sqrt{\pi}}\int_0^{-\infty+i\infty} e^{-w^2}dw-<br /> \frac{1}{i}\frac{2}{\sqrt{\pi}}\int_0^{\infty+i\infty} e^{-w^2}dw[/tex]
I choose to evaluate these line integrals along a parametric path using the following relation from Complex Analysis:
[tex] \int_C f(z)dz=\int_c udx-vdy+i\int_c vdx+udy[/tex]
Thus the first step is to express the integrand in terms of a real and imaginary part. Relying on Euler's formula:
[tex] \begin{align*}<br /> e^{-w^2}&=e^{-(x+yi)^2} \\<br /> &=e^{-(x^2-y^2+2xyi)} \\<br /> &=e^{-(x^2-y^2)}Cos(2xy)+ie^{-(x^2-y^2)}(-Sin(2xy))\\<br /> &=u+iv<br /> \end{align}[/tex]
For the first integral, the path goes from the origin to the point [itex]-\infty+i\infty[/itex]
Thus if I let:
[tex]x=-t,\quad dx=-dt[/tex]
[tex]y=t,\quad dy=dt[/tex]
Making this substitution, note that the exponent term is 1 and we have for the first integral (with t going from 0 to infinity):
[tex] \begin{align*}<br /> \int_0^{-\infty+i\infty} e^{-w^2}dw&=\left(-\int_0^\infty Cos(-2t^2)dt+\int_0^\infty Sin(-2t^2)dt\right) \\<br /> &+i\left(\int_0^\infty Cos(-2t^2)dt+\int_0^\infty Sin(-2t^2)dt\right) \\<br /> &=\left(-\frac{\sqrt{\pi}}{4}-\frac{\sqrt{\pi}}{4}\right)+i\left(\frac{\sqrt{\pi}}{4}-\frac{\sqrt{\pi}}{4}\right) \\<br /> &=-\frac{\sqrt{\pi}}{2}<br /> \end{align}[/tex]
The same analysis can be made (with the path going from 0 to [itex]\infty+i\infty[/itex] for the second integral with the results that the line integral in that case is:
[tex] \frac{\sqrt{\pi}}{2}[/tex]
Substituting these values into (1), I obtain:
[tex] \frac{1}{i}\frac{2}{\sqrt{\pi}}\int_0^{-\infty+i\infty} e^{-w^2}dw=\frac{1}{i}\frac{2}{\sqrt{\pi}}(-\frac{\sqrt{\pi}}{2})=-\frac{1}{i}\frac{i}{i}=i[/tex]
and:
[tex] \frac{1}{i}\frac{2}{\sqrt{\pi}}\int_0^{\infty+i\infty} e^{-w^2}dw=\frac{1}{i}\frac{2}{\sqrt{\pi}}(\frac{\sqrt{\pi}}{2})=\frac{1}{i}\frac{i}{i}=-i[/tex]
So that the limit is [itex]2i[/itex].
Substituting this into the last expression of the inverse transform, we obtain:
[tex] y(x)=\frac{e^{-\pi/t}}{2i t^{3/2}}(2i)=\frac{e^{-\pi/t}}{t^{3/2}}[/tex]
[tex] y(t)=\frac{e^{-\pi/t}}{2 i t^{3/2}}\mathop\lim\limits_{R\to\infty}\left\{\text{Erfi}\left(v-\sqrt{\pi/t}\right)_{\sqrt{Rt}e^{-\pi i/4}}^{\sqrt{Rt}e^{\pi i/4}}\right\}\tag{1}[/tex]
First note:
[tex] \text{Erfi}(z)=\frac{\text{Erf(iz)}}{i}=\frac{1}{i}\frac{2}{\sqrt{\pi}}\int_0^{iz} e^{-w^2}dw[/tex]
Using this relation, I can re-write the above limit as:
[tex] \frac{2}{i\sqrt{\pi}}\left\{\mathop\lim\limits_{R\to\infty}\int_0^{i(\sqrt{Rt}e^{\pi i/4})} e^{-w^2}dw-\mathop\lim\limits_{R\to\infty}\int_0^{i(\sqrt{Rt}e^{-\pi i/4})} e^{-w^2}dw\right\}[/tex]
Using Euler's relation, I convert the upper limits to real and imaginary parts and multiply by i:
[tex] \begin{align*}<br /> i\sqrt{Rt}e^{\pi i/4}&=i\left[\left(\sqrt{Rt/2}-\sqrt{\pi/2}\right)+i\sqrt{Rt/2}\right] \\<br /> &=-\sqrt{Rt/2}+i\left(\sqrt{Rt/2}-\sqrt{\pi/t}\right)<br /> \end{align}[/tex]
Now, as R goes to infinity, this quantity goes to [itex]-\infty+i\infty[/itex].
And the lower limit:
[tex] \begin{align*}<br /> i\sqrt{Rt}e^{-\pi i/4}&=i\left[\sqrt{Rt/2}+\left(\sqrt{Rt/2}-\sqrt{\pi/2}\right)\right] \\<br /> &=\sqrt{Rt/2}+i\left(\sqrt{Rt/2}-\sqrt{\pi/2}\right)<br /> \end{align}[/tex]
as R goes to infinity, this limit goes to [itex]\infty+i\infty[/itex].
Thus the integrals can be rewritten as:
[tex] \frac{1}{i}\frac{2}{\sqrt{\pi}}\int_0^{-\infty+i\infty} e^{-w^2}dw-<br /> \frac{1}{i}\frac{2}{\sqrt{\pi}}\int_0^{\infty+i\infty} e^{-w^2}dw[/tex]
I choose to evaluate these line integrals along a parametric path using the following relation from Complex Analysis:
[tex] \int_C f(z)dz=\int_c udx-vdy+i\int_c vdx+udy[/tex]
Thus the first step is to express the integrand in terms of a real and imaginary part. Relying on Euler's formula:
[tex] \begin{align*}<br /> e^{-w^2}&=e^{-(x+yi)^2} \\<br /> &=e^{-(x^2-y^2+2xyi)} \\<br /> &=e^{-(x^2-y^2)}Cos(2xy)+ie^{-(x^2-y^2)}(-Sin(2xy))\\<br /> &=u+iv<br /> \end{align}[/tex]
For the first integral, the path goes from the origin to the point [itex]-\infty+i\infty[/itex]
Thus if I let:
[tex]x=-t,\quad dx=-dt[/tex]
[tex]y=t,\quad dy=dt[/tex]
Making this substitution, note that the exponent term is 1 and we have for the first integral (with t going from 0 to infinity):
[tex] \begin{align*}<br /> \int_0^{-\infty+i\infty} e^{-w^2}dw&=\left(-\int_0^\infty Cos(-2t^2)dt+\int_0^\infty Sin(-2t^2)dt\right) \\<br /> &+i\left(\int_0^\infty Cos(-2t^2)dt+\int_0^\infty Sin(-2t^2)dt\right) \\<br /> &=\left(-\frac{\sqrt{\pi}}{4}-\frac{\sqrt{\pi}}{4}\right)+i\left(\frac{\sqrt{\pi}}{4}-\frac{\sqrt{\pi}}{4}\right) \\<br /> &=-\frac{\sqrt{\pi}}{2}<br /> \end{align}[/tex]
The same analysis can be made (with the path going from 0 to [itex]\infty+i\infty[/itex] for the second integral with the results that the line integral in that case is:
[tex] \frac{\sqrt{\pi}}{2}[/tex]
Substituting these values into (1), I obtain:
[tex] \frac{1}{i}\frac{2}{\sqrt{\pi}}\int_0^{-\infty+i\infty} e^{-w^2}dw=\frac{1}{i}\frac{2}{\sqrt{\pi}}(-\frac{\sqrt{\pi}}{2})=-\frac{1}{i}\frac{i}{i}=i[/tex]
and:
[tex] \frac{1}{i}\frac{2}{\sqrt{\pi}}\int_0^{\infty+i\infty} e^{-w^2}dw=\frac{1}{i}\frac{2}{\sqrt{\pi}}(\frac{\sqrt{\pi}}{2})=\frac{1}{i}\frac{i}{i}=-i[/tex]
So that the limit is [itex]2i[/itex].
Substituting this into the last expression of the inverse transform, we obtain:
[tex] y(x)=\frac{e^{-\pi/t}}{2i t^{3/2}}(2i)=\frac{e^{-\pi/t}}{t^{3/2}}[/tex]
Last edited: