13.07.2015 Views

Mathematical Methods for Physics and Engineering - Matematica.NET

Mathematical Methods for Physics and Engineering - Matematica.NET

Mathematical Methods for Physics and Engineering - Matematica.NET

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

APPLICATIONS OF COMPLEX VARIABLESFinally, putting the various values into the <strong>for</strong>mula yields( ) 1/2 2πF(x) ∼ + g(i)exp[f(i)]exp[ 1 i(π − α)]2A( ) 1/2 (2π x 1/2=+2x 3/2 2π exp − 2 )3 x3/2 exp[ 1 i(π − π)]2(1=2 √ πx exp − 2 )1/4 3 x3/2 .This is the leading term in the asymptotic expansion of F(x),which,asshowninequation(25.39), is a particular contour integral solution of Stokes’ equation. The fact that it tendsto zero in a monotonic way as x → +∞ allows it to be identified with the Airy function,Ai(x).We may ask why the saddle point at t = −i was not used. The answer to this is asfollows. Of course, any path that starts <strong>and</strong> ends in the right sectors will suffice, butif another saddle point exists close to the one used, then the Taylor expansion actuallyemployed is likely to be less effective than if there were no other saddle points or if therewere only distant ones.An investigation of the same <strong>for</strong>m as that used at t =+i shows that the saddle at t = −iis higher by a factor of exp( 4 3 x3/2 ) <strong>and</strong> that its l.s.d. is orientated parallel to the imaginaryt-axis. Thus a path that went through it would need to go via a region of largish negativeimaginary t, over the saddle at t = −i, <strong>and</strong> then, when it reached the col at t =+i, bendsharply <strong>and</strong> follow part of the same l.s.d. as considered earlier. Thus the contribution fromthe t = −i saddle would be incomplete <strong>and</strong> roughly half of that from the t =+i saddlewould still have to be included. The more serious error would come from the first of these,as, clearly, the part of the path that lies in the plane Re t = 0 is not symmetric <strong>and</strong> is farfrom Guassian-like on the side nearer the origin. The Gaussian-path approximation usedwill there<strong>for</strong>e not be a good one, <strong>and</strong>, what is more, the resulting error will be magnified bya factor exp( 4 3 x3/2 ) compared with the best estimate. So, both on the grounds of simplicity<strong>and</strong> because the effect of the other (neglected) saddle point is likely to be less severe, wechoose to use the one at t =+i. ◭25.8.3 Stationary phase methodIn the previous subsection we showed how to use the saddle points of anexponential function of a complex variable to evaluate approximately a contourintegral of that function. This was done by following the lines of steepest descentthat passed through the saddle point; these are lines on which the phase of theexponential is constant but its amplitude is varying at the maximum possiblerate <strong>for</strong> that function. We now introduce an alternative method, one that entirelyreverses the roles of amplitude <strong>and</strong> phase. To see how such an alternative approachmight work, it is useful to study how the integral of an exponential function ofa complex variablecan be represented as the sum of infinitesimal vectors in thecomplex plane.We start by studying the familiar integralI 0 =∫ ∞−∞exp(−z 2 ) dz, (25.67)912

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!