Exercise 7.   [Graphics:Images/IntegralsBranchPointsModHome_gr_513.gif].  

Solution 7.

See text and/or instructor's solution manual.

Answer.   [Graphics:../Images/IntegralsBranchPointsModHome_gr_514.gif].  

Solution.  We use the function   [Graphics:../Images/IntegralsBranchPointsModHome_gr_515.gif],   where the branch of the logarithm is    

            [Graphics:../Images/IntegralsBranchPointsModHome_gr_516.gif],  

where  [Graphics:../Images/IntegralsBranchPointsModHome_gr_517.gif]  and  [Graphics:../Images/IntegralsBranchPointsModHome_gr_518.gif],    (see Equation (5-20) in Section 5.2),  (notice that the branch cut is the negative y-axis).  

      The path  C  of integration will consist of the segments  [Graphics:../Images/IntegralsBranchPointsModHome_gr_519.gif]  and  [Graphics:../Images/IntegralsBranchPointsModHome_gr_520.gif] of the x-axis

together with the upper semicircles  [Graphics:../Images/IntegralsBranchPointsModHome_gr_521.gif]  and  [Graphics:../Images/IntegralsBranchPointsModHome_gr_522.gif],  for  [Graphics:../Images/IntegralsBranchPointsModHome_gr_523.gif],  as shown in the figure below.

                                         [Graphics:../Images/IntegralsBranchPointsModHome_gr_524.gif]

                    The indented contour  [Graphics:../Images/IntegralsBranchPointsModHome_gr_525.gif]  for the integrand  [Graphics:../Images/IntegralsBranchPointsModHome_gr_526.gif]

                    consists of the large semi-circle [Graphics:../Images/IntegralsBranchPointsModHome_gr_527.gif], and a small semi-circle [Graphics:../Images/IntegralsBranchPointsModHome_gr_528.gif], and segments  [Graphics:../Images/IntegralsBranchPointsModHome_gr_529.gif]  along the real axis.  

                    The point  [Graphics:../Images/IntegralsBranchPointsModHome_gr_530.gif]  that lies inside the contour  [Graphics:../Images/IntegralsBranchPointsModHome_gr_531.gif].  

      Factor the denominator and get

                    [Graphics:../Images/IntegralsBranchPointsModHome_gr_532.gif].  

The zeros of  [Graphics:../Images/IntegralsBranchPointsModHome_gr_533.gif] are  [Graphics:../Images/IntegralsBranchPointsModHome_gr_534.gif].

It follows that  [Graphics:../Images/IntegralsBranchPointsModHome_gr_535.gif]  has simple poles at  [Graphics:../Images/IntegralsBranchPointsModHome_gr_536.gif].

      Now evaluate the contour integral of  [Graphics:../Images/IntegralsBranchPointsModHome_gr_537.gif]  over  [Graphics:../Images/IntegralsBranchPointsModHome_gr_538.gif].  

We chose the branch  [Graphics:../Images/IntegralsBranchPointsModHome_gr_539.gif]  because it is analytic on C and its interior.  

Applying the Cauchy Residue Theorem (see Theorem 8.1 in Section 8.1), we obtain:  

            [Graphics:../Images/IntegralsBranchPointsModHome_gr_540.gif]  

Here the denominator of  [Graphics:../Images/IntegralsBranchPointsModHome_gr_541.gif]  has a factor of the form  [Graphics:../Images/IntegralsBranchPointsModHome_gr_542.gif],  and  [Graphics:../Images/IntegralsBranchPointsModHome_gr_543.gif].  

            [Graphics:../Images/IntegralsBranchPointsModHome_gr_544.gif]  

      We can write the contour integral in the following form:

                    [Graphics:../Images/IntegralsBranchPointsModHome_gr_545.gif],  
                    
                    or
                    
                    [Graphics:../Images/IntegralsBranchPointsModHome_gr_546.gif].  

      The ML-Inequality and a proof similar to the one used to establish equation (8-11) can be used to prove that:

                    [Graphics:../Images/IntegralsBranchPointsModHome_gr_547.gif]     and     [Graphics:../Images/IntegralsBranchPointsModHome_gr_548.gif]     (see the complete details below).

Taking the limits as  [Graphics:../Images/IntegralsBranchPointsModHome_gr_549.gif]  we obtain

                    [Graphics:../Images/IntegralsBranchPointsModHome_gr_550.gif],

                    then
                    
                    [Graphics:../Images/IntegralsBranchPointsModHome_gr_551.gif],

                    then

                    [Graphics:../Images/IntegralsBranchPointsModHome_gr_552.gif],

                    then

                    [Graphics:../Images/IntegralsBranchPointsModHome_gr_553.gif],
                    
                    then

                    [Graphics:../Images/IntegralsBranchPointsModHome_gr_554.gif].  

Equating the real parts in this equation gives  

                    [Graphics:../Images/IntegralsBranchPointsModHome_gr_555.gif].  

From ordinary calculus we know that

                    [Graphics:../Images/IntegralsBranchPointsModHome_gr_556.gif].
                    
Making this substitution, we get   

                    [Graphics:../Images/IntegralsBranchPointsModHome_gr_557.gif],

which is easy to solve for the desired result

                    [Graphics:../Images/IntegralsBranchPointsModHome_gr_558.gif]
                    

We are done.   

      For curiosity sake, the imaginary part is

                    [Graphics:../Images/IntegralsBranchPointsModHome_gr_559.gif],
                    
                    or
                    
                    [Graphics:../Images/IntegralsBranchPointsModHome_gr_560.gif].  

Aside.  We can let Mathematica double check our work.

[Graphics:../Images/IntegralsBranchPointsModHome_gr_561.gif]

[Graphics:../Images/IntegralsBranchPointsModHome_gr_562.gif]


[Graphics:../Images/IntegralsBranchPointsModHome_gr_563.gif]

[Graphics:../Images/IntegralsBranchPointsModHome_gr_564.gif]


[Graphics:../Images/IntegralsBranchPointsModHome_gr_565.gif]

[Graphics:../Images/IntegralsBranchPointsModHome_gr_566.gif]

Maple can check our work too!

     > V1 := evalc( residue( (log(z))^2/(z^2+1), z=I ));

               [Graphics:../Images/IntegralsBranchPointsModHome_gr_567.gif]

     > V := evalc( 2*Pi*I*V1 );

               [Graphics:../Images/IntegralsBranchPointsModHome_gr_568.gif]

     > ( Pi^3/2 - Pi^3/4 )/2;

               [Graphics:../Images/IntegralsBranchPointsModHome_gr_569.gif]

We are really done.   

                     [Graphics:../Images/IntegralsBranchPointsModHome_gr_570.gif]

                    A portion of the "area under the curve"  [Graphics:../Images/IntegralsBranchPointsModHome_gr_571.gif].  

We are really really done.   

The complete details.

     For completeness, we now give the details for establishing  [Graphics:../Images/IntegralsBranchPointsModHome_gr_572.gif]  and  [Graphics:../Images/IntegralsBranchPointsModHome_gr_573.gif].  

Use the ML-Inequality (Theorem 6.3 in Section 6.2) and a proof similar to the one used to establish equation (8-11) of Theorem 8.3 to prove that:

                    [Graphics:../Images/IntegralsBranchPointsModHome_gr_574.gif].  

For this example, for  [Graphics:../Images/IntegralsBranchPointsModHome_gr_575.gif],  we have the relations  

                    [Graphics:../Images/IntegralsBranchPointsModHome_gr_576.gif]

Now taking limits, we get

                    [Graphics:../Images/IntegralsBranchPointsModHome_gr_577.gif].

Therefore, it follows that  [Graphics:../Images/IntegralsBranchPointsModHome_gr_578.gif].  

A similar proof shows that:

                    [Graphics:../Images/IntegralsBranchPointsModHome_gr_579.gif].  

For this example, for  [Graphics:../Images/IntegralsBranchPointsModHome_gr_580.gif],  we have the relations  

                    [Graphics:../Images/IntegralsBranchPointsModHome_gr_581.gif]  

Now taking limits, we get

                    [Graphics:../Images/IntegralsBranchPointsModHome_gr_582.gif].

Therefore, it follows that  [Graphics:../Images/IntegralsBranchPointsModHome_gr_583.gif].  

We are done with the details.

Aside.  Both [Graphics:../Images/IntegralsBranchPointsModHome_gr_584.gif] and [Graphics:../Images/IntegralsBranchPointsModHome_gr_585.gif] are capable of finding the definite integral.

[Graphics:../Images/IntegralsBranchPointsModHome_gr_586.gif]

[Graphics:../Images/IntegralsBranchPointsModHome_gr_587.gif]


     > int( (log(x))^2/(x^2+1), x=0..infinity );

               [Graphics:../Images/IntegralsBranchPointsModHome_gr_588.gif]

We are really really really done.

The above computation illustrates how the Residue Calculus can be used to find Contour Integrals with Integrands with Branch Points.  

Aside.  For comparison, we show how to find the value of the integral using the indefinite integral.

                               [Graphics:../Images/IntegralsBranchPointsModHome_gr_589.gif]

                    The indefinite integral  [Graphics:../Images/IntegralsBranchPointsModHome_gr_590.gif].   From the graph we see that  

                    [Graphics:../Images/IntegralsBranchPointsModHome_gr_591.gif].

Caveat.  If you do not trust the graphical explanation then you must trust the Residue Calculus solution.  

Or you could compute the Cauchy Principal value, but this will require some very difficult limits.

                    The Cauchy Principal Value (P.V.) of the integral is  [Graphics:../Images/IntegralsBranchPointsModHome_gr_592.gif]  and it can be computed as follows:

                    [Graphics:../Images/IntegralsBranchPointsModHome_gr_593.gif]  

      The indefinite integral will involve the special function "PolyLog."  

      We do not want to develop this special functions at this time and leave it for the reader to investigate.  It is interesting to observe that the current software includes this function.  For this type of problem the Residue Calculus can get the answer without resorting to special functions.

[Graphics:../Images/IntegralsBranchPointsModHome_gr_594.gif]

[Graphics:../Images/IntegralsBranchPointsModHome_gr_595.gif]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 



This solution is complements of the authors.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) 2008 John H. Mathews, Russell W. Howell