Exercise 5.   [Graphics:Images/IntegralsIndentedContourModHome_gr_168.gif]   and   [Graphics:Images/IntegralsIndentedContourModHome_gr_169.gif].

Solution 5.

See text and/or instructor's solution manual.

Answer.   [Graphics:../Images/IntegralsIndentedContourModHome_gr_170.gif]   and   [Graphics:../Images/IntegralsIndentedContourModHome_gr_171.gif].  

Solution.  The complex integrand is  [Graphics:../Images/IntegralsIndentedContourModHome_gr_172.gif].  

The denominator is   [Graphics:../Images/IntegralsIndentedContourModHome_gr_173.gif][Graphics:../Images/IntegralsIndentedContourModHome_gr_174.gif].

The complex zeros are  [Graphics:../Images/IntegralsIndentedContourModHome_gr_175.gif],  and the real zeros are  [Graphics:../Images/IntegralsIndentedContourModHome_gr_176.gif].

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

                    The indented contour C that consists of the semi-circle [Graphics:../Images/IntegralsIndentedContourModHome_gr_178.gif], and segments along the real axis and small semi-circles

                    above the points  [Graphics:../Images/IntegralsIndentedContourModHome_gr_179.gif]  that lie on the real axis.  The point  [Graphics:../Images/IntegralsIndentedContourModHome_gr_180.gif]  lies inside  C.  

Using Theorem 8.1 (Cauchy's Residue Theorem), and Theorem 8.5 (Contour Integration for Indented Contour Integrals), the value of the integral is computed   

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

We are done.   

Aside.  We can let Mathematica double check our work.

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

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


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

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


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

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


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

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

Maple can check our work too!

     > V1 := residue( z^2/(z^4-1), z=I );

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

     > V2 := residue( z^2/(z^4-1), z=-1 );

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

     > V3 := residue( z^2/(z^4-1), z=1 );

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

     > 2*Pi*I*V1 + Pi*I*(V2+V3);

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

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

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

We are really done.   

Aside.  Both [Graphics:../Images/IntegralsIndentedContourModHome_gr_196.gif] and [Graphics:../Images/IntegralsIndentedContourModHome_gr_197.gif] are capable of finding the definite integral, but the "Principal Value" option must be used.

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

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


     > f := proc(x) x^2/(x^4-1) end proc;

                
f := proc(x) x^2/(x^4 - 1) end proc

     > simplify(convert(int(f(x), x=-infinity..infinity, 'CauchyPrincipalValue'),exp));

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

We are really really done.

The above computation illustrates how the Residue Calculus can be used to find Indented Contour Integrals.

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

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

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

                    [Graphics:../Images/IntegralsIndentedContourModHome_gr_203.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 difficult limits.

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

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

The details are tedious and we let Mathematica do the work.

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

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


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

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

Alternately, we have the computation:

[Graphics:../Images/IntegralsIndentedContourModHome_gr_210.gif]
[Graphics:../Images/IntegralsIndentedContourModHome_gr_211.gif]

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


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

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 



This solution is complements of the authors.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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