Showing posts with label singular points. Show all posts
Showing posts with label singular points. Show all posts

Monday, January 27, 2014

More Complex Integrals


We continue now looking at more examples of complex integration and application of the residue theorem.

Example 1:

-¥  ¥  (1 + x2 )  dx / 1 + x4 


Find the residues at exp(pi/4)  and exp(3pi/4)  and hence integrate the preceding:




Solution:



 Begin by identifying the singular points where:

f(z) = 1 + x2  / 1 + x4 

Now:  We note that  1 + x4     =   0 when   x = (-1) ¼   

Or the fourth root of (-1). The modulus of (-1) = 1 and the argument is pTherefore, by de Moivre’s theorem:

(-1  + 0i) ¼     =   4Ö(-1) [ cos (p + 2kp/ 4) + i sin (p + 2kp/ 4)]
k=   0, 1, 2, 3…….

At k= 0:     (-1) ¼      =  cis (p/4)   =   exp (pi/4)

At k=1:   (-1) ¼      =  cis (3p/4)    =      exp (3pi/4) 


(We’re  not concerned with k=2, 3……Why?)

We therefore take the semi-circular contour as shown in the graphic. So that:

òCR   f(z) dz  +   -R  R  f(x) dx  =  2 pi  (sum of residues inside C)

As R increases without bound we have:
-¥  ¥  (1 + x2 )  dx / 1 + x4     =


2 pi [Res f(z) at exp (pi/4) and exp (3pi/4) ]

For the case f(z) = p(z)/ q(z)   Then:

Res f(z) z = a      =  lim z ® a   [p(z)/ q’(z)]

Then for z = exp (pi/4):

Res f(z) z = exp pi/4       =   lim z ® exp pi/4      1 + z2  / 4 z3     =

1 +  exp(pi/2) / 4 exp(3pi/4) =  1 + cis (p/2) / 4 cis (3p/4)

Res f(z) z = exp pi/4       =    1 + i /  4(- Ö2/ 2   + iÖ2/ 2   )

=  1/ 2Ö2 (i+ 1/ i- 1)

Rationalize the denominator to get:  2i/ -2  =  - Ö2i/ 4
Then for z = exp (3pi/4):

Res f(z) z = exp 3pi/4       =   lim z ® exp 3pi/4      1 + z2  / 4 z3     =

1 +  exp(3pi/4) / 4 exp(9pi/4) = 

¼ [exp(-3pi/4)  +  exp(- 9pi/4)] = 


¼ [ cis (-3pi/4)  + cis (- 9pi/4)]

Therefore:
-¥  ¥  (1 + x2 )  dx / 1 + x4     =


2 pi [Res f(z) at exp (pi/4) and exp (3pi/4) }

=  2 pi {- Ö2i/ 4   +  ¼ [ cis (-3pi/4)  + cis (- 9pi/4)]}

Problem for Serious Math Mavens:

Obtain the integral for:

  -¥  ¥    x  dx / (x2   - 2x + 2)
 (Hint: Remember  the residue formula (i.e. for a given pole of order m) and your complex numbers!)

Thursday, December 26, 2013

More Laurent Series Examples


No, we’re not yet finished with these exotic mathematical denizens! We want to give yet more examples, and especially showing different ways to obtain the Laurent equivalent of different functions.

 
Ex. (1):   exp(2z)/ (z - 1)3  about  z 0  =   1

 
Let u = (z – 1)  so that z = (u -1)


Then:

 exp ( 2z) / (z - 1)3     =  exp(2 + 2u) /  u 3  =   [( e) 2  /  u 3   ] /  e 2u 

 
=   e 2  /  u 3      [ 1  + 2u   + (2u) 2 / 2!   +   (2u) 3 / 3!   + 

 (2u) 4 / 4! +     …….   

=  e 2  /  (z - 1)3     +   2 e 2  /  (z - 1)2      +   2 e 2  /  (z – 1)  +  

4 e 2  / 3  +  2 e 2  / 3 (z – 1)  +   …….

 
Note that z is a pole of order 3 or a “triple pole” (Why?)    Note also the series converges for all values of z  ¹  1



Ex. (2):   f(z) =  z – sin z /  z 3

We’ve already seen the series for sin z:

 
sin(z) = z -  z3 / 3!  +  z5 / 5!  …..  -  (-1) n- 1 z 2n -1/  (2n -1)!  +
 
Then:
z – sin z /  z 3     =  1/  z 3    [z -  (z -  z3 / 3!  +  z5 / 5!  …..  ) ] 
 
 =  1/  z 3    [z3 / 3!   -  z5 / 5!    +   z 7 / 7!  -    ……….
 
 
=   1/3!  -  z 2 / 5!     +   z 4 / 7!  -    z 6 / 9!   + ……..
 

Ex. (3)
 Find a Laurent series for f(z) =  sin (4z) / z 4
 
Using the sine series we can write:
 
sin (4z) / z 4     =  1/   z 4   [4z -    4z3 / 3!  +  4 z5 / 5!  +  ……..]
 
=     4 /   z 3    -    43 / 3!  z  +   45  z  / 5!       +   …..
 
= 4 z  -3    -  43  z -1 /  3!   +   45  z  / 5!        +   …..
 
 
 
Which on inspection is found to have the form:
 
å¥ n = 0    (-1) n   4 2n+ 1  z 2n -3  /   (2n +  1)!
Check term number 3 (n = 2):
 
4 2n+ 1     = 4 2(2)+ 1      =      45
 
z 2n -3  =  z 2(2) - 3      = z 1=  z
 
(2n +  1)!   =  (2(2) +  1)!  = 5!
 
Then:  we get  45  z  / 5!        
 
A removable singularity occurs at z = 0 so we expect the region of convergence to be:
 
÷ z  ÷      >   0   
 

Problems for Math Mavens:
 
Let f(z) = exp (-1/z 2) /  z5
 
a) Show the Laurent series can be written: 
 
å¥ n = 0    (-1) n   /   n !  z 2n + 5   

b)Specify where the singularity would occur and the precise region of convergence
 
 
Challenge Problem: For über Math Mavens only –

 
Let f(z) = 7 z 2  + 9 z  - 18 /  z 3 -  9z

 
Find Laurent series for the convergence regions:
 
a)  0  < ÷ z ÷   <    3  and
b)    ÷ z  ÷      >   3  
 
(Hint: Approach the function f(z) by resort to partial fraction decomposition)



 
 

Sunday, December 22, 2013

More on Laurent Series


Note to readers: I reiterate once more that this set of blog posts isn’t intended to be comprehensive nor should they be seen as replacing the hard work of diligently working problems on your own using available text books. They are merely intended to spur on curiosity and perhaps offer a few different perspectives from what one might have seen already.


Included here are the coefficients in typical Laurent series which are generally obtained by other means (i.e. than appealing directly to integral representations). Some examples I have included below might suffice to expose this.

 

Ex. (1):  You have the Maclaurin series:

 

exp (z) = å¥ n = 0   (z) n  / n!  = 1 + z + z2 / 2!  + z3 / 3!  +…..

 

 
Now, replace z by 1/z in the expansion:
 
 
exp( 1/z) = å¥ n = 0   1  / n!   z n  = 1 + 1/ 1! z   + 1/ 2!  z2 

 +  1/ 3! z3   +…..

 
for: 0   <   ÷ z ÷    <     ¥ 

 
Which then becomes a Laurent series expansion.
 
To establish this on inspection note that no positive powers of z appear, only negative, i.e. then 1/z is 
z - 1   

 So, in effect we can say that the coefficients of the positive powers are zero.  It’s also important to note here that the coefficient of : 
 
1/ 1! z     = 1/ z
 
is unity, so according to Laurent’s theorem we designate that the coefficient:

 
c n  =  1/2 pi  òC  exp( 1/z)  dz

 
where C is any positively oriented simple closed contour around the origin.  Now, since from our prior  Dec. 12, Introducing Laurent Series) examination c n  =  1 then:
 
=1/2 pi  òC  exp( 1/z)  dz   and:  òC  exp( 1/z)  dz    =  2 pi 

Ex. (2):  Consider the function:

f(z) = 1/ (z – i) 2
 
This is already in the form of a Laurent series , where z 0  =   i.  That is,
 
å¥ n = -¥   c n (z – i) n     (For:  0   <   ÷ z - i ÷    <     ¥  )

Where in this case,  c - 2  =  1  and all other coefficients are zero.  Then from the previous theorems, relations we’ve explored (blog post of December 12, ‘Introduction to Laurent Series’, sub-header: ‘More intricacies and singularities’):
 
c n  =  1/2 pi  òC   dz/ ( z – i) n +3      n = 0, +1, +2, +3 …..
 
 
where C denotes any positively oriented circle ÷ z - i ÷   = R  about the point: z 0  =   i
 
 
Then it follows from this:
 
 
a)    òC   dz/ ( z – i) n +3      =  0   (when n ¹   2)

 
b) òC   dz/ ( z – i) n +3      =  1    (when n=    2)
 
Problems for the Math Maven:
 
1)     The function: f(z) =  -1/ (z – 1) (z – 2)
 
a) Rewrite it using the partial fraction form

 
b) Identify the two singular points.
 
 
c) If  the function is analytic in the domains:

 
÷ z ÷    <   1;   1     < ÷ z ÷     <      2   and:   2    ÷ z ÷    <     ¥,  
 
 
Draw a sketch showing the different domains.
 
 
2) Find the Laurent series that represents the function:
 
 
f(z) =  z2  sin (1 / z2 )
 
In the domain:  : 0   ÷ z ÷   <     ¥ 
 
 
Hint:  Recall the series for sin(z) = z -  z3 / 3!  +  z5 / 5!  …..-  
   (-1) n- 1 z 2n -1/  (2n -1)!  +   ……   (For ÷ z ÷   <    ¥  )