Complex Variables. Chapter 18 Integration in the Complex Plane. March 12, 2013 Lecturer: Shih-Yuan Chen

Similar documents
One Dimension Again. Chapter Fourteen

Chapter 3 Differentiation and Integration

Solutions to Problem Set 6

APPENDIX A Some Linear Algebra

COMPLEX NUMBERS AND QUADRATIC EQUATIONS

PES 1120 Spring 2014, Spendier Lecture 6/Page 1

PHYS 1101 Practice problem set 12, Chapter 32: 21, 22, 24, 57, 61, 83 Chapter 33: 7, 12, 32, 38, 44, 49, 76

Section 8.3 Polar Form of Complex Numbers

Math1110 (Spring 2009) Prelim 3 - Solutions

A Quantum Gauss-Bonnet Theorem

ˆ (0.10 m) E ( N m /C ) 36 ˆj ( j C m)

Mathematical Economics MEMF e ME. Filomena Garcia. Topic 2 Calculus

University of Washington Department of Chemistry Chemistry 452/456 Summer Quarter 2014

3 Basic boundary value problems for analytic function in the upper half plane

MTH 263 Practice Test #1 Spring 1999

OPTIMISATION. Introduction Single Variable Unconstrained Optimisation Multivariable Unconstrained Optimisation Linear Programming

Week 9 Chapter 10 Section 1-5

Phys102 General Physics II

A Simple Research of Divisor Graphs

Chapter Twelve. Integration. We now turn our attention to the idea of an integral in dimensions higher than one. Consider a real-valued function f : D

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0

More metrics on cartesian products

Module 2. Random Processes. Version 2 ECE IIT, Kharagpur

Work is the change in energy of a system (neglecting heat transfer). To examine what could

Spring Force and Power

Lecture 12: Discrete Laplacian

Bezier curves. Michael S. Floater. August 25, These notes provide an introduction to Bezier curves. i=0

University of Washington Department of Chemistry Chemistry 452/456 Summer Quarter 2014

Complex Numbers. x = B B 2 4AC 2A. or x = x = 2 ± 4 4 (1) (5) 2 (1)

Digital Signal Processing

U.C. Berkeley CS294: Beyond Worst-Case Analysis Luca Trevisan September 5, 2017

Lecture 3. Ax x i a i. i i

Modeling curves. Graphs: y = ax+b, y = sin(x) Implicit ax + by + c = 0, x 2 +y 2 =r 2 Parametric:

NUMERICAL DIFFERENTIATION

8.6 The Complex Number System

Modelli Clamfim Equazioni differenziali 7 ottobre 2013

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity

36.1 Why is it important to be able to find roots to systems of equations? Up to this point, we have discussed how to find the solution to

princeton univ. F 17 cos 521: Advanced Algorithm Design Lecture 7: LP Duality Lecturer: Matt Weinberg

Affine transformations and convexity

Implicit Integration Henyey Method

Example: (13320, 22140) =? Solution #1: The divisors of are 1, 2, 3, 4, 5, 6, 9, 10, 12, 15, 18, 20, 27, 30, 36, 41,

Exercise Solutions to Real Analysis

PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 125, Number 7, July 1997, Pages 2119{2125 S (97) THE STRONG OPEN SET CONDITION

Modelli Clamfim Equazione del Calore Lezione ottobre 2014

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X

SELECTED PROOFS. DeMorgan s formulas: The first one is clear from Venn diagram, or the following truth table:

Curvature and isoperimetric inequality

ECEN 5005 Crystals, Nanocrystals and Device Applications Class 19 Group Theory For Crystals

PHYS 705: Classical Mechanics. Calculus of Variations II

Integrals and Invariants of Euler-Lagrange Equations

Strong Markov property: Same assertion holds for stopping times τ.

C. J. PAPACHRISTOU ASPECTS OF

Appendix B. Criterion of Riemann-Stieltjes Integrability

Causal Diamonds. M. Aghili, L. Bombelli, B. Pilgrim

The Schrödinger Equation

9. Complex Numbers. 1. Numbers revisited. 2. Imaginary number i: General form of complex numbers. 3. Manipulation of complex numbers

Finding Dense Subgraphs in G(n, 1/2)

Spectral Graph Theory and its Applications September 16, Lecture 5

1 Matrix representations of canonical matrices

Appendix B. The Finite Difference Scheme

Bézier curves. Michael S. Floater. September 10, These notes provide an introduction to Bézier curves. i=0

Modelli Clamfim Equazioni differenziali 22 settembre 2016

Physics 5153 Classical Mechanics. Principle of Virtual Work-1

Difference Equations

Finite Difference Method

CHAPTER 4d. ROOTS OF EQUATIONS

Linear, affine, and convex sets and hulls In the sequel, unless otherwise specified, X will denote a real vector space.

Complex Numbers Alpha, Round 1 Test #123

Salmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2

General Tips on How to Do Well in Physics Exams. 1. Establish a good habit in keeping track of your steps. For example, when you use the equation

Module 1 : The equation of continuity. Lecture 1: Equation of Continuity

Force = F Piston area = A

Bernoulli Numbers and Polynomials

The Second Eigenvalue of Planar Graphs

Physics for Scientists & Engineers 2

8.022 (E&M) Lecture 4

NP-Completeness : Proofs

Lecture Note 3. Eshelby s Inclusion II

PHYS 1443 Section 004 Lecture #12 Thursday, Oct. 2, 2014

Indeterminate pin-jointed frames (trusses)

Calculus of Variations Basics

Conservation of Energy

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system

Affine and Riemannian Connections

Formulas for the Determinant

Lecture 21: Numerical methods for pricing American type derivatives

Lecture 10: Euler s Equations for Multivariable

b ), which stands for uniform distribution on the interval a x< b. = 0 elsewhere

Lecture 2 Solution of Nonlinear Equations ( Root Finding Problems )

1 (1 + ( )) = 1 8 ( ) = (c) Carrying out the Taylor expansion, in this case, the series truncates at second order:

Statistical Mechanics and Combinatorics : Lecture III

Problem Set 9 Solutions

PHYS 705: Classical Mechanics. Newtonian Mechanics

2.3 Nilpotent endomorphisms

9 Characteristic classes

Here is the rationale: If X and y have a strong positive relationship to one another, then ( x x) will tend to be positive when ( y y)

From Biot-Savart Law to Divergence of B (1)

STUDY PACKAGE. Subject : Mathematics Topic : COMPLEX NUMBER Available Online :

Transcription:

omplex Varables hapter 8 Integraton n the omplex Plane March, Lecturer: Shh-Yuan hen Except where otherwse noted, content s lcensed under a BY-N-SA. TW Lcense.

ontents ontour ntegrals auchy-goursat theorem Independence o path auchy s ntegral ormulas

ontour Integrals Recall that the de. o ( x) starts wth a real uncton y (x) that s dened on [a, b] on the x-axs. Ths can be generaled to ntegrals o real unctons o two varables dened on a curve n the artesan plane. b a dx

ontour Integrals For the ntegral o a complex x () that s dened along a curve n the complex plane. could be smooth curve, pecewse smooth curve, closed curve, or smple closed curve, and s called a contour or path. s dened n terms o parametrc equatons x x(t), y y(t), a t b, where t s a real parameter. Usng x(t) & y(t) as real & magnary parts, we can descrbe n the complex plane by means o a complex-valued uncton o a real varable t: (t) x(t) y(t), a t b. 4

ontour Integrals An ntegral o () on s denoted by or ( ) d s closed. It s reerred to as a contour ntegral or a complex lne ntegral. ( ) d 5

ontour Integrals () u(x, y) v(x, y) Let be dened at all ponts on a smooth curve dened by x x(t), y y(t), a t b. Dvde nto n subarcs accordng to the partton a t < t < < t n b. orrespondng ponts on : x y x(t ) y(t ), x y x(t ) y(t ),, n x n y n x(t n ) y(t n ). Let k k k, k,,, n & let P be the norm o the partton. hoose a pont on each subarc. Form the sum * * k xk n ( ) k y k * k * 6 k

ontour Integrals De. Let be dened at ponts o a smooth curve dened by x x(t), y y(t), a t b. The contour ntegral o along s ( ) ( * d lm ) ( ) P k n The lmt exsts s contnuous at all ponts on & s ether smooth or pecewse smooth. k k 7

ontour Integrals Method o evaluaton Wrte () n the abbrevated orm ( ) d lm ( u v)( x y) lm { ( u x v y) ( v x u y) } u dx v dy Thus, a contour ntegral s a combnaton o two real-lne ntegrals. v dx u dy ( ) 8

ontour Integrals Snce x x(t), y y(t), a t b, the RHS o () s b [ u ( x( t), y( t) ) x ( t) v( x( t), y( t) ) y ( t) ] a b a I we use (t) x(t) y(t) to descrbe, the last result s the same as b ( ( t) ) ( t) separated nto two ntegrals. a dt dt [ v( x( t), y( t) ) x ( t) u( x( t) y( t) ) y ( t) ], when dt 9

ontour Integrals Thm. I s contnuous on a smooth curve gven by (t) x(t) y(t), a t b, then ( ) d ( ( t) ) ( t) dt ( ) b a I s expressed n terms o, then to evaluate () we smply replace by (t). I not, we replace x & y wherever they appear by x(t) & y(t), respectvely.

ontour Integrals Ex. Evaluate d, where s gven by x t, (sol) y t, t 4. ( ) t t t ( t) ( ( t) ) d t 4 4 t t t ( t t )( t) dt t ( ) 4 9t t dt t dt 95 65

ontour Integrals Ex. Evaluate d, where s the crcle x cos t, y sn t, t π. (sol) ( ) t t cost sn t e ( t) ( ( t) ) d π e t e ( ) π t t e e dt dt π t

ontour Integrals Propertes o contour ntegrals Suppose & g are contnuous n a doman D & s a smooth curve lyng entrely n D. ( ) k ( ) d k ( ) ( ) [ ( ) g( ) ] d ( ) d g( ) ( ) ( ) d ( ) d ( ) ( v) ( ) d ( ) d, d k a constant. d, d

ontour Integrals ( x y ) Ex. Evaluate, where. d ( ) ( ) ( x y d x y d x y ) s dened by y x t, t. ( ) ( x y d t t )( ) s dened by x, y t, t. ( ) ( ) ( ) x y d t dt t dt dt t dt dt d 7

( ) ML d ontour Integrals Thm. I s contnuous on a smooth curve & () M or all on, then, where L s the length o. (proo) From (5) o h.7, we can wrte k : length o the chord jonng the ponts k & k. ( ) ( ) n k k n k k k n k k k M * * ( ) ML M n k k n k k k * 5

ontour Integrals e Ex. Fnd an upper bound or d, where s the crcle 4. The length o the crcle s 8π. From the nequalty o h.7, e e e e x ( cos y sn y) e x e 4 e d 8πe 4 6

ontour Integrals rculaton & net lux T & N denote the unt tangent vector & unt normal vector to a postvely orented smple closed contour. Interpret () u(x,y) v(x,y) as a vector, the lne ntegrals T ds u dx v dy ( 4) N ds u dy v dx ( 5) The lne ntegral n (4) s called crculaton around & measures the tendency o the low to rotate.

ontour Integrals The net lux across s the derence between the rate at whch lud enters & the rate at whch lud leaves the regon bounded by. It s gven by the lne ntegral n (5). Note that ( ) ( ) T ds N ds ( u v)( dx dy) ( ) ( ( ) ) rculaton Re d ( 6) ( ( ) ) Net lux Im d ( 7) d 8

ontour Integrals Ex. Gven the low () ( ), compute the crculaton around & the net lux across the crcle :. ( ) ( ) ( t) e t, ( ) d ( ) π t e t e dt ( ) dt π ( ) π π t 9

auchy-goursat Theorem Focus on contour ntegrals where s a postve smple closed curve Smply & multply connected domans Smply connected every smple closed lyng entrely n D can be shrunk to a pont wthout leavng D. Ex. The entre complex plane. Multply connected a doman that s not smply connected.

auchy-goursat Theorem auchy s Theorem Suppose that s analytc n a smply connected doman D & that s contnuous n D. For every smple closed contour n D, ( ) d (proo) Snce s contnuous throughout D, the real & magnary parts o () u v & ther st partal dervatves are contnuous n D. From (), ( ) d u( x, y) dx v( x, y) dy v( x, y) dx u( x, y) v x u y da D D u x v y da dy

auchy-goursat Theorem auchy-goursat Theorem Suppose s analytc n a smply connected doman D. Then or every smple closed contour n D, ( ) d ( 8) In other words, s analytc at all ponts wthn & on a smple closed contour, then e d Ex. Evaluate. ( ) d Snce e s entre & s a smple closed contour, t ollows rom the auchy- Goursat theorem that e d

auchy-goursat Theorem d ( x ) Ex. Evaluate, where : () / s analytc everywhere except at. But s not nteror to or on. ( ) cos Ex. Gven the low, compute the crculaton around & net lux across, where s the square wth vertces ± & ±. rculaton & net lux are both ero. ( ) d cos d ( y 5) 4 d

auchy-goursat Theorem For multply connected domans I s analytc n a multply connected doman D, then we cannot conclude that ( ) d or every smple closed n D. Suppose D s doubly connected and & are smple closed contours such that surrounds the hole n D & s nteror to. Also, suppose s analytc on each contour & at each pont nteror to but exteror to.

auchy-goursat Theorem As we ntroduce the cut AB, the regon bounded by the curves s smply connected. ( ) d ( ) d ( ) d ( ) ( ) d ( ) d ( 9) AB (9) s called the prncple o deormaton o contours. Thus, one can evaluate an ntegral over a complcated smple closed contour by replacng that contour wth one that s more convenent. BA d 5

auchy-goursat Theorem Ex. Evaluate shown n the gure. From (9), we choose. Takng the radus o the crcle to be r, les wthn. : Let e t, t π e d t d d, where s the outer contour π e e t t dt π

auchy-goursat Theorem I s any constant complex number nteror to any smple closed contour, then d π, n ( ) ( ) n, n For n s ero or negatve, the ntegral s ero. For n s postve nteger derent rom one, t re t :, π and t d d π re dt r n n nt r e ( ) n ( ) n re e t ( n ) t ( n) π 7

auchy-goursat Theorem Ex. Evaluate wth :. d 7 5 ( )( ) d d d π 6π 7 5 7 5 7 5 8

auchy-goursat Theorem I,, are smple closed contours & s analytc on them & at ponts nteror to but exteror to both &, then by ntroducng two cuts, ( ) d ( ) d ( ) d

auchy-goursat Theorem G Thm or Multply onnected Domans Suppose,,, n are postve smple closed contours such that,,, n are nteror to but ther nteror regons have no ponts n common. I s analytc on each contour & at each pont nteror to but exteror to all k, k n, n ( ) d ( ) d ( ) k k

auchy-goursat Theorem Ex. Evaluate, where :. Surround ± by crcle contours & that le entrely wthn. d d d : and :

d d d d d d d π π auchy-goursat Theorem

auchy-goursat Theorem It can be shown that the auchy-goursat theorem s vald or any closed n a smply connected doman D. As shown n the gure, s analytc n D, then ( ) d

Independence o Path De. Let & be ponts n D. A contour ntegral ( ) d s sad to be ndependent o the path ts value s the same or all contours n D wth an ntal pont & a termnal pont. Suppose & are contours n a smply connected doman D, both rom to. Note that & orm a closed contour. Thus, s analytc n D, ( ) d ( ) d ( ) d ( ) d ( )

Independence o Path Analytcty mples path ndependence Thm. I s an analytc uncton n a smply connected doman D, then ( ) d s ndependent o. Ex. Evaluate d where s the contour rom to shown n the gure. Snce () s entre, we can replace by any, say : x, y t, t d d ( t) dt 5

Independence o Path ( ) A contour ntegral that s ndependent o rom to s wrtten. De. Suppose s contnuous n D. I there exsts a uncton F such that F () () or n D, then F s called an antdervatve o. The most general antdervatve o () s wrtten ( ) d F( ) Snce F has a dervatve at each pont n D, t s necessarly analytc & hence contnuous n D. d ( ) d 6

Independence o Path Fundamental Theorem o ontour Integrals Thm. Suppose s contnuous n D & F s an antdervatve o n D. Then or any n D rom to, (proo) ( ) d F( ) F( ) ( ) As s a smooth curve dened by (t), a t b. Use () & the act that F () () or n D, ( ) d ( ( t) ) ( t) dt F ( ( t) ) ( t) dt F( ( t) ) F b a b d dt ( ( t) ) F( ( b) ) F( ( a) ) F( ) F( ) a b a b a 7 dt

Independence o Path Ex. Use () to solve the prevous example. Ex. Evaluate d rom to. d ( ) ( ) cos cos d sn d ( ).4, where s any contour cos d sn sn.489 8

Independence o Path I s closed, then & consequently ( ) d ( 4) ( ) d ( ) Snce the value o d &, t s the same or any n D connectng these ponts. In other words, depends on only I s contnuous & has an antdervatve F n D, then s ndependent o. ( ) I s contnuous & s ndependent o n D, then has an antdervatve everywhere n D. d

Independence o Path (proo) Assume that s contnuous, ( ) ndependent o n D, & F s a uncton dened by F ( ) ( s) ds s where s s a complex varable, s a xed pont n D, & s any pont n D. F hoose so that s n D and & can be joned by a straght segment n D. d ( ) F( ) ( s) ds ( s) ( s) ds ds 4

Independence o Path (proo) As xed, we can wrte ( ) ( ) ds ( ) ( ) ( ) ( ) ( ) ds F F ( ) [ ( s) ( ) ] ds Now s contnuous at the pont. It means that or any ε > there exsts a δ > so that (s) () < ε whenever s < δ. ds 4

Independence o Path F (proo) onsequently, we choose so that < δ, F ( ) [ ( s) ( ) ] ds ( ) ( ) lm F ( ) F( ) [ ( s) ( ) ] ds ε ε ( ) or F ( ) ( ) 4

Independence o Path It s known that I s analytc n a smply connected doman D, s necessarly contnuous throughout D. Its contour ntegral s ndependent o the path. ombne wth the result just obtaned, Thm. I s analytc n a smply connected doman D, then has an antdervatve F everywhere n D or there exsts a uncton F so that F () () or all n D. 4

Independence o Path Recall rom h.7 that / s the dervatve o Ln. It means that under some crcumstances Ln s an antdervatve o /. Suppose D s the entre complex plane wthout the orgn. / s analytc n D. I s any smple closed contour contanng the orgn, In ths case, Ln s NOT an antdervatve o / n D, snce Ln s not analytc n D. (Ln als to be analytc on the nonpostve real axs!) d π 44

Independence o Path Ex. Evaluate wth shown n the gure. D s the smply connected doman dened by x d ( ) >, y Im( ) Re > Here, Ln s an antdervatve o /, snce both Ln & / are analytc n D. d Ln Ln ( ) π loge loge.455.578 Ln 45

Independence o Path I & g are analytc n a smply connected doman D contanng rom to, then the ntegraton by parts ormula s vald n D: ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) d g g d g d g d g d g d d g d g d d 46

auchy s Integral Formulas auchy s ntegral ormula Thm. Let be analytc n a smply connected doman D and be a smple closed contour lyng entrely wthn D. I s any pont wthn, then ( ) ( ) ( 5) d π (proo) Let be a crcle centered at wth radus small enough that t s nteror to. 47

auchy s Integral Formulas By the prncple o deormaton o contours, ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) d d d d d d π 48

auchy s Integral Formulas Snce s contnuous at, or any ε >, there exsts a δ > such that () ( ) < ε whenever < δ. Thus, we choose to be δ/ < δ, then by the ML-nequalty, ( ) ( ) d ε δ δ π πε In other words, the absolute value o the ntegral can be made arbtrarly small by takng the radus o to be sucently small. Ths can happen only the ntegral s ero. (5) can thus be obtaned. 49

auchy s Integral Formulas Snce most problems do not have smply connected domans, a more practcal restatement o auchy s ntegral theorem: ( ) I s analytc wthn & on a smple closed contour & s nteror to, ( ) 4 4 d ( ) d π ( ) d π Ex. Evaluate wth :. 4 4 ( 5) ( ) π ( 4 )

auchy s Integral Formulas Ex. Evaluate wth : 4. ( ) ( ) d d π π π 6 9 and 9 d 9 5

auchy s Integral Formulas ( ) k ( ) Ex., where k a b & are complex numbers, gves rse to a low n the doman. I s a smple closed contour contanng n ts nteror, then a b ( ) d d π ( a b) Thus, the crculaton around s πb, & the net lux across s πa. I were n the exteror o, both the crculaton & the net lux would be ero by auchy s theorem. 5

auchy s Integral Formulas When k s real, the crculaton around s ero but the net lux across s πk. The complex number s called a source or the low when k > & a snk when k <. 5

auchy s Integral Formulas auchy s ntegral ormula or dervatves Thm. Let be analytc n a smply connected doman D and be a smple closed contour lyng entrely wthn D. I s any pont wthn, then ( n ) n! ( ) ( 6) d n π (proo) ( ) lm lm ( ) ( ) ( ) ( ) π ( ) ( ) d ( ) d

auchy s Integral Formulas Snce s contnuous on, there exsts a real number M so that () M or on. Besdes, let L be the length o and δ be the shortest dstance between ponts on &. Thus, or all on, I we choose δ/, then ( ) ( ) ( )( ) d lm π or δ δ δ δ δ and

auchy s Integral Formulas ( ) ( )( ) ( ) ( ) ( ) d ( ) ( ) δ ML Snce the last expresson goes to ero as, d ( ) ( ) ( ) π d d 56

auchy s Integral Formulas I () u(x,y) v(x,y) s analytc at a pont, then ts dervatves o all orders exst at that pont & are contnuous. onsequently, rom u & v have contnuous partal dervatves o all orders at a pont o analytcty. ( ) ( ) x y u x y v x v x u y u y v x v x u 57

auchy s Integral Formulas Ex. Evaluate wth :. The ntegrand s not analytc at & 4, but only les wthn. By (6) we have d 4 4 ( ),, 4 4 4 4 n ( ) ( ) d 4 6! 4 4 π π π 58

auchy s Integral Formulas Ex. Evaluate wth shown n the gure. Thnk o as the unon o two smple closed contours &. ( ) d ( ) ( ) ( ) ( ) ( ) I I d d d d d 59

auchy s Integral Formulas Use (5) & (6) or I & I, respectvely. ( ) ( ) ( ) ( ) ( ) I I d d I d I π π π π π π π π 4 6 4! 6 6

auchy s Integral Formulas Louvlle s Theorem The only bounded entre unctons are constants. Suppose s entre & bounded, that s, () M or all. Then or any pont, ( ) M/r. ( n ) ( ) n! ( ) n! M n! d πr ( ) n n n π r r π By takng r arbtrarly large, we can make ( ) as small as we wsh. ( ) or all ponts n the complex plane. Hence must be a constant. M 6

auchy s Integral Formulas Fundamental theorem o algebra I P() s a nonconstant polynomal, then the equaton P() has at least one root. Suppose that P() or all. Ths mples that the recprocal o P, () /P(), s an entre uncton. Snce () as, must be bounded or all nte. It ollows rom Louvlle s thm that s a constant & thereore P s a constant. ontradcton! We conclude that there must exst at least one 6 number or whch P().

opyrght Page Work Lcense Author/Source 4 Shh-Yuan hen 9 Shh-Yuan hen Shh-Yuan hen Shh-Yuan hen 6

opyrght Page Work Lcense Author/Source Shh-Yuan hen Shh-Yuan hen 4 Shh-Yuan hen 5 Shh-Yuan hen 64

opyrght Page Work Lcense Author/Source 6 Shh-Yuan hen 9 Shh-Yuan hen Shh-Yuan hen Shh-Yuan hen 65

opyrght Page Work Lcense Author/Source 4 Shh-Yuan hen 5 Shh-Yuan hen 4 Shh-Yuan hen 45 Shh-Yuan hen 66

opyrght Page Work Lcense Author/Source 5 Shh-Yuan hen 5 Shh-Yuan hen 59 6 Shh-Yuan hen 67