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

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

25.8 APPROXIMATIONS TO INTEGRALSFinally, we should mention that the lines in the z-plane on which the exponentsin the WKB solutions are purely imaginary, <strong>and</strong> the two solutions haveequal amplitudes, are usually called the anti-Stokes lines. For the general Bessel’sequation they are the real positive <strong>and</strong> real negative axes.25.8 Approximations to integralsIn this section we will investigate a method of finding approximations to thevalues or <strong>for</strong>ms of certain types of infinite integrals. The class of integrals tobe considered is that containing integr<strong>and</strong>s that are, or can be, represented byexponential functions of the general <strong>for</strong>m g(z) exp[ f(z) ]. The exponents f(z) maybe complex, <strong>and</strong> so integrals of sinusoids can be h<strong>and</strong>led as well as those withmore obvious exponential properties. We will be using the analyticity propertiesof the functions of a complex variable to move the integration path to a partof the complex plane where a general integr<strong>and</strong> can be approximated well by ast<strong>and</strong>ard <strong>for</strong>m; the st<strong>and</strong>ard <strong>for</strong>m is then integrated explicitly.The particular st<strong>and</strong>ard <strong>for</strong>m to be employed is that of a Gausssian function ofa real variable, <strong>for</strong> which the integral between infinite limits is well known. This<strong>for</strong>m will be generated by expressing f(z) as a Taylor series expansion about apoint z 0 , at which the linear term in the expansion vanishes, i.e. where f ′ (z) =0.Then, apart from a constant multiplier, the exponential function will behave likeexp[ 1 2 f′′ (z 0 )(z − z 0 ) 2 ] <strong>and</strong>, by choosing an appropriate direction <strong>for</strong> the contourto take as it passes through the point, this can be made into a normal Gaussianfunction of a real variable <strong>and</strong> its integral may then be found.25.8.1 Level lines <strong>and</strong> saddle pointsBe<strong>for</strong>e we can discuss the method outlined above in more detail, a number ofobservations about functions of a complex variable <strong>and</strong>, in particular, about theproperties of the exponential function need to be made. For a general analyticfunction,f(z) =φ(x, y)+iψ(x, y), (25.58)of the complex variable z = x + iy, we recall that, not only do both φ <strong>and</strong> ψsatisfy Laplace’s equation, but ∇φ <strong>and</strong> ∇ψ are orthogonal. This means that thelines on which one of φ <strong>and</strong> ψ is constant are exactly the lines on which the otheris changing most rapidly.Let us apply these observations to the functionh(z) ≡ exp[ f(z) ] = exp(φ) exp(iψ), (25.59)recalling that the functions φ <strong>and</strong> ψ are themselves real. The magnitude of h(z),given by exp(φ), is constant on the lines of constant φ, which are known as the905

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

Saved successfully!

Ooh no, something went wrong!