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.

SERIES SOLUTIONS OF ORDINARY DIFFERENTIAL EQUATIONSWe may check that y 1 (z) <strong>and</strong>y 2 (z) are indeed linearly independent by computing theWronskian as follows:W = y 1 y 2 ′ − y 2 y 1′=sin √ (z − 12 √ z sin √ )z − cos √ ( 1z2 √ z cos √ )z= − 12 √ ( √ sin2z +cos √ 2 z ) = − 1z2 √ z ≠0.Since W ≠0,thesolutionsy 1 (z) <strong>and</strong>y 2 (z) are linearly independent. Hence, the generalsolution to (16.17) is given byy(z) =c 1 sin √ z + c 2 cos √ z. ◭16.3.2 Repeated root of the indicial equationIf the indicial equation has a repeated root, so that σ 1 = σ 2 = σ, then obviouslyonly one solution in the <strong>for</strong>m of a Frobenius series (16.12) may be found asdescribed above, i.e.y 1 (z) =z σ∞ ∑n=0a n z n .<strong>Methods</strong> <strong>for</strong> obtaining a second, linearly independent, solution are discussed insection 16.4.16.3.3 Distinct roots differing by an integerWhatever the roots of the indicial equation, the recurrence relation correspondingto the larger of the two always leads to a solution of the ODE. However, if theroots of the indicial equation differ by an integer then the recurrence relationcorresponding to the smaller root may or may not lead to a second linearlyindependent solution, depending on the ODE under consideration. Note that <strong>for</strong>complex roots of the indicial equation, the ‘larger’ root is taken to be the onewith the larger real part.◮Find the power series solutions about z =0ofz(z − 1)y ′′ +3zy ′ + y =0. (16.20)Dividing through by z(z − 1) to put the equation into st<strong>and</strong>ard <strong>for</strong>m, we obtainy ′′ + 3 1(z − 1) y′ + y =0, (16.21)z(z − 1)<strong>and</strong> on comparing with (16.7) we identify p(z) =3/(z − 1) <strong>and</strong> q(z) =1/[z(z − 1)]. Weimmediately see that z = 0 is a singular point of (16.21), but since zp(z) =3z/(z − 1) <strong>and</strong>z 2 q(z) =z/(z−1) are finite there, it is a regular singular point <strong>and</strong> we expect to find at least542

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

Saved successfully!

Ooh no, something went wrong!