AST-V57-Sem1-S24014 Differential Equations
AST-V57-Sem1-S24014 Differential Equations
AST-V57-Sem1-S24014 Differential Equations
Yashwantrao Chavan
Maharashtra Open University
Semester I
Code : S24014
Differential Equations
For
Developed by
K J
Copyright © Registrar,
Shivaji University,
Kolhapur. (Maharashtra)
First Edition 2008
Second Edition 2010
All rights reserved, No part of this work may be reproduced in any form by mimeography
or any other means without permission in writing from the Shivaji University, Kolhapur
(MS)
Copies : 1000
Published by:
Dr. D. V. Muley
Registrar,
Shivaji University,
Kolhapur-416 004
H Further information about the Centre for Distance Education & Shivaji University may be
obtained from the University Office at Vidyanagar, Kolhapur-416 004, India.
H This material has been produced with the developmental grant from DEC-IGNOU, New
Delhi.
ISBN- 978-81-8486-012-2
p Production :
Shri. Anand Yadav
Head, Print Production Centre
Yashwantrao Chavan Maharashtra Open University
Nashik-422 222
p Republished by :
Dr. Prakash Atkare
Registrar, Yashwantrao Chavan Maharashtra Open University
Nashik-422 222
p Reprinted by :
Printed by : M/s. Mr. Shashikant Ahirrao,
Abhiyankit Printers, B-11, NICE Area,
Satpur MIDC, Nashik - 422 007.
Anurath/AB15-204 (ii)
Centre for Distance Education
Shivaji University, Kolhapur
n ADVISORY COMMITTEE n
Prof. (Dr.) N. J. Pawar Dr. J. S. Patil
Vice-Chancellor, Dean, Faculty of Social Sciences,
Shivaji University, Kolhapur Shivaji University, Kolhapur
Prof. P. Prakash Dr. C. J. Khilare
Hon'ble Vice-Chancellor, Dean, Faculty of Science,
Dr. Babasaheb Ambedkar Open University, Shivaji University, Kolhapur
Hyderabad. Dr. R. G. Phadatare
Prof. (Dr.) K. S. Rangappa Dean, Faculty of Commerce,
Hon. Vice-Chancellor, Shivaji University, Kolhapur
University of Mysore Prof. (Dr.) A. B. Rajage
Prof. (Dr.) R. Krishna Kumar Director, B.C.U.D.,
Hon. Vice-Chancellor, Shivaji University, Kolhapur
Yashwantrao Chavan Maharashtra Prof. (Dr.) D. V. Muley
Open University, Dnyangangotri, Registrar,
Near Gangapur Dam, Nasik Shivaji University, Kolhapur
Prin. (Dr.) A. S. Bhoite Dr. B. M. Hirdekar
Pro-Vice-Chancellor, Controller of Examinations,
Shivaji University, Kolhapur Shivaji University, Kolhapur
Prof. (Dr.) Cima Yeole Shri. V. T. Patil,
Git Govind, Flat No. 2, Finance and Accounts Officer,
1139 Sykes Extension, Shivaji University, Kolhapur
Kolhapur-416001 Prof. (Dr.) A. R. Bhosale
Dr. A. P. Gavali (Member Secretary)
Dean, Faculty of Arts and Fine Arts, Director, Centre for Distance Education,
Shivaji University, Kolhapur Shivaji University, Kolhapur.
n B. O. S. MEMBERS OF MATHEMATICS n
Chairman- Dr. H. T. Dinde
Karmveer Bhaurao Patil College, Urun-Islampur, Tal. Walwa, Dist. Sangli
l Dr. L. N. Katkar l Dr. S. D. Thikane
Head, Dept. of Mathematics, Jaysingpur College, Jaysingpur, Tal. Shirol,
Shivaji University, Kolhapur. Dist. Kolhapur
l Dr. H. G. Datar l Dr. S. P. Thorat
Willingdon College, Sangli S. G. M. College, Karad
l Dr. D B. Dhaigude
l Dr. T. B. Jagtap
Department of Mathematics, Dr. Babasaheb Ambedkar
Yashwantrao Chavan Institute of Science, Satara.
Marathwada University, Aurangabad
l Shri. L. B. Jamale
l Dr. Mrs. S. H. Thakar
Krishna Mahavidyalaya, Rethare Bk., Karad, Dist.
Department of Mathematics,
Satara.
Shivaji University, Kolhapur
l Dr. A. D. Lokhande
l Mr. Sarjerao R. Patil
Yashavantrao Chavan Warana Mahavidyalaya,
Symantec Software India Pvt. Ltd., West End
Warananagar, Dist. Kolhapur
Centre three-survey No. 169/1 Aundh,
Pune-411007.
(iii)
Centre for Distance Education Differential Equations
Shivaji University,
Kolhapur.
n Editor n
Dr. (Mrs.) Sarita Thakar
Department of Mathematics,
Shivaji University, Kolhapur.
Maharashtra.
(iv)
Preface
Large numbers of students appear for M.A./M. Sc. Examinations externally every
year. In view of this, Shivaji University has introduced the Distance Education Mode for
external students from the year 2007-2008, and entrusted the task to us to prepare the
Self Instructional Material (SIM) for aspirants.
It is hoped that students must learn Mathematics not only to become competent
mathematicians but also skilled users of Mathematics in the solution of problems in the
real world. They must learn how to use their Mathematical knowledge in solving the
problems of the real world. Differential equations usually are description of physical
systems. This book on Differential Equations consists of four chapters. Chapter one
contains the complete discussion of linear equations with constant coefficients, including
the uniqueness theorem. In chapter two linear equations with variable coefficients are
trea. Equations with analytic coefficients are introduced and series solutions are
obtained by a simple formal process. A detailed treatment of linear equations with
regular singular points is discussed in chapter four. Classification of regular singular
points and regular singular points at infinity is studied. In chapter five existence and
uniqueness of solutions of first order initial value problem are established. The
innumerable examples and exercises are given at the end of each unit.
The book introduces the students to some of the abstract topics that pervade
modern analysis. The first chapter deals with the Riemann Stieltjes integration. The
problems in Physics and Chemistry which involve mass distribution that are partly
discrete and partly continuous can be solved by using Riemann Stietjes integrations.
The Chapter 2 deals with convergence and uniform convergence of sequences of
functions and series where as the Chapter 3 consists of multidimensional calculus.
The Chapter 4 deals with implicit functions and extremum problems which have wide
applications in optimization theory. Line integrals, surface integrals and Volume integrals
are the subject matter of Chapter 5. This provides sufficient background to study the
Gauss divergence Theorem and Stokes Theorem.
Editor
(v)
(vi)
M. Sc. (Mathematics)
Differential Equations
Contents
Chapter 1 : Linear Equations with Constant Coefficients 1
(vii)
Each Unit begins with the section objectives -
Dear Students
The SIM is simply a supporting material for the study of this paper.
It is also advised to see the new syllabus 2013-14 and study the
reference books & other related material for the detailed study of the
paper.
(viii)
p p
M. Sc. (Mathematics)
Paper III
Differential Equations
❁ D r. Sarita Thakar ❁
Department of Mathematics
Shivaji University, Kolhapur (M.S.)
p p
Differential Equations
Chapter
1
Linear Equations with
Constant Coefficients
Contents :
Unit 1 : Initial value problems for second order equations.
Unit 2 : Linear dependence and independencce
Unit 3 : The homogenous equation of order n
Unit 4 : The non-homogeneous equation of order n
Introduction :
We live in a world of interrelated changing entities. The position of the earth changes with
time, the velocity of falling body changes with distance, the bending of a beam changes with the
weight of the load placed on it, the area of circle changes with the size of the radius, the path of
projectile changes with the velocity and angle at which it is fired.
In the language of mathematics changing entities are called variables and the rate of change
of one variable with respect to another is called derivative. Equations which express a relation
among these variables and their derivatives are called differential equations.
A Linear differential equation of order n with constant coefficients is an equation of the
form
a0 y ( n ) + a1 y ( n –1) + a2 y ( n –2) + ⋅⋅⋅ + an y = b( x),
where, a0 ≠ 0, a1, a2 , ⋅⋅⋅, an are complex constants
and b is complex valued function on an interval I : a < x < b .
The operator L defined by
Theorem 1.1.1
Let, a1, a2 be constants and consider the equation L(y) = y¢¢ + a1y¢ + a2 y = 0
1. If r1, r2 are distinct roots of the characteristic polynomial
p(r) = r2 + a1r + a2
then the functions φ1 ( x ) = e r1x and φ2 ( x ) = e r2 x are solutions of L(y) = 0.
2. If r1 is a repeated root of the characteristic polynomial p(r), then the functions φ1 ( x) = e r1x
rx
and φ2 ( x ) = xe 1 are solutions of L(y) = 0.
Proof : Let f (x) = erx be a solutions of L(y) = 0.
L(e rx ) = (e rx )′′ + a1 (e rx )′ + a2 e rx
= (r 2 + a1r + a2 )e rx
L (erx) = 0 if and only if p(r) = r2 + a1r + a2 = 0.
1. If r 1 and r 2 are distinct roots of p(r) then L(e r1x ) = L(e r2 x ) = 0 and φ1 ( x ) = e r1x and
φ2 ( x) = e r2 x are solutions of L(y) = 0.
2. If r1 is a repeated root of p(r) then
P(r ) = (r – r1 ) 2 and p′(r ) = 2(r – r1 )
Theorem 1.1.2 :
If f1 and f2 are two solutions of L(y) = 0 then C1 f1 + Cf2 is also a solution of L(y) = 0.
Where, C1 and C2 are any two constants.
Proof : Let f1 and f2 be two solutions of L(y) = 0
L (φ1 ) = φ1′′ + a1 φ1′ + a2 φ1 = 0
Differential Equations (2)
L(φ2 ) = φ2′′ + a1 φ2′ + a2 φ1 = 0
Suppose C1 and C2 are any two constants then the function f defined by f = C1 f1 + C2 f2
is also a solution of L(y) = 0.
L(φ ) = (aφ1 + c2φ 2 )′′ + a1 (aφ1 + c2φ 2 ) + a2 (aφ1 + c2φ 2 )
= c1L(φ1 ) + c2 L (φ2 )
=0
The function f which is zero for all x is also a solution called the trivial solution
of L(y) = 0.
The results of above two theorems allow us to solve all homogeneous linear second order
differential equations with constant coefficients.
Definition 1.1 :
An initial value problem L(y) = 0 is a problem of finding a solution f satisfying
φ ( x0 ) = α 0 and φ ′( x0 ) = β 0 where, x0 is some real number and a0, b0 are given constants.
L( y ) = y ′′ + a1 y ′ + a2 y = 0, y ( x0 ) = α , y ′( x0 ) β , – ∞ < x < ∞.
Proof : By theorem 1.1.1 there exist two solutions f1 and f2 that satisfy L(f1) = L(f2) = 0. From
theorem 1.1.2 we know that c1 f1 + c2 f2 is a solution of L(y) = 0. We show that there are
unique constants c1, c2 such that φ = c1φ1 + c2φ2 satisfies φ ( x0 ) = α and φ ′( x0 ) = β .
φ ( x0 ) = c1φ1 ( x0 ) + c2φ 2 ( x0 ) = α
φ ′( x0 ) = c1φ1′ ( x0 ) + c2φ 2′ ( x0 ) = β
Above system of equations will have a unique solution c1, c2 if the determinant
φ1 ( x0 ) φ 2 ( x0 )
∆= = φ1 ( x0 ) φ2′ ( x0 ) – φ2 ( x0 ) φ1′ ( x0 ) ≠ 0.
φ1′ ( x0 ) φ 2′ ( x0 )
By theorem 1.1.1 (1), φ1 ( x ) e r1 x and φ2 ( x) e r2 , x are two solution of L( y ) = 0 for r1 ≠ r2
and
∆ = e r1 x0 r2 e r2 x0 – e r2 x0 r1 e r1x0
= (r2 – r1 )e(r1 +r2 ) x0 ≠ 0.
Defination 1.2 :
A solution of a differential equation will be called a particular solution if it satisfies the
equation and does not contain arbitrary constants.
Theorem 1.1.4 :
Let, f be any solution of
L( y ) = y′′ + a1 y ′ + a2 y = 0
on an interval I containing a point x0, Then for all x in I.
|| φ ( x0 ) || e – k | x – x0 | ≤ || φ ( x) || ≤ || φ ( x0 ) || e k | x – x0 |
Where,
1
φ ( x ) = | φ ( x ) |2 + | φ ′( x ) |2
2
and k = 1+ | a1 | + | a2 | .
Proof : Let,
u ( x) = || φ ( x) ||2
= | φ ( x ) |2 + | φ ′( x ) |2
= φ ( x) φ ( x) + φ ′( x) φ ′( x)
Then, u ′( x) = φ ′( x) φ ( x) + φ ( x) φ ′( x ) + φ ′′( x) φ ′( x) + φ ′( x) φ ′′( x)
and | u ′( x) | ≤ 2 | φ ( x ) | | φ ′( x) | +2 | φ ′( x ) | | φ ′′( x ) |
as | φ ( x ) | = | φ ( x ) |
Since f is a solution of L(y) = 0, L(φ ) = φ ′′ + a1φ ′ + a2φ = 0
But, 2| φ ( x) || φ ′( x) | ≤ | φ ( x) |2 + | φ ′( x) |2
Therefore,
| u ′( x) | ≤ 2 (1+ | a1 | + | a2 | )| φ ′( x ) |2 + 2 (1 + | a2 | )| φ ( x ) |2
≤ 2 (1+ | a1 | + | a2 |) | φ ′( x) |2 + | φ ( x) |2
≤ 2k u ( x)
Thus, we get
–2u ( x) ≤ u ′( x) ≤ 2ku ( x )
u ′( x) ≤ 2ku ( x) is equivalent to u′( x ) – 2k u ( x) ≤ 0 since exponential functions are positive
on multiplying above inequality by e–2kx we get
|| θ ( x0 ) ||2 = | θ ( x0 ) |2 + | θ ′( x0 ) 2 | = 0
Theorem 1.1.6 :
Let f 1, f 2 be two solutions of L(y) = 0 given by theorem 1.1.1. If c1, c2 are any two
constants the function f = c1 f 1 + c2 f 2 is a solution of L(y) = 0 on – ¥ < x <¥.
Conversely, if f is any solution of L(y) = 0 on – ¥ < x <¥, then there are unique constants
C1 and C2 such that f = C1 f 1+ C2 f 2.
Proof : First part of the theorem follows from theorem 1.1.2.
Conversely suppose f is any solution of L(y) = 0. Let φ ( x0 ) = α and φ ′( x0 ) = β for
some constants a and b. In the proof of existence theorem 1.1.3 we showed that there is a
solution y of L(y) = 0 satisfying.
ψ ( x0 ) = α , ψ ′( x0 ) = β of the form
ψ ( x ) = c1 φ1 ( x ) + c2φ2 ( x )1 where c1 and c2 are uniquely determined by a and b . By
uniqueness theorem 1.1.5 φ = ψ , for all x.
Examples :
1. Find all solutions of the following equations.
(a) y²– 4 y = 0 (b) y²+ 2 i y¢ + y = 0 (c) y²– 4 y¢ + 5y = 0
Answer :
(a) The characteristic polynomial is p(r) = r2 – 4. r1 = 2 and r2 = – 2 are two distinct roots of
p (r) = 0.
Therefore φ1 ( s ) = e 2 x and φ2 ( x) = e –2 x are two solutions. For any constants c1 and
c2, c1e2x + c2e–2x is a solution. Thus the general solution is φ1 ( x ) = c1e 2 x + c2e –2 x .
(b) The characteristic polynomial p(r) = r2 + 2ir + 1
1
p( r ) = 0 ⇒ r= –2i ± (2i) 2 – 4
2
1
= –2i ± –8
2
= –i ± 2 i
(
= –1 ± 2 i )
Therefore φ1 ( x) = e
( –1+ 2 )ix and φ2 ( x) = e
( –1– 2 )ix are two solutions. Thus, for any
constants A & B.
π
(b) φ ′′ + φ = 0, φ (0) = 1, φ = 0
2
Answer :
(a) The characteristic polynomial p(r) = r2 + r – 6. r1 = 2 and r2 = – 3 are distinct roots
φ ( x ) = c1e 2 x + c2e –3x is a general solution.
φ (0) = 1 ⇒ c1 + c2 = 1 .......... (1)
φ ′ (0) = 0 ⇒ φ ′ ( x) = 2c1e 2 x – 3 c2e –3 x at x = 0, gives φ ′ (0) = 2c1 – 3c2 = 0 .......... (2)
solving equation (1) and (2) for c1 and c2 we get c1 = 3/5 and c2 = +2/5.
3e 2 x 2 e3 x
Thus, the required solution is φ ( x ) = + .
5 5
Case 1. k > 0
k ix
∴ φ ( x) = c1 e + c2 e – k ix is a general solution
In general φ ( x) = A cos k x + B sin k x is a solution.
φ (0) = 0 ⇒ A cos 0 + B sin 0 = 0 i.e. A = 0
φ (π ) = 0 ⇒ A cos π + B sin π = 0 i.e. A = 0
Case 2. k = 0
p(r ) = r 2 = 0 ⇒ r = 0 a repeated root.
∴ φ ( x ) = c1e 0 + c2 xe 0 = c1 + c2 x is a solution
φ (0) = 0 ⇒ c1 = 0
φ (π ) = 0 ⇒ c1 + c2π = 0 ⇒ c2 = 0
Therefore there is no nontrivial solution corresponding to k = 0.
Case 3. k < 0
for k = 0, p (r) = r2 + k has distinct roots
r1 = – k & r2 = – – k ( Since k < 0, – k > 0)
–k x
φ (x) = c1e + c 2e – –k x
φ (0) = c1 + c2 = 0
–k π –k π
φ (π ) = c1e + c2e – =0
Simultaneous evaluation of above two equations give c1 = c2 = 0.
Thus, there is no non-trival solution corresponding to k < 0.
The only non-trivial solution for the given equation is φ ( x ) = B sin k x.
(d) The characteristic polynomial p(r) = r2 – 2r – 3
r1 = 3, r2 = 1 are two distinct roots.
Differential Equations (8)
\ φ ( x) = c1e3 x + c2e – x is a general solution
φ (0) = 0 ⇒ φ (0) = c1 + c2 = 0
φ ′ ( x) = 3c1 e3 x – c2 e – x
φ ′ (0) = 1 ⇒ φ ′ (0) = 1 = 3c1 – c2
Thus, c1 + c2 = 0 and 3c1 – c2 = 1 gives
1 1
c1 =
and c2 = –
4 4
1 3x 1 – x
Therefore φ ( x) = e – e is the required solution.
4 4
EXERCISES
Defination 1.4 : Given the functions f1, f 2, f3, ⋅⋅⋅, f n if constants c1, c2 , c3 , ⋅⋅⋅, cn not all zero
exist such c1 f1 ( x ) + c2 f 2 ( x) + c3 f 3 ( x) + ⋅⋅⋅ + cn f n ( x) = 0 for every x in (a, b), then these
functions are linearly dependent.
A set which is not linearly independent is said to be linearly dependent.
There are two notions of linear independence, according as we allow the coefficients
ck, k = 1, 2, 3, ...., n to assume only real values or also complex values. In the first case, one says
that the functions are linearly independent over the field of reals; in the second case, that they
are linearly independent over the complex field.
Lemma 1.2.1 : A set of real valued functions on an interval (a, b) is linearly independent over
the complex field if and only if it is linearly independent over the real field.
Proof : If the set of real valued functions on an interval (a, b) is linearly independent over
the complex field then it is linearly independent over the field of reals.
Conversely suppose the set is linearly independent over the real field. Therefore for
n
α j ∈ R, ∑ α j f j ( x ) = α1 f1 ( x ) + α 2 f 2 ( x ) + α3 f 3 ( x ) + ⋅⋅⋅ + α n f n ( x ) = 0 for all x in (a, b)
j =1
implies aj = 0 for all j = 1, 2, 3...., n. Let ∑ c j f j ( x ) = 0 for all x in (a, b) and for some
n
j =1
c j ∈ C , j = 1, 2,3 ⋅⋅⋅, n. Since the function f j are real valued and ∑ c j f j ( x) = 0,
c j – c j*
( )
* n n
∑ c j f j ( x) = 0 . implies ∑ c j f j ( x ) = 0 . Thus, ∑
* f j ( x ) = 0 . But c j – c j* / i
j =1 j =1 i
are all real and the set is linearly independent over the real field therefore c j = c j*. But then cj’s
Differential Equations (10)
n
are all real therefore ∑ c j f j ( x) = 0 implies cj = 0 for j = 1, 2,....n.
j =1
A set of functions which is linearly dependent on a given domain may become linearly
independent when the functions are extended to a larger domain. However, a linearly independent
set of functions clearly remain linearly independent on the restricted domain.
Illustration 1 : The functions f1 and f2 define by f1(x) = Cos x and f2(x) = Sin x are linearly
independet on the real line IR and therefore are linearly independent on (0, 2 p).
Illustration 2 : The functions f1 and f2 define by f1(x) = x, f2(x) = | x | are linearly indepent on
the interval (–1, 1) but is not linearly independent on the interval (0, 1) as on the interval
(0, 1), f1(x) = f2(x).
Theorem 1.2.1 :
Let a1, a2 be constants and consider the equation L( y ) = y ′′ + a1 y ′ + a2 y = 0. The two
solutions of L (y) = 0 given in the theorem 1.1.1 are linearly independent on any interval I.
Proof : Let r1, r2 be the roots of characteristic polynomial p(r) = r2 + a1 r + a2.
Case 1.
Differentiation of above equation with respect to x gives c2 (r2 – r1 )e(r2 – r1 ) x = 0 for all x in
I.
Since, r2 ¹ r1 and exponential function in non-zero, c2 is zero. But if c2 is zero then
c1 + c2 e( r2 – r1 ) x = 0 implies c1 is zero. Thus, c1e r1x + c2 e r2 x = 0 implies c1 = c2 = 0.
Therefore φ1 ( x ) = e r1 x and φ2 ( x ) = e r2 x are linearly independent.
Case 2.
Defination 1.5 : Assume that each of the functions f1 ( x ), f 2 ( x), f3 ( x ), ⋅⋅⋅, f n ( x ) are
differentiable atleast (n – 1) times in the interval (a, b). Then the determinant
Differential Equations (11)
f1 ( x ) f 2 ( x) f3 ( x) L fn ( x)
f1′ ( x) f 2′ ( x ) f3′ ( x) L f n′ ( x )
f1′′ ( x) f 2′′ ( x) f3′′ ( x) L f n′′ ( x)
M M M M
f1(n –1) ( x ) f 2(n –1) ( x ) f 3(n –1) ( x ) L f n (n –1) ( x )
denoted by W( f1, f 2, f 3, ...., f n ) ( x) is called the wronskian of the n functions f1, f 2, f 3, ...., f n .
Theorem 1.2.2 :
Two solutions f1, f2 of L (y) = 0 are linearly independent on an interval I if and only if
W(φ1, φ2 ) ( x) ≠ 0 for all x in I.
φ1 ( x0 ) φ 2 ( x0 ) c1 0
c = 0
φ1′ ( x0 ) φ2′ ( x0 ) 2
has a solution c1, c2 where at least one of these numbers is not zero. Let c1, c2, be such a solution
and consider the function ψ ( x ) = c1 φ1 ( x) + c2 φ2 ( x ). Now L(ψ ) = 0 and ψ ( x0 ) = 0, ψ ′( x0 ) = 0.
1
Therefore || ψ ( x0 ) ||= | ψ ( x0 ) |2 + | ψ ′( x0 ) |2 2 = 0. By theorem 1.1.4 || ψ ( x) || = 0. But
|| ψ ( x) ||= | ψ ( x) |2 + | ψ ′( x) |2 = 0. Therefore ψ ( x ) = 0 for all x in I and thus
c1 φ1( x) + c2 φ 2 ( x) = 0 for all x in I. But then f1 and f2 are linearly dependent. Thus, the
supposition W(φ1, φ2 ) ( x0 ) = 0 must be false and therefore W(φ1, φ 2 ) ( x) ≠ 0 for all x in I.
In the next theorem we will prove that we need to compute W(φ1, φ 2 ) at only one point to
test the linear independence of the solutions f1 and f2 .
Differential Equations (12)
Theorem 1.2.3 :
Let f1, f2 be two solution of L(y) = 0 on an interval I and let x0 be any point in I. Then two
solutions f1 and f2 are linearly independent on I if and only if W(φ1, φ2 ) ( x0 ) ≠ 0.
Proof : If f1 and f2 are linearly independent on I then by theorem 1.2.2, W(φ1, φ 2 ) ( x) ≠ 0 for
all x in I. In particular W(φ1, φ2 ) ( x0 ) ≠ 0 conversely, suppose W(φ1, φ2 ) ( x0 ) ≠ 0 and
suppose c1, c2 are constants such that c1 φ1 ( x ) + c2 φ 2 ( x ) = 0 for all x in I. Then
c1 φ1 ( x0 ) + c2 φ 2 ( x0 ) = 0 and c1 φ1′ ( x0 ) + c2 φ 2′ ( x0 ) = 0.
φ1 ( x0 ) φ 2 ( x0 ) c1 0
i.e. c = 0
φ1′ ( x0 ) φ2′ ( x0 ) 2
Theorem 1.2.4 :
Let f1, f2 be any two linearly independent solutions of L(y) = 0 on an interval I. Every
solution f of L(y) = 0 can be written uniquely as
φ = c1 φ1 + c2 φ2 where c1, c2 are constants.
φ1 ( x0 ) φ 2 ( x0 ) c1 α
c = β
φ1′ ( x0 ) φ 2′ ( x0 ) 2
Since W(φ1, φ2 ) ( x0 ) ≠ 0, above system of equations has a unique solution c1, c2 . For
this choice of c1, c2 the functionψ ( x ) = c1 φ1 ( x ) + c2 φ2 ( x ) satisfies ψ ( x0 ) = c1 φ1 ( x0 ) + c2 φ2 ( x0 )
= α = φ ( x0 ) i.e. ψ ( x0 ) = φ ( x0 ) similarly ψ ′( x0 ) = φ ′( x0 ) and L(ψ ) = 0 . From the uniqueness
theorem 1.1.5 it follows that ψ = φ on I i.e. φ = c1 φ1 + c1 φ2 .
Examples :
Q1. Show that the functions ex, e2x, e3x are linearly independent.
Ans. :
Method 1 :
Let c1e x + c2 e 2 x + c3e3 x = 0
then c1 + c2 e x + c3e 2 x = 0 .......... (1)
Differential Equations (13)
Differentiate above equation (1) with respect to x then c2ex + 2 c3 e2x = 0 implies
c2 + 2c3e x = 0 .......... (2)
Let φ1 ( x) = e x , φ2 ( x) = e 2 x , φ3 ( x) = e3x
ex e2 x e3 x 1 1 1
x
W(φ1, φ2 , φ3 ) = e x 2e 2 x 3e3 x = e e 2 x e3 x 1 2 3
ex 4e 2 x 9e 3 x 1 4 9
Again differentiate above equation with respect to x then r 2c2e rx = 0. But r ≠ 0 and e rx ≠ 0
therefore c2 = 0 and from equation (1) we get c1 = 0. Thus f1, f2 are linearly independent.
Method 2 :
x e rx x 1
W(φ1, φ2 ) = = e rx
1 re rx 1 r
Ans. (ii) :
Let c1 φ1 + c2 φ 2 = 0
i.e. c1 x 2 + c2 5 x 2 = 0
if (c1 + 5c2 ) x 2 = 0
Ans. (iii) :
For x > 0 c1 φ1 + c2 φ 2 = (c1 + c2 ) x as | x |= x
and for x < 0 c1 φ1 + c2 φ 2 = (c1 – c2 ) x as | x |= – x
Thus, c1 φ1 + c2 φ 2 = 0 for x ∈ R
⇒ (c1 + c2 ) x = 0 and (c1 – c2 ) x = 0
for every x ∈ R above two equations hold true if and only if c1 = c2 = 0. Thus f1, f2 defined by
φ1 ( x ) = x and φ 2 ( x ) =| x | are linearly independent.
Ans. (iv) :
φ1 ( x) = cos x; φ2 ( x) = sin x
cos x sin x
W(φ1, φ2 ) ( x) = =1
– sin x cos x
Q W (φ1 , φ2 ) ( x) = 1 ≠ 0, φ1 , φ 2 are linearly independent.
Ans. :
The characteristic polynomial p (r ) = r 2 + n2 has roots r1 = in, r2 = –in and therefore the
But φn (0) = cn, φn (2π ) = cn ; φn′ (0) = ndn, φn′ (2π ) = ndn
Q4. (a) : Show that fn (x) = Sin nx satisfies the boundary value problem y² + n2y = 0,
y ( 0 ) = 0 , y ( p ) = 0, n = 1, 2.....
Ans. 4(a) :
Method 1 :
Method 2 :
φn ( x ) = sin nx, φn′ ( x ) = n cos nx
φn′′ ( x ) = – n 2 sin nx
Thus, φn′′ ( x ) + n 2φn ( x ) = – n 2 sin nx + n 2 sin nx = 0 .
Ans. 4(b) :
Working on the similar line as in example 2 we get,
π π
(n2 – m2 ) ∫ φn ( x) φm ( x) dx = (n2 – m2 ) ∫ sin nx sin mx dx
0 0
π
= [sin nx (– m cos mx) – sin mx(– n cos nx)]0
= 0 (as sin 0 = sin np = 0)
π
Since n ≠ m, ∫ φ n ( x) φ m ( x) dx = 0.
0
Q5 : Suppose f1, f2 are linearly independent solutions of the constant coefficient equation
y ′′ + a1 y ′ + a2 y = 0, Let W (f1, f2 ) be abbreviated to W. Show that W is constant if and
only if a1 = 0.
Ans. :
φ1 φ2
W = W(φ1, φ2 ) =
φ1′ φ2′
(
= φ1 φ2′ – φ2 φ1′ )
Then (
W ′ = φ1 φ 2′ – φ 2 φ1′
′
)
= φ1 φ2′′ + φ1′ φ2′ – φ2′ φ1′ – φ2 φ1′′
= φ1 φ2′′ – φ2 φ1′′
But f1 and f2 are solutions of y ′′ + a1 y ′ + a2 y = 0.
= – a1 ( φ1 φ 2′ – φ 2 φ1′ )
= – a1W
Thus, W ′ = 0 iff a1 = 0
Therefore W = constant if and only if a1 = 0
Q6 : Let f1, f2 be two different function on an interval I, which are not necessarily
solutions of an equation L(y) = 0. Prove the following
(a) If f1, f2 are linearly dependent on I then W(f1, f2 ) (x) = 0 for all x in I
(b) If W(f1, f2 ) (x0) ¹ 0 for some x0 in I, then f1, f2 are linearly independent on I.
(c) W(f1, f2 )(x) = 0 for all x in I does not imply that f1, f2 are linearly dependent on I.
(d) W(f1, f2 ) (x) = 0 for all x in I and f2 (x) ¹ 0 on I, imply that are f1, f2 linearly
dependent.
Ans. 6(a) :
Suppose f1, f2 are linearly dependent on I then c1 φ1 ( x ) + c2 φ 2 ( x ) = 0 for some non-zero
c1 and c2.
c2
i.e. φ1 ( x) = – φ2 ( x).
c1
φ1 φ2
W(φ1, φ2 ) ( x) = = φ1 ( x) φ2′ ( x) – φ2 ( x) φ1′ ( x)
φ1′ φ2′
c c
∴ W(φ1 , φ2 ) ( x) = – 2 φ2 ( x ) φ2′ ( x) – φ 2 ( x) – 2 φ2′ ( x ) = 0
c1 c1
∴ W(φ1 , φ2 ) ( x ) = 0 for all x ∈ I.
Ans. 6(b) :
Suppose c1 φ1 ( x ) + c2 φ 2 ( x ) = 0 then
c1 φ1′ ( x) + c2 φ 2′ ( x) = 0
Thus we have a system of equation
φ1 ( x ) φ 2 ( x ) c1 0
=
φ1′ ( x ) φ2′ ( x ) c2 0
φ1 ( x0 ) φ 2 ( x0 ) c1 0
c = 0
φ1′ ( x0 ) φ2′ ( x0 ) 2
Thus, c1 = c2 = 0 if and only if the coefficient matrix is invertible i.e. the determinant of
coefficient matrix is non-zero
φ1 ( x0 ) φ 2 ( x0 )
But = W(φ1 , φ2 ) ( x0 ) ≠ 0
φ ′
1 0( x ) φ ′ (
2 0 x )
Since, W(φ1, φ2 ) ( x0 ) ≠ 0 ⇒ c1 = c2 = 0
∴ c1 φ1 ( x ) + c2 φ 2 ( x ) = 0 ⇒ c1 = c2 = 0.
Hence f1 and f2 are linearly independent on I.
Ans. 6(c) :
Define φ1 ( x) = x 2 , φ2 ( x) = x | x |
for x > 0, | x |= x ∴ φ1 ( x) = x 2 , φ2 ( x) = x 2
x2 x2
∴ W(φ1, φ2 ) = = 0.
2x 2x
for x = 0, φ1 ( x ) = φ 2 ( x ) = 0 ∴ W(φ1 , φ 2 ) = 0
x2 – x2
∴ W(φ1, φ 2 ) = = 0.
2 x –2 x
Let c1 φ1 ( x ) + c2 φ 2 ( x ) = 0
⇒ c1 + c2 = 0 .......... (i)
for x < 0, c1 φ1 ( x ) + c2 φ2 ( x ) = c1 x 2 – c2 x 2 = 0.
⇒ c1 – c2 = 0 .......... (ii)
But c1 + c2 = 0 and c1 – c2 = 0 ⇒ c1 = c2 = 0
Thus, c1 φ1 + c2 φ 2 = 0 ⇒ c1 = c2 = 0
Therefore f1, f2 are linearly independent.
⇒ φ1 ( x ) φ2′ ( x ) – φ2 ( x ) φ1′ ( x ) = 0
⇒ φ2 ( x ) φ1′ ( x ) – φ1 ( x ) φ2′ ( x ) = 0
Since φ 2 ( x ) ≠ 0 ∀x ∈ I
φ2 ( x) φ1′ ( x) – φ1 ( x) φ2′ ( x)
∴ =0
φ2 2 ( x )
φ ′ φ
⇒ 1 = 0 ⇒ 1 = constant = k (say)
φ2 φ2
Q7 : If f1, f2 are two solution of L(y) = 0 on an interval I containing a point x0, then
W(φ1 , φ 2 )( x ) = e – a1 ( x – x0 ) W(φ1 , φ 2 )( x0 ).
Ans. :
Since f1, f2 are solution of L(y) = 0,
φ1′′ + a1 φ1′ + a2φ1 = 0
φ2′′ + a1 φ2′ + a2 φ2 = 0
On multiplying the first equation by –f2, second equation by f1 and adding we obtain
φ1 φ2′′ – φ2 φ1′′ + a1 (φ1 φ2′ – φ2 φ1′ ) + a2 (φ1 φ 2 – φ2 φ1 ) = 0
(φ1 φ2′′ – φ2 φ1′′ ) + a1 (φ1 φ2′ – φ2 φ1′ ) = 0 .......... (i)
φ1 ( x) φ2 ( x)
Let W = W (φ1, φ2 ) ( x ) =
φ1′ ( x ) φ 2′ ( x )
Then W = φ1 ( x) φ 2′ ( x) – φ 2 ( x) φ1′ ( x)
and W ′ = φ1 ( x) φ2′′ ( x) + φ1′ ( x) φ2′ ( x) – φ 2′ ( x) φ1′ ( x) – φ2 ( x) φ1′ ( x)
= φ1 ( x ) φ 2′′ ( x ) – φ2 ( x ) φ1′′ ( x )
Thus, equation (i) becomes
W ′ +a1W = 0.
Thus W satisfies the first order differential equation
W ′ +a1W = 0
= e – a1 ( x – x0 ) W( x0 )
Therefore W(φ1, φ2 ) ( x) = e – a1 ( x – x0 ) W(φ1 , φ 2 ) ( x0 )
EXERCISES
(ii) φ1 ( x ) = sin x, φ2 ( x ) = e ix
(iv) φ1 ( x) = 1, φ2 ( x) = cos x
(vii) φ1 ( x) = cos x, φ2 ( x) = ei x + e – i x
2. State whether the following statements are true or false.
(a) If f1, f2 are linearly independent functions on an interval I, they are linearly independent
on any interval J contained inside I.
(b) If f1, f2 are linearly dependent on an internal I, they are linearly dependent on any internal
J contained inside I.
(c) If f1, f2 are linearly independent solutions of L (y) = 0 on an internal I, they are linearly
independent an any internal J contained inside I.
(d) If f1, f2 are linearly dependent solutions of L (y) = 0 on an interval I, they are linearly
dependent on any internal J contained inside I.
Ans. : 1.
(i) independent (ii) independent (iii) independent
(iv) independent (v) independent (vi) dependent
(vii) dependent.
Ans. : 2.
(a) false (b) true (c) true (d) true
S
Differential Equations (21)
Unit 3 : The Homogeneous Equation of Order n
Everything we have done for the second order equation can be carried over to the case of
the equation of order n. Here, we are concerned with the equation
L( y ) = y (n ) + a1 y (n –1) + a2 y (n –2) + ⋅⋅⋅ + an y = 0,
where, a1, a1, a3 ,......, an are constants.
Theorem 1.3.1 :
Let r 1 , r 2 , r 3 ,......., r s be the distinct roots of the characteristic polynomial
p(r ) = r n + a1r n –1 + a2 r n –2 + ⋅⋅⋅ + an and suppose r i has multiplicity mi (m1 + m2 + m3 + ⋅⋅⋅
+ ms = n). Then n functions
[Polynomial of order n – mi ]
mi –( mi –1)
p ( mi –1) (r ) = (r – ri )
= (r – ri ) [Polynomial of order n – mi ]
Therefore, p( ri ) = p′( ri ) = p′′( ri ) = ⋅⋅⋅ = p ( mi –1) ( ri ) = 0.
Let erx be a solution of L(y) = 0. We see that L(e rx ) = p (r )e rx where p(r ) = r n + a1r n –1
+ a2r n –2 + ⋅⋅⋅ + an.
∂ k rx
∂k
∂r k
L ( e rx
) = L
∂r
k rx
k e =L x e ( )
Differential Equations (22)
k ( k – 1) ( k –2)
= p ( k ) (r ) + kp ( k –1) (r ) x + p (r ) x 2 + ⋅⋅⋅ + p (r ) x k e rx
2!
Thus for r = r i and k = 0, 1, 2,.....m i – 1 we get L( x k e ri x ) = 0. Therefore
x k e ri x , k = 0,1, 2,......mi –1 , are solutions of L(y) = 0. This is true for every characteristic root ri
with multiplicity mi. i.e. x k e ri x , k = 0,1, 2,....mi – 1, i = 1, 2,3,....s are solutions of L(y) = 0 and
the result follows.
Theorem 1.3.2 :
The n solutions of L(y) = 0 given in theorem 1.3.1 are linearly independent on any
interval I.
Proof : We prove that functions given in theorem 1.3.1 satisfy the condition given in defination
1.3.
+.... + (c e + c x e + c x e + .... + c
s0
rs x
s1
rs x
s2
2 rs x
s ( ms –1) x
ms –1 rs x
e ) = 0.
Define pi ( x) = ci 0 + ci1x + ci 2 x 2 + .... + ci (mi –1) x mi –1
Theorem 1.3.3 :
Let f be any solution of
L( y ) = y ( n ) + a1 y ( n –1) + a2 y ( n –2) + .... + an y = 0
on an interval I containing a point x0. Then for all x in I
|| φ ( x0 ) || e – k|x – x0| ≤ || φ ( x) || ≤ || φ ( x0 ) || e k ( x – x0 )
where, k = 1+ | a1 | + | a2 | + | a3 | +....+ | an |
1
and || φ ( x) || = | φ ( x) |2 + | φ ′ ( x) |2 +....+ | φ ( n –1) ( x) |2 2
Proof : This proof is similar to the proof of theorem 1.1.4.
Let u ( x) = || φ ( x) ||2
= | φ |2 + | φ ′ |2 +....+ | φ (n –1) |2
= φ φ + φ ′ φ ′ + .... + φ (n –1)φ (n –1)
Hence u ′( x) = φ ′φ + φ φ ′ + φ ′′ φ ′ + φ ′φ ′′ + .... + φ ( n –1) φ ( n) + φ ( n) φ ( n –1)
Therefore | u′( x) | ≤ 2 | φ ( x) | | φ ′( x) | + 2 | φ ′ | | φ ′′ | +.... + 2 | φ ( n –1) | | φ ( n) |
Since f is solution of L(y) = 0, L(f ) = 0 and therefore
φ ( n ) = – a1 φ ( n –1) – a2 φ ( n –2) – a3 φ ( n –3) – .... – an φ
On substituting the expression for φ (n) we get
Theorem 1.3.5
If φ1, φ2 , φ3 ,....φn , are n solutions of L(y) = 0 on an interval I, they are linearly independent
if and only if W(φ1, φ2 , φ3 ,....φn ) ( x) ≠ 0 for all x in I. (definition 1.5)
Proof : The proof is entirely similar to the proof of theorem 1.2.2
Suppose W(φ1, φ2 , φ3 ,....φn ) ( x) ≠ 0 for all x in I. Let c1, c2, c3,...., cn be constants such
Conversely, assume that φ1, φ2 ,....φ n are linearly independent on I. Suppose there is a
point x0 in I such that W(φ1, φ2 , φ3 ,....φn ) ( x0 ) = 0. Then the system of equations
φ1 ( x0 ) φ2 ( x0 ) φ3 ( x0 ) L φn ( x0 )
c1 0
φ1′ ( x0 ) φ2′ ( x0 ) φ3′ ( x0 ) L φn′ ( x0 )
c2 0
φ1′′ ( x0 ) φ 2′′ ( x0 ) φ3′′ ( x0 ) L φ n′′ ( x0 ) c3 = 0
M M M M
M M
φ1(n –1) ( x0 ) φ2 (n –1) ( x0 ) φ3(n –1) ( x0 ) cn 0
φn (n –1) ( x0 )
has a solution c1 , c2 , c3, ....., cn where at least one of these numbers is not zero. Let c1, c2 ,....., cn
be such a solution and consider a function
ψ ( x ) = c1 φ1 ( x ) + c2 φ2 ( x ) + .... + cn φ n ( x ).
Now L(ψ ) = 0 and ψ ′( x0 ) = ψ ′′( x0 ) = .... = ψ ( n –1) ( x0 ) = 0.
is a solution of L(y) = 0 satisfying the given initial conditions φ (i ) ( x0 ) = αi , i = 0,1, 2,...., n –1.
These constants c1, c2 , c3 ,....., cn would have to sartisfy
φ1 ( x0 ) φ 2 ( x0 ) φ3 ( x0 ) L φn ( x0 )
c1 α1
α
φ1′ ( x0 ) φ 2′ ( x0 ) φ3′ ( x0 ) L φn′ ( x0 )
c2 2
c3 = α3
φ1′′ ( x0 ) φ 2′′ ( x0 ) φ3′′ ( x0 ) L φ n′′ ( x0 )
M M M M M M
cn α n
φ1(n –1) ( x0 ) φ 2(n –1) ( x0 ) φ3(n –1) ( x0 ) φn (n –1) ( x0 )
Since φ1, φ2 , φ3 ,....φ n are linearly independent, by theorem 1.3.5, the determinant of the
coefficients i.e. W(φ1, φ 2 , φ3 ,....φ n ) ( x0 ) ≠ 0. Thus the coefficient matrix is invertible. Therefore
there is a unique set of constants c1, c2 , c3 ,....., cn satisfying above system of equations. For this
choice of c1, c2 , c3 ,....., cn the function
φ ( x ) = c1 φ1 ( x ) + c2 φ 2 ( x ) + c3 φ3 ( x ) + .... + cn φ n ( x )
will be the desired solution.
Theorem 1.3.7 :
Let φ1, φ2 , φ3 ,....φ n be n linearly independent solutions of L(y) = 0 on an interval I. If
c1, c2 , c3 ,....., cn are any constants
φ ( x ) = c1 φ1 ( x ) + c2 φ 2 ( x ) + c3 φ3 ( x ) + .... + cn φ n ( x )
is a solution and every solution may be represented in this form.
Proof : Since φi , i = 1, 2, 3.....n is solution of L(y) = 0, L(φi ) = 0, i = 1, 2, 3.....n .
Theorem 1.3.8
Let φ1, φ2 , φ3 ,....φn be n solutions of L(y) = 0 on an interval I constaining a point x0. Then
φ1 ( x ) φ2 ( x) φ3 ( x ) L φn ( x)
φ1′ ( x) φ 2′ ( x ) φ3′ ( x) L φ n′ ( x )
W ′(φ1, φ2 , φ3 ,..., φn ) ( x) = M M M L M
φ1(n –2) ( x) φ2 (n –2) ( x) φ3(n –2) ( x) φ n (n –2) ( x)
– a1 φ1(n –1) ( x) – a1 φ 2(n –1) ( x) – a1 φ3(n –1) ( x) L – a1 φ n (n –1) ( x)
Since, φ1 ( x ) φ2 ( x) φ3 ( x ) L φn ( x)
φ1′ ( x ) φ 2′ ( x ) φ3′ ( x ) L φ n′ ( x )
M M M L M =0
φ1( n –2) ( x) φ2( n –2) ( x) φ3( n –2) ( x) φ n( n –2) ( x)
– ak φ1(n – k ) ( x ) – ak φ2 (n –k ) ( x ) – ak φ3(n – k ) ( x ) L – ak φ n (n – k ) ( x )
for k = 2, 3, 4,....n, as two rows of the determinant are constant multiplies of each other are
Thus,
φ1 ( x) φ 2 ( x) φ3 ( x) L φ n ( x)
φ1′ ( x ) φ 2′ ( x) φ3′ ( x) L φ n′ ( x )
W ′(φ1, φ2 , φ3 ,..., φn ) ( x) = – a1 φ1′′ ( x ) φ 2′′ ( x) φ3′′ ( x) L φ n′′ ( x )
M M M M
φ1(n –1) ( x ) φ 2 (n –1) ( x) φ3(n –1) ( x) L φ n (n –1) ( x )
= – a1 W(φ1, φ 2 , φ3 ,..., φ n ) ( x)
e a1 x W ( x ) = e a1x0 W ( x0 )
or W ( x) = e – a1 ( x – x0 ) W ( x0 )
Theorem 1.3.9
Let φ1, φ2 , φ3 ,....φn be n solutions of L(y) = 0 on an interval I containing x0. Then they are
linearly independent on I if and only if W(φ1, φ 2 , φ3 ,..., φ n ) ( x0 ) ≠ 0
EXAMPLES
ex x ex e– x sin x cos x
ex (1 + x) e x –e– x cos x – sin x
W (φ1 , φ 2 , φ3 , φ 4 , φ5 ) ( x) = e x (2 + x) e x e– x – sin x – cos x
ex (3 + x) e x –e– x – cos x sin x
ex (4 + x) e x e– x sin x cos x
Differential Equations (30)
1 0 1 0 1
1 1 –1 1 0
W (φ1, φ2 , φ3 , φ4 , φ5 ) (0) = 1 2 1 0 –1
1 3 –1 –1 0
1 4 1 0 1
The row transformations
R2 – R1 , R3 – R1 , R4 – R1 , R5 – R1 gives
1 0 1 0 1
0 1 –2 1 –1
W(φ1, φ2 , φ3 , φ4 , φ5 ) (0) = 0 2 0 0 –2
0 3 –2 –1 –1
0 4 0 0 0
1 –2 1 –1
2 0 0 –2
=
3 –2 –1 –1
4 0 0 0
0 0 –2 2 0 –2 2 0 0
= –2 –1 –1 + 2 3 –1 –1 + 3 –2 –1
0 0 0 4 0 0 4 0 0
2 0 0
+ 3 –2 –1 = – 32
4 0 0
Thus, W(φ1, φ2 , φ3 , φ4 , φ5 ) = e x W(φ1, φ2 , φ3 , φ4 , φ5 ) (0) = –32e x
Ans (c) :
The general solution f is φ ( x) = c1e x + c2 x e x + c3 e – x + c4 sin x + c5 cos x
The initial conditions φ (0) = 1, φ ′ (0) = φ ′′(0) = φ ′′′(0) = φ (iv) (0) = 0 gives the following
system of equations.
1 0 +1 1 0 c1 1
1 1 –1 1 0 c 0
2
1 2 1 0 –1 c3 = 0
1 3 –1 –1 0 c4 0
1 4 1 0 1 c5 0
Ans. (a) :
p( r ) = r 4 – (2 i ) 4 = ( r 2 + (2 i ) 2 ) (r 2 – (2 i ) 2 )
= (r 2 – i 2 (2 i )2 ) (r 2 – ( i 2) 2 )
= ( r + 2i i ) ( r – 2i i ) ( r + 2 i ) ( r – 2 i )
Thus, p ( r ) = ( r + 2i i ) ( r – 2i i ) ( r + 2 i ) ( r – 2 i )
π
π π i
i = cos + i sin = e 2
2 2
1
iπ
2 π
π π
\
i
i = e 2
= e 4 = cos + i sin
4 4
1+ i i(1 + i) –1 + i
Therefore i= , i i= =
2 2 2
The roots of characteristic polynomial are – 2(–1 + i), 2(–1 + i ), 2(1 + i ), – 2(1 + i)
Thus four linearly independent solutions are
2 + 2 i) x 2 +i 2 ) x
e( 2 –i 2 ) x
, e(– , e( , e(– 2 –i 2 ) x
3 2
Ans. (c) : The characteristic polynomial is p (r ) = r – 5r + 6r and its roots are 0, 3, 2. Thus
three linearly independent solutions are given by 1, e3x, e2x and any solution f has the
form φ ( x) = c1 e3 x + c2 e2 x + c3
Ans. (a) : The characteristic polynomial p (r ) = r 3 – 4r and its roots are 0, 2, –2. Thus, three
linearly independent solution are e° = 1, e 2 x , e –2 x and every solution f has the form
φ ( x) = c1 + c2e 2 x + c3e –2 x
Ans. (b) :
W(φ1, φ2 , φ3 ) ( x) = e0( x –0) W(φ1, φ 2 , φ3 ) (0)
1 e2x e –2 x
W(φ1, φ2 , φ3 ) ( x) = 0 2e 2 x –2e –2 x
0 4e 2 x 4e –2 x
1 1 1
W(φ1, φ2 , φ3 ) (0) = 0 2 –2
0 4 4
Thus, W(φ1, φ2 , φ3 ) ( x) = 16 .
Ans. (c) :
φ (0) = 0, φ ′ (0) = 1, φ ′′ (0) = 0,
φ ( x) = c1 + c2e 2 x + c3e –2 x , φ (0) = c1 + c2 + c3 = 0 and so on
1 1 1 c1 0
0 2 –2 c = 1
2
0 4 4 c3 0
R3 – 2 R2 gives
1 1 1 c1 0
0 2 –2 c = 1
2
0 0 8 c3 –2
Thus, φ ( x ) = c1 + c2 e 2 x + c3e –2 x =
4
(
1 2x
)
e – e –2 x is the required solution.
EXERCISE
(c) φ1 ( x ) = x, φ2 ( x ) = e 2 x , φ3 ( x) = | x |
Ans. 1 :
(a) True (b) false
Ans. 2 :
(a) independent (b) dependent (iii) independent
Ans. 3 :
(a) φ1 ( x ) = e –4 x , φ2 ( x ) = e – x
(b) φ1 ( x) = e –2 x , φ2 ( x) = xe –2 x , φ3 ( x) = x 2 e –2 x
Theorem 1.4.1
Let b(x) be continuous on an interval I. Every solution y of L(y) = b (x) on I can be
written as ψ = ψ p + c1 φ1 + c2 φ2 .
Where yp is a particular solution, f1, f2 are two linearly independent solutions of L(y) = 0
and c1, c2 are constants. A particular solution yp is given by
x
[φ1 (t ) φ 2 ( x) – φ1 ( x) φ 2 (t )] b (t )
ψ p ( x) = ∫ dt.
x0 W(φ1 , φ2 ) (t )
Conversely every such y is a solutions of L(y) = b (x)
Proof :
Let y and yp be two solutions of
L( y ) = y′′ + a1 y ′ + a2 y = b
φ1′u1′ + φ 2′u 2′ = b( x )
Then u1φ1 + u2φ2 will satisfy L(y) = b(x).
x φ (t ) b(t )
φ ( x ) = φ1 ( x ) – ∫ 2
dt + u1 ( x0 )
x0 W(φ1 , φ2 ) (t )
x φ (t ) b(t )
+ φ2 ( x) + ∫ 1
dt + u2 ( x0 )
x0 W(φ1, φ2 ) (t )
x
[φ1 (t ) φ 2 ( x ) – φ 2 (t ) φ1( x)] b(t )
ψ p ( x) = ∫ dt
x0 W(φ1, φ 2 ) (t )
The function yp(x) is a solution of L(y) = b (x).
Theorem 1.4.1 provides a method to find a solution of second order non-homogeneous
differential equation with constant coefficients. The same procedure can be generalized for the
non-homogeneous equation of order n.
ψ p (n ) = ∑ ui′φi (n –1) + ∑ ui φi (n )
EXAMPLES
cos x sin x e– x
W(φ1, φ2 , φ3 ) ( x) = – sin x cos x – e – x
– cos x – sin x e– x
W(φ1, φ2 , φ3 ) ( x) = e – a1 ( x – x0 ) W(φ1, φ2 , φ3 ) ( x0 )
= e – x W(φ1 , φ 2 , φ3 ) (0)
1 0 1
W(φ1, φ2 , φ3 ) (0) = 0 1 –1
–1 0 1
= 1[1 – 0] + 1 [0 + 1] = 2
Thus W(φ1 , φ2 , φ3 ) ( x) = 2e – x
cos x 0 e– x
W2 ( x) = – sin x 0 – e – x = e – x (cos x – sin x)
– cos x 1 e– x
cos x sin x 0
W3 ( x) = – sin x cos x 0 = 1
– cos x – sin x 1
x x
W1 (t ) b(t ) e – t (cos t + sin t )
u1 ( x ) = ∫ dt = ∫ – dt
W(φ1, φ2 , φ3 ) (t ) 2e – t
1
= – [ + sin x – cos x]
2
1
Thus, u1 ( x ) = [cos x – sin x]
2
x x –t
W2 (t ) b(t ) e (cos t – sin t )
u2 ( x ) = ∫ dt = ∫ dt
W(φ1, φ2 , φ3 ) (t ) 2e – t
1
= [+ sin x + cos x ]
2
x x
W3 (t ) b(t ) dt dt 1
u3 ( x) = ∫ =∫ – t = e x
W(φ1, φ2 ,φ3 )(t ) 2e 2
Therefore a particular solution is given by
ψ p = u1 ( x ) φ1 ( x ) + u2 ( x ) φ2 ( x ) + u3 ( x ) φ3 ( x)
1 1 1
= – (cos x – sin x) cos x + (cos x + sin x )sin x + e x e – x = 1
2 2 2
The most general solution is
ψ ( x ) = ψ p + c1 φ1 + c2 φ 2 + c3 φ3
= 1 + c1 cos x + c2 sin x + c3 e – x
ψ (0) = 0 ⇒ 1 + c1 + c3 = 0
ψ ′(0) = 1 ⇒ ψ ′( x ) = – c1 sin x + c2 cos x – c3 e – x
Thus, y¢ (0) = c2 – c3 = 1
Ans. (a) : The characteristic polynimial p (r ) = r 3 – r has roots 0, 1, –1 and the linearly
independent solution of the related homogeneous equation are φ1 ( x) = 1, φ2 ( x) = e x ,
φ3 ( x ) = e – x
Let ψ p = u1 ( x ) φ1 ( x ) + u2 ( x ) φ2 ( x ) + u3 ( x ) φ3 ( x)
x
Wk (t ) b (t ) ds
uk ( x ) = ∫ k = 1, 2,3.
W(φ1, φ2 ,φ3 ) (t )
1 ex e– x
W(φ1, φ2 , φ3 ) ( x) = 0 e x –e – x = 2
0 ex e– x
0 ex e– x
b ( x)W1 ( x ) = 0 e x – e – x = –2 x
x ex e– x
1 0 e– x
b ( x)W2 ( x) = 0 0 – e x = xe – x
0 x e–x
1 ex 0
b ( x )W3 ( x) = 0 e x 0 = xe x
0 ex x
ex e2 x
where, W1 (φ1, φ2 ) ( x) = = e3 x ,
x 2x
e 2e
0 e2x
b ( x)W1 ( x) = = – e 2 x sin e – x ,
–x 2x
sin e 2e
ex 0
b ( x)W2 ( x ) = = e x sin e – x .
x –x
e sin e
W1 ( x) b( x) dx W2 ( x) b( x ) dx
u1 ( x) = ∫ , u2 ( x ) = ∫
W(φ1, φ2 )( x) W(φ1 , φ2 ) ( x)
–2 x e –2 x xe –2 x
where, b ( x) = 3 xe , W1 (φ1, φ2 ) ( x) = = e –4 x ,
–2 x –2 x
–2e (1 – 2 x )e
0 xe –2 x
b ( x)W1 ( x) = = –3 x 2 e –4 x ,
–2 x –2 x
3xe (1– 2 x )e
–e 2 x 0
b ( x)W2 ( x) = = 3 xe –4 x .
–2 x –2 x
–2e 3 xe
3x 2e –4 x
Thus, u1 ( x) = ∫ – –4 x
dx = – x3 and
e
3x e –4 x 3 2
u2 ( x ) = ∫ + –4 x
dx = x
e 2
3 3 –2 x
Therefore ψ p = u1 φ1 + u2 φ 2 = – x 3e –2 x + xe
2
1 3 –2 x
= x e
2
The general solution
1 3 –2 x
ψ = c1e –2 x + c2 xe –2 x + xe
2
Q.3. Find the general solution of
π π
y ′′ + y = tan x , – < x< .
2 2
2
Ans : The characteristic polynimial p( r ) = r + 1 has roots +i, –i and the two linearly independent
solutions are
φ1 ( x ) = cos x and φ2 ( x ) = sin x
u2 ( x ) = ∫ sin x d x = – cos x
π
and ψ p = – cos x log (sec x + tan x ), – < x<π2
2
The general solution
Note. The formula for a particular a solution yp of L(y) = b(x) makes sense for some
discontinuous functions b (x). Then yp will be a solution of L(y) = b(x) at the continuity
points of b .
Q 4. Find a particular solution of the equation.
y ′′ + y = b ( x ),
Where, b( x ) = –1 (– Π ≤ x < 0) ,
=1 (0 ≤ x ≤ Π ) ,
=0 (1x1 > Π ).
Ans : Let us find out the particular solution of y ′′ + y = α where a is a constant.
The characteristic polynomial is p( r ) = r 2 + 1 and has roots +i, –i. Therefore the basic
solutions (linearly independent solutions) are
φ1 ( x) = cos x, φ2 ( x) = sin x
Let ψ p = u1 ( x) φ1 ( x) + u2 φ2 ( x) be a particular solution of the equation y ′′ + y = α . By
method of separation of parameters we get,
W1 ( x) b( x) dx W2 ( x) b( x ) dx
u1 ( x) = ∫ , u2 ( x ) = ∫
W(φ1, φ2 )( x) W(φ1 , φ2 ) ( x)
cos x sin x
where, W(φ1, φ2 ) ( x) = =1
– sin x cos x
0 sin x cos x 0
W1 ( x ) = = – sin x, W2 ( x ) = = cos x
1 cos x – sin x 1
W1 ( x) b( x) dx –α sin x
Then, u1 ( x ) = ∫ =∫ dx = α cos x
W(φ1 ,φ 2 ) ( x) 1
c3 = c1* + 1
c5 = c3 – 2 = c1* + 1 – 2 = c1* –1
c7 = c5 – 1 = c1* – 1– 1 = c1* – 2
Thus, the particular solution becomes
ψ ( x ) = 2 cos x –∞ < x < π
= cos x – 1 –π ≤ x < 0
= – cos x + 1 0 ≤ x ≤π
= –2 cos x π <x<∞
If we choose c3 = c1* + 2 then c1 = c3 + 1 = c1* + 3,
c3 = c1* + 2,
c5 = c3 – 2 = c1* ,
c7 = c5 – 1 = c1* –1,
and the particular solution becomes
ψ ( x ) = 3cos x ; – ∞ < x < –π
= 2cos x –1 ; –π ≤ x < 0
=1 ; 0 ≤ x ≤π
= – cos x ; π < x < ∞.
Differential Equations (47)
Thus, we can generate infinitely many particular solutions that are piecewise continuous.
Method of undetermined coefficients :
The method described sofar is called the method of variation of parameters. Although this
method yields a solution of the non-homogeneous equation it sometimes require more labor
than necessary. We now explain a method which is often faster than a method of variation of
parameters. This method is useful to solve the non-homogeneous equation L (y) = b (x), when
b(x) is a solution of some homogeneous equation with constant coefficients. The procedure we
are about to describe is called the method of undetermined coefficients.
For the given different equation L (y) = b (x), suppose b(x) is a solution of some
homogeneous equation M (y) = 0 with constant coefficients. Then M (b(x)) = 0. If y is a solution
of L (y) = b (x) and M (b) = 0 then
M [ L(ψ )] = M (b) = 0.
Therefore y is a solution of the homogeneous equation M (L(y)) = 0 with constant
coefficients. If the order of differential operator Lis n and that of M is m then M(L(y)) = 0 is a
homogeneous differential equation of order m + n and therefore there are m + n linearly
independent solutions of M (L(x)) = 0. Since b(x) is a particular solution of M (y) = 0 every
linear combination of these n + m linearly independent solution will not be a solution of
L (y) = b (x). Thus, to fine the solution of L (y) = b (x) we substitute the linear combination of
solutions into L (y) = b (x) and determine the set of coefficients other than the coefficients of the
solutions corresponding to the homogeneous equation L (y) = 0.
We give an example to show the usefulness of this method. Suppose we consider
L( y) = y′′ – 3 y′ + 2 y = x 2
Since ( x 2 )′′′ = 0, x2 is a solution of M ( y ) = y′′′ = 0.
Every solution y of L (y) = x2 is a solution of
M ( L( y)) = M ( y′′′ – 3 y′ + 2 y) = M ( x 2 ) = 0.
But M ( y′′′ – 3 y ′ + 2 y ) = ( y′′′ – 3 y ′ + 2 y )′′′ = 0
i.e. y (v) – 3 y (iv) + 2 y′′′ = 0.
5 4 3
The characteristic polynomial of this equation is p(r ) = r – 3 y + 2r (just the product
of characteristics polynomials of L and M ). The roots of p (r) are 0, 0, 0, 1, 2 and hence y must
have the form ψ = c0 + c1 x + c2 x 2 + c3e x + c4e 2 x observe that c3e x + c4e 2 x is a solution
of L (y) = 0.
ψ p = c0 + c1 x + c2 x 2
Since y p is a solution, it should satisfy the differential equation L (y) = x2.
L(ψ p ) = ψ p′′ – 3ψ p′ + 2ψ p = 2 c0 + c1x + c2 x 2 – 3[c1 + 2c2 x] + [2c2 ]
EXAMPLES
Exp. 1. Using the annihilator method find a particular solution of each of the following
equations.
(a) y ′′ + 4 y = cos x (b) y′′ – 4 y = 3e2 x + 4e – x
(c) y ′′ – y ′ – 2 y = x 2 + cos x
= y (iv ) + 5 y′′ + 4 y
ψ p′ = [2c1 + (1 + 2 x) c2 ]e 2 x – c3e – x
But yp satisfies y ′′ – 4 y = 3e 2 x + 4e – x
Therefore M [ L( y )] = [ y ′′ – y ′ – 2 y ]( ) + [ y ′′ – y ′ – 2 y ]′′′
v
The problem is to determine the constants c0, c1, c2, c3, c4 so that L(ψ p ) = x 2 + cos x.
L(ψ p ) = ψ p′′ – ψ p′ – 2ψ p
Exp. 1. Use the method of variation of parameters and find the general solution of each of
the following equation.
(i) 6 y ′′ + 5 y ′ – 6 y = x
Exp. 2. Find the particular solution of each of the following equation using the method of
undetermined coefficients.
x –x 1 2 2
Ans.(1): (a) y = c1e + c2e – sin x – (b) c1 cos x + c2 sin x – x 2 cos x + x sin x
5 5
x 4e – x
(c) c1e –2 x + c2 c – x + 2e x (d) c1e – x + c2 xe – x +
12
1 1
(e) c1 cos 2 x + c2 sin 2 x + cos x (f) c1 cos 3x + c2 sin 3 x – x cos 3 x
3 6
x x
6x 1x – 1
(g) c1e + c2 e + (7 cos x + 5sin x ) (b) c1e 2 + c2 e 2 + ex
74 3
2x 3x
– 1 5
(i) c1e 3 + c2 e 2 – x– .
6 36
2 3 2 x 3x 3 3x 3 9
Ans.(2): (a) x – 2 x + + e (b) xe – e + sin x + cos x
2 3 2 10 10
1 3 –2 x
(c) x e (e) e – x
2
S
Differential Equations (52)
Chapter
2
Linear Equations with
Variable Coefficients
Contents :
Unit 1 : Homogenous equations with variable coefficients.
(a) Initial value problems for the homogeneous equation.
(b) Solutions of homogenous equation
(c) Reduction of an order of a homogeneous equation
Unit 2 : Basis
(a) Linear independence and Wronskian
(b) Solution of non-homogeneous equations
Unit 3 : Homogenous equations with analytic coefficients.
Introduction
Solutions to linear equations with variables coefficients are necessary to analysis most of
the situations in science and technology. In the last chapter we have studied linear equations
with constant coefficients. In this chapter we are going to study linear equations with variable
coefficients. There is no standard procedure to find all possible solutions of a given equation.
However it is possible to construct series solution if the coefficient functions and the control
function are analytic on some open set.
Theorem 2.1.1 :
Let b1, b2 , b3 , ...., bn be non-negative constants such that for all x in I
| ai ( x) | ≤ bi i = 1, 2, 3, ...., n and define k by
k = 1 + b1 + b2 + b3 + .... + bn .
If x0 is a point in I and f is a solution of L(y) = 0 on I then
|| φ ( x0 ) || e – k |x – x0 | ≤ || φ ( x) || ≤ || φ ( x0 ) || e k |x – x0 | for all x in I.
Proof : The proof of this theorem is similar to the proof of theorem 1.3.3.
Let u ( x) = || φ ( x) ||2 =| φ |2 + | φ ′( x) |2 + | φ ′′ |2 +....+ | φ (n –1) |2
= φ φ + φ ′φ ′ + φ ′′ φ ′′ + .... + φ ( n)φ ( n–1) + φ ( n –1)φ ( n –1)
Hence u ′( x ) = φ ′ φ + φ φ ′ + φ ′′φ ′ + φ ′ φ ′′ + .... + φ ( n)φ ( n –1) + φ ( n –1)φ ( n)
Theorem 2.1.4
There exist n linearly independent solutions (definition 1.3) of L(y) = 0 on I.
Proof : Let x0 be a point in I. According to theorem 2.1.3 and theorem 2.1.2, there is a unique
solution of L(y) = 0 satisfying given initial conditions at x0.
Let f 1 be a solution of L(y) = 0 satisfying
φ1 ( x0 ) = 0, φ1′ ( x0 ) = 0, φ1′′ ( x0 ) = 0,......, φ1( n –1) ( x0 ) = 0
Let f 2 be a solution of L(y) = 0 satisfying
φ2 ( x0 ) = 0, φ2′ ( x0 ) = 1, φ2′′ ( x0 ) = 0,......, φ2 ( n –1) ( x0 ) = 0
In general Let f i be a solution of L(y) = 0 with
We will prove that these solutions φ1, φ2 , φ3 ,...., φn are linearly independent on I. Suppose
there are constants c1, c2, c3,..... cn such that
c1 φ1 ( x ) + c2 φ 2 ( x) + c3 φ3 ( x) + ..... + cn φ n ( x ) = 0 for all x in I.
Differentiating above equation (n – 1) times we get,
c1 φ1′ ( x) + c2 φ 2′ ( x) + c3 φ3′ ( x) + ..... + cn φ n′ ( x) = 0
Theorem 2.1.5
If f1(x) is a solution of L( y ) = y ′′ + a1 ( x) y′ + a2 ( x) y = 0 on an interval I and f1(x) ¹ 0 on
I, the second solution f2 (x) is given by
x
1 s
φ2 ( x) = φ1 ( x) ∫ 2
exp – ∫ 1
a ( t ) dt ds.
[
x0 1 φ ( s )]
x0
The functions f1 and f2 are linearly independent.
Proof : Since f1 is a solution of L(y) = 0, L(f1) = 0.
Let φ2 ( x ) = u ( x ) φ1 ( x ) be second solution of L(y) = 0.
L(φ 2 ) = L(u φ1 ) = (u φ1 )′′ + a1 ( x) (u φ1 )′ + a2 ( x) (u φ1) = 0.
i.e. u ′′( x ) φ1 ( x ) + 2u ′( x ) φ1′ ( x ) + u ( x ) φ1′′ ( x )
Differential Equations (57)
+ a1 ( x ) u′ ( x) φ1 ( x ) + u ( x) φ1′ ( x) + a2 ( x ) (u ( x )φ1 ( x)) = 0.
Since L(φ1 ) = φ1′′ + a1 ( x ) φ1′ + a2 ( x )φ1 = 0,
φ1 ( x )v′( x ) + 2φ1′ ( x) + a1 ( x ) φ1 ( x ) v ( x ) = 0
v′( x) 2φ1′ ( x)
+ + a1 ( x) = 0
v ( x ) φ1 ( x )
On integrating above equation between the limits x0 to x, we get
x 2φ ′ (t )
log v ( x) – log v ( x0 ) + ∫ φ (t ) + a1 (t ) dt = 0
1
1
x0
x
log v( x) – log v ( x0 ) + 2 [log φ1 ( x ) – log φ1 ( x0 )] + ∫ a1 (t ) dt = 0.
x0
v ( x )φ12 ( x ) x
log = – ∫ a1 (t ) dt
v ( x0 )φ12 ( x0 ) x0
x
v( x)φ12 ( x)
– ∫ a1 ( t ) dt x
=e x0
= exp – ∫ a1 (t ) dt
v( x0 )φ12 ( x0 ) x0
v ( x0 )φ12 ( x0 ) x
i.e. v( x) = exp – ∫ a1 (t ) dt
φ12 ( x) x0
But v ( x0 ) φ12 ( x0 ) are the values of v ( x ) φ12 ( x ) evaluated at point x0 and therefore is constant
c x
v ( x) = exp – ∫ a1 (t ) dt
φ12 ( x ) x0
But v ( x ) = u′( x ) and therefore
x
c s
u ( x) = ∫ 2 exp – ∫ a1 (t ) dt ds.
φ1 ( s ) x0
Since, φ2 ( x ) = u ( x ) φ1 ( x ) we get the required result.
We can generalize above theorem for linear differential equation
Differential Equations (58)
L( y ) = y (n ) + a1 ( x ) y (n –1) + a2 ( x ) y (n –2) + ..... + an ( x) y = 0.
Theorem 2.1.6 :
Let φ1 be a solution of L(y) = 0 on an interval I and suppose φ1 ( x) ≠ 0 on I. Then we can
reduce the order of equation L(y) = 0 by one. If v2, v3,...., vn are linearly independent solutions of
the reduced differential equation of order n – 1 and if vk = u′k , k = 1, 2, 3,..., n, then
φ1, u1φ1, u2 φ 2 ,...., un φ n are linearly independent solutions of L(y) = 0 on I.
Proof : Let φ1 be solution of L(y) = 0 on I. we try to find a solution f of L(y) = 0 of the form
φ = u ( x ) φ1 ( x ), where u(x) is n times differentiable function defined on an interval I. If
φ ( x ) = u ( x) φ1 ( x) is a solution of L(y) = 0 then L(u ( x ) φ1 ( x )) = 0 .
φ1v ( n –1) + n φ1′ + a1 ( x ) φ1 v ( n –2) + .... + n φ1( n –1) + (n –1) a1 ( x ) φ1( n –2) + .... + an –1 φ1 v = 0
Since, φ1 ( x) ≠ 0 on I we can divide above equation byf1. Thus, we can reduce the order of
differential equation by one. Suppose v2, v3, v4,..., vn are linearly independent solutions of the
differential equation in v of order n – 1. Then
φ1vk ( n –1) + nφ1′ + a1 φ1 vk ( n –2) + .... + n φ1( n –1) + ( n – 1) a1 φ1( n –2) + .... + an –1φ1 vk = 0
But then vk ( x ) = uk ′ ( x ) for k = 2,3, 4,...., n
x
and uk ( x ) = ∫ vk (t ) dt k = 2, 3, 4,...., n
x0
EXAMPLES
y ′ – 2 y = 0 Since f is a solution
1 1
Let φ ( x ) = x be a solution to L( y ) = y ′′ +
r
x x
L(y) = 0.
Therefore r ( r – 1) x r –2 + r x r –2 – x r –2 = 0
that is (r 2 –1) x r –2 = 0 for x>0
Thus, r 2 – 1 = 0 or r = +1, – 1.
1
Therefore φ1 ( x) = x and φ2 ( x) = are two solutions of L(y) = 0.
x
Ans (b) :
c2
Let c1φ1 + c2φ 2 = 0 then c1x + = 0 . Differentiate this equation twice with respect to
x
Differential Equations (60)
= 0 implies c2 = 0 and therefore c1 = 0. Thus, f1 and f2 are linearly independent.
2c2
x we get
x3
Ans (c) :
1
φ1 ( x) = c1 x + c2
x
φ1 (1) = 1 and φ1′ (1) = 0 gives
1
c1x + c2 = 1 at x = 1 i.e. c1 + c2 = 1
x
c2
c1 x – = 0 at x = 1 i.e. c1 – c 2 = 0
x2
and φ1 ( x) = x +
1 1 1
Thus, c1 = c2 =
2 2 x
1
Let, φ2 ( x ) = d1x + d 2
x
φ2 (1) = 0 and φ2′ (1) = 1 gives
1 1
d1 + d 2 = 0 and d1 – d 2 = 1. Then d1 = and d2 = –
2 2
φ 2 ( x ) = x – .
1 1
and
2 x
Q. 2. Find two linearly independent solutions of the equation
1
(3 x – 1)2 y′′ + (9 x – 3) y′ – 9 y = 0 for x>
2
dy dy dt
Ans. : Put t = 3 x – 1 then = ⋅ = y& ⋅ 3
dx dt dx
where . represents derivative with respect to t.
d2y d dy d d dt
2
= = y& ⋅ 3 = (3 y& ) = 3&&y ⋅ 3 = 9 &&
y
dx dx dx dx dt dx
Therefore 9t 2 &&
y + 9t y& – 9 y = 0
or t 2 &&
y + ty& – y = 0
Let y = tr be a solution then
r ( r – 1)t r + rt + r – t r = 0 implies (r 2 –1)t r = 0.
But t > 0 therefore r = +1, –1
1
and φ1 (t ) = t and φ2 (t ) = are solutions
t
But t = 3x – 1 and therefore
1
φ1 ( x ) = 3 x –1 and φ2 ( x) = are two solutions of given equation.
3x –1
x2
But then u = k Let k = 2 Then φ ( x ) = u ( x ) φ1 ( x ) = x 2 ( x 3 ) = x 5 is the
2
second solution independent of f1 as c1x 3 + c2 x 5 = 0 implies c1 = c2 = 0.
1 1
Ans (b) : φ1 ( x) = e x = φ1′ ( x) = φ1′′ ( x) Let L( y ) = y ′′ – (1 + ) y ′ + y = 0 . (We can divide the
x x
given equation by x as x > 0.)
1 1
L(φ1 ) = e x – (1 + ) e x + e x = 0 ∴ φ1 is a solution.
x x
To determine second solution, let φ ( x ) = u ( x) φ1 ( x) be a solution then by theorem 2.1.5
1
exp – ∫ a1 (t ) dt dx
x
u ( x) = ∫
φ12 ( x )
= ∫ e –2 x xe x dx = ∫ xe – x dx = – (1 + x) e–x
Thus, f2 (x) = – (1 + x) is a second solution of the equation L(y) = 0
c1φ1 ( x) + c2φ 2 ( x) = c1e x + c2 (–1) (1 + x) = 0 ⇒ c1 + c2 (–1) (1 + x)e – x = 0
But then c2 = c1 = 0 therefore f1 and f2 are linearly independent solutions.
2x 2y
Ans (c) : φ1 ( x) = x, φ1′ ( x) = 1, φ1′′ ( x) = 0. Let L( y ) = y ′′ – 2
y′ + = 0.
1– x 1– x 2
By Theorem 2.1.5.
1
exp – ∫ a1 (t ) dt dx
x
u ( x) = ∫
φ12
1 x –2t dt 1 1
=∫ 2
exp – ∫ 2
dx = ∫ 2 dx
x 1– t x 1– x 2
dx dx 1 dx 1 dx
=∫ =∫ 2 + ∫ + ∫
2 2
x (1– x ) x 2 1– x 2 1 + x
1+ x
= – x –1 + log
1
.
2 1– x
x 1+ x
Then φ2 ( x ) = φ1 ( x ) u ( x ) = –1 + log is a second solution.
2 1– x
Q. 4. One solution of x 3 y′′′ – 3x 2 y′′ – 6 xy′ – 6 y = 0 for x > 0 is f1 (x) = x find the remaining
two independent solutions for x > 0.
Ans : Let f = xu be a solution of L (y) = 0. Then φ ′ = xu ′ + u , φ ′′ = xu ′′ + 2u ′, φ ′′′ = xu ′′′ + 3u ′′ .
φ ′ = u ′φ1 + u φ1′ , φ ′′ = u ′′φ1 + 2u ′φ1′ + uφ1′′ , φ ′′′ = u ′′′φ1 + 3u ′′φ1′ + +3u ′φ1′′ + uφ1′′′
+ a 2 ( x ) u ′φ1 + uφ1′ + a3 ( x ) u φ1 = 0.
Since f1 is a solution φ1′′′ + a1φ1′′ + a2φ1′ + a3φ1 = 0.
φ ′ φ ′ φ φ ′ φ ′′ φ ′ φ ′ φ ′ φ ′ φ φ ′′ 2φ φ ′2
v = 2 = 2 – 2 21 , v′ = 2 – 1 22 – 1 22 – 2 21 + 2 31
φ1 φ1 φ1 φ1 φ1 φ1 φ1 φ1
( φ
) (
= φ2′′′ + a1φ2′′ + a2φ2′ – 2 φ1′′′ + a1φ1′′ + a2φ1′
φ1
)
Since f 2 and f 1 are solutions, φ2′′′ + a1φ2′′ + a2φ2′ + a3φ2 = 0 and φ1′′′ + a1φ1′′ + a2φ1′
+ a3φ1 = 0 and therefore
φ2
L (v ) = – a3 φ 2 – ( – a3φ1 ) = 0.
φ1
φ2 ′
Thus, v = is a solution of reduced equation.
φ1
Differential Equations (64)
EXERCISE
Use the reduction of order method and find the general solution of each of the following
equations. Verify that f1 satisfies the equation.
2 ′ 2
(b) y′′ – y + 2 y = 0, φ1 = x (Ans. y = c1x + c2 x 2 )
x x
x3 x 2
– +
e 3 2
(d) y ′′ + ( x 2 – x ) y ′ – ( x – 1) y = 0, φ1 = x (Ans. y = c1x + c2 x ∫ dx )
x2
– x2
x 1
(e) y ′′ + – y ′ – y = 0, φ1 = x
2 e 4
2 x (Ans. y = c1 x 2 + c2 x 2 ∫ dx)
x3
1
(f) 2 x y ′′ + 3xy′ – y = 0, φ1 = x
2 12
(Ans. y = c1 x2 + c2 x –1 )
Unit 2 : Basis
In the course on linear algebra we learn about a vector space also called as linear space
and the basis of a linear space. Suppose S is a set of functions with the following property.
If f1, f2 Î S, c1f1 + c2f2 Î S for any two constants c1, c2. Then the set S is called a linear
space of functions. If a linear space of functions S contains n functions φ1, φ2 , φ3 ,...., φn which
are linearly independent and every function from S can be represented as a linear combination
of these functions φ1, φ2 , φ3 ,...., φn then the set { φ1, φ2 , φ3 ,...., φn } is called a basis for the linear
space S. The number n is called dimension of S.
For a given linear differential equation L( y ) = y ( n ) + a1 ( x ) y ( n –1) + a2 ( x) y ( n –1) + ⋅⋅⋅
+ an ( x ) y = 0, the collection of all solutions denoted by S of L(y) = 0 is a linear space. Every
basis of S contains n linearly independent functions and therefore dimension of solution space S
is n.
To check the linear independence of functions φ1, φ2 , φ3 ,...., φn , we consider the wronskian
W( φ1, φ2 , φ3 ,...., φn ). There is a relation between the linear indepence of functions and the
Wronskian W( φ1, φ2 , φ3 ,...., φn ). In chapter I we have proved this result for the linear differential
equation with constant coefficients.
Differential Equations (65)
A. Linear Independence and Wronskian
In section 1(B) we have seen that for the differential equation L( y ) = y ( n ) + a1 ( x) y ( n–1)
+ a2 ( x ) y ( n –2) + ⋅⋅⋅ + an y = 0 there are n linearly independent solutions φ1, φ2 , φ3 ,...., φn satisfying
the initial conditions φi(i –1) ( x0 ) = 1, φi( j –1) ( x0 ) = 0, j ≠ i. These linearly independent solutions
is a basis of solution space of L(y) = 0. Every solution of L(y) = 0 can be represented as a linear
combination of these functions φ1, φ2 , φ3 ,...., φn .
Theorem 2.2.1
Let φ1, φ2 , φ3 ,...., φn be n solutions of L(y) = 0 on I satisfying the initial conditions.
Since φ1, φ2 ,...., φ n are solutions of L(y) = 0, by superposition principle (chapter 2 unit
1(B)) y is also a solution of L(y) = 0 and clearly
ψ ( x0 ) = α1 φ1 ( x0 ) + α 2 φ2 ( x0 ) + α3 φ3 ( x0 ) + ⋅⋅⋅ + α n φn ( x0 ) = α1
as φ1 ( x0 ) = 1 and φi ( x0 ) = 0 for i = 2,3, 4,..., n.
Recall that the Wronkian of n functions φ1, φ2 , φ3 ,...., φn is defined as the determinant
φ1 φ2 φ3 L φn
φ1′ φ 2′ φ3′ L φ n′
W (φ1, φ 2 , φ3 ,...., φn ) = φ1′′ φ 2′′ φ3′′ L φ n′′
M M M M
φ1(n –1) φ 2(n –1) φ3(n –1) φ n (n –1)
Theorem 2.2.2 :
If φ1, φ2 , φ3 ,...., φn are n solutions of L(y) = 0 where L( y ) = y ( n ) + a1 ( x ) y ( n–1) +
a2 ( x) y ( n –2) + ⋅⋅⋅ + an ( x) y, on an interval I, then they are linearly independent on I if and only
if W (φ1, φ 2 , φ3 ,...., φn ) ( x) ≠ 0 for all x in I.
Proof :
Suppose W (φ1, φ 2 , φ3 ,...., φn ) ( x) ≠ 0 for all x in I. We show that φ1, φ2 , φ3 ,...., φn are
n
linearly independent on I. i.e. ∑ ci φi = 0 ⇒ c1 = c2 = c3 = .... = cn = 0
i=1
If there are constants c1, c2, c3,....,cn such that
∑ ci φi ( x) = c1 φ1 ( x) + c2 φ2 ( x) + c3 φ3 ( x) + .... + cn φ n ( x) = 0 for all x in I then clearly,
φ1 ( x ) φ2 ( x ) φ3 ( x ) L φn ( x )
c1
0
φ1′ ( x ) φ 2′ ( x ) φ3′ ( x ) L φ n′ ( x )
c2 0
φ1′′ ( x ) φ 2′′ ( x ) φ3′′ ( x ) L φ n′′ ( x ) c3 = 0
M M M M
M M
φ1(n –1) ( x) φ2 (n –1) ( x) φ3(n –1) ( x) cn 0
φn (n –1) ( x)
φ1 ( x0 ) φ 2 ( x0 ) φ3 ( x0 ) L φn ( x0 )
c1 α1
α
φ1′ ( x0 ) φ 2′ ( x0 ) φ3′ ( x0 ) L φn′ ( x0 )
c2 2
c3 = α3
φ1′′ ( x0 ) φ 2′′ ( x0 ) φ3′′ ( x0 ) L φ n′′ ( x0 )
M M M M M M
cn α n
φ1(n –1) ( x0 ) φ 2(n –1) ( x0 ) φ3(n –1) ( x0 ) φn (n –1) ( x0 )
Since φ1, φ2 , φ3 ,...., φn are linearly independent by theorem 2.2.2, W (φ1, φ 2 , φ3 ,...., φ n )
( x0 ) ≠ 0. Therefore the coefficient matrix is inversible and there is a unique solution
c1, c2, c3,..... cn of the above system of equations.
Thus we have a unique solution
ψ = c1 φ1 ( x ) + c2 φ 2 ( x) + c3 φ3 ( x) + ⋅⋅⋅ + cn φ n ( x )
In theorem 2.2.2 we have seen that the function φ1, φ2 , φ3 ,...., φn are linearly independent
solutions of L(y) = 0 if and only if the Wronskian W (φ1, φ 2 , φ3 ,...., φn ) ( x) ≠ 0 for all x in I. In the
next theorem we show that it is sufficient to calculate the Wronskian W (φ1, φ 2 , φ3 ,...., φ n ) at
some point x0 in I.
φ1 φ2 φ3 L φn
φ1′ φ 2′ φ3′ L φ n′
Let W = W (φ1, φ 2 , φ3 ,...., φn ) = φ1′′ φ 2′′ φ3′′ L φ n′′
M M M M
φ1(n –1) φ 2(n –1) φ3(n –1) φ n (n –1)
φ′ φ 2′ φ3′ L φ n′ φ1 φ2 φ3 L φn
1
φ1′ φ 2′ φ3′ L φ n′ φ1′′ φ 2′′ φ3′′ L φ n′′
W ′ = φ ′′ φ ′′ φ ′′ L φ ′′ + φ ′′ φ ′′ φ ′′ L φ ′′ + ...
1 2 3 n 1 2 3 n
M
M M M M M M M
( n –1) ( n–1) ( n–1) ( n –1) φ ( n –1)
φ ( n –1)
φ ( n –1)
φ ( n –1)
φ1 φ2 φ3 φn 1 2 3 n
φ1 φ2 φ3 L φn
φ1′ φ 2′ φ3′ L φ n′
+ M M M L M
φ (n –2) φ (n –2) φ (n –2) φn ( n –2)
1 2 3
φ1(n ) φ 2( n ) φ3(n ) φ n (n )
= V1 + V2 + V3 + ..... + Vn (say)
Where Vk differs from W only in its kth row and the kth row of Vk is obtained by
differentiating the kth row of W. The first n – 1 determinants are all zero, since they each have
two identical rows. Observe that Vk has kth and (k + 1)th row identical.
Since φ1, φ2 ,...., φ n are solution of L(y) = 0, we have
φ1 φ2 φ3 L φn
φ1′ φ2′ φ3′ L φn′
W ′ = M M M M
φ ( n –2) φ2( n –2) φ3( n –2) L φ1( n–2)
1
– a1φ1( n –1) – a1φ2( n –1) – a1φ3( n –1) L a1 n
– φ ( n –1)
φ1 φ2 φ3 L φn
φ1′ φ 2′ φ3′ L φ n′
= – a1 φ1′′ φ 2′′ φ3′′ L φ n′′ = – a1W(φ1, φ2 , φ3 ,...., φ n )
M M M M
φ1(n –1) φ 2(n –1) φ3(n –1) φn ( n –1)
Therefore W¢ + a1 W = 0 and we get,
x
– ∫ a1 ( t ) dt
W( x) = e x0
W( x0 )
x
i.e. W(φ , φ
1 2 3, φ ,...., φ n ) ( x ) = exp – ∫ 1
a (t ) dt W(φ1, φ2 , φ3 ,...., φn ) ( x0 ).
x0
Corollatory : If the coefficient a1 is constant then
W(φ1 , φ2 , φ3 ,...., φ n ) ( x) = e – a1 ( x – x0 ) W(φ1, φ 2 , φ3 ,...., φ n ) ( x0 ).
From theorem 2.2.2 and theorem 2.2.4 it follows that n solutions φ1, φ2 , φ3 ,...., φn of
L(y) = 0 on I are linearly independent if and only if W(φ1, φ 2 , φ3 ,..., φ n ) ( x0 ) ≠ 0 for some point
x0 in I.
Differential Equations (71)
B. Solutions of non-homogeneous equation
Theorem 2.2.5
Let b (x) be a continuous function on an interval I and let φ1, φ2 , φ3 ,...., φn be a basis for the
solutions of L(y) = 0 on I. Every solution y of L( y ) = y ( n ) + a1 ( x) y ( n –1) + a2 ( x) y ( n –2) +
... + an ( x) y = b( x) can be written as ψ = ψ p + c1φ1 + c2φ2 + c3φ3 + ... + cnφn where is a yp is a
particular solution of L(y) = b (x) and c1, c2, c3,..... cn are constants. Every such y is a solution of
L(y) = b(x). A particular solution yp is given by
n x
Wk (t ) b(t )
ψ p = ∑ φk ( x ) ∫ dt
k =1 x0 W (φ1, φ 2 , φ3 ,..., φ n ) (t )
where W(φ1, φ 2 , φ3 ,..., φ n ) is a wronkian of φ1, φ2 , φ3 ,...., φn and Wk is the determinant obtained
Choose u1, u2, u3,....un such that u1′φ1 + u2′ φ 2 + u3′ φ3 + ... + un′ φn = 0
and ψ p (n ) = u1′φ1(n –1) + u2′ φ2(n –1) + u3′ φ3(n –1) + ... + un′ φn (n –1) + u1φ1(n ) + u2 φ 2(n ) + ... + un φn (n)
If we choose u1′φ1( n –1) + u2′ φ2( n –1) + u3′ φ3( n –1) + ... + un′ φn( n –1) = b( x ) . Then
ψ p (n ) = u1φ1(n ) + u2 φ2 (n) + u3 φ3(n) + ... + un φn (n ) + b ( x)
Thus we have the following equations
ψ p = u1φ1 + u2 φ2 + u3 φ3 + ... + un φn
ψ p′ = u1φ1′ + u2 φ2′ + u3 φ3′ + ... + un φn′ ; u1′φ1 + u2′ φ 2 + u3′ φ3 + ... + u′n φn = 0
ψ p′′ = u1φ1′′ + u2 φ 2′′ + u3 φ3′′ + ... + un φn′′ ; u1′φ1′ + u2′ φ 2′ + u3′ φ3′ + ... + un′ φn′ = 0
ψ p′′′ = u1φ1′′′ + u2 φ 2′′′ + u3 φ3′′′ + ... + un φn′′′ ; u1′φ1′′ + u2′ φ 2′′ + u3′ φ3′′ + ... + un′ φn′′ = 0
M
ψ p (n –1) = u1φ1(n –1) + u 2 φ2 (n –1) + ... + un φn (n –1) ; u1′φ1(n –2) + u2′ φ2(n –2)
φ1 φ2 L φ k –1 0 φ k +1 L φn
φ1′ φ 2′ φ k′ –1 0 φ k′ +1 L φ n′
= b ( x) φ1′′ φ 2′′ φ k′′–1 0 φ k′′+1 L φ n′′ = b( x ) Wk
M M M M M M
φ1( n –1) φ2( n–1) φ k –1( n –1) 1 φk(n+1–1) φ n( n–1)
x
b ( x )Wk b(t ) Wk (t )
Thus, uk′ = i.e. uk = ∫ dt.
W(φ1 , φ 2 , φ3 ,...,φ n ) x0 W (φ1,φ 2 , φ3 ,...,φ n ) (t )
n
and ψ p = u1φ1 + u2 φ 2 + u3 φ3 + ... + un φn = ∑ φk uk
k =1
n x
b(t )Wk (t )
= ∑ φk ( x) ∫ dt.
k =1 x0 W (φ1, φ 2 , φ3 ,..., φ n ) (t )
EXAMPLES
φ1 φ2
= (c1d 2 – c2d1 ) = (c1d 2 – c2 d1 ) W(φ1, φ 2 ) ( x)
φ1′ φ2′
Thus, W(ψ 1,ψ 2 ) ( x) = (c1d 2 – c2 d1 ) W(φ1, φ 2 ) ( x)
Since ψ 1,ψ 2 are independent c1d 2 – c2 d1 ≠ 0.
Therefore there is a non-zero constant k = c1 d2 – c2 d1 such that W(ψ 1,ψ 2 ) ( x) =
k W(φ1 , φ2 ) ( x ).
φ 2′ φ 2
φ 2′ φ 2
[ We use the elementory properties of determinants det A = det AT and if we interchange
row / column, the value of det change its sign.]
2
Q. 4. Find a general solution of y′′ – y = x (0 < x < a)
x2
2
Ans. : Assume that the solution of homogeneous equation L( y ) = y ′′ – y = 0 is of the form
x2
xr. Then y = xr implies L( x r ) = r (r – 1) x r –2 – 2 x r –2 = 0 gives r (r – 1) – 2 = 0. Then
2 1 1
x – 2 – x 2 1 1 1
u1 ( x ) =
′ x 2 1
= 2– 3 and u1 ( x ) = – – – 2
–3 3x 3x 3 x 3 2x
2 1 2
x – 2 (x )
u2′ ( x ) = x 2 1 1 1
= – x+ and u2 ( x) = – x 2 + x
–3 3 3 3 3
Thus ψ p = u1 ( x ) φ1 ( x ) + u2 ( x ) φ2 ( x )
2 1 x2 x 1 1
= – + 2 x2 + – + = –x + .
3x 6 x 3 3 x 2
1 c
The general solution of given non-homogeneous equation is ψ = – x + + c1x 2 + 2 , where
2 x
c1 and c2 are constants.
Ans. : The given non-homogeneous equation is y ′′ – 1 y ′ + 1 y = 1. (We can divide the equation
x x2
by x2 as x2 is positive)
Let φ2 ( x) = u ( x) φ1 ( x) = u ( x) x be an other solution.
1 1 c1
L(φ 2 ) = [u ′′x + 2u ′] – [u ′x + u ] + 2 u ( x) x = 0 gives u′′x + u′ = 0. Therefore u ′ = x
x x
and u ( x) = c1 log x.
φ2 ( x) = u ( x) φ1 ( x) = c1 x log x is second solution. Without loss of generality we choose
c1 = 1.
Thus, φ1 ( x ) = x and φ2 ( x ) = x log x are two solutions of homogeneous equation
1 ′ 1
y′′ –
y + 2 y = 0.
x x
A solution yp of a non-homogeneous equation
1 ′ 1
y′′ – y + 2 y = 1 has the form
x x
ψ p = u1 ( x ) φ1 ( x ) + u2 ( x ) φ2 ( x )
b( x )W1 b( x)W2
Then u1′ ( x ) = , u2′ ( x ) =
W(φ1 , φ2 ) W(φ1 , φ2 )
Q. 7.
(a) Show that there is a basis f1, f2 for the solutions of x 2 y ′′ + 4 xy ′ + (2 + x 2 ) y = 0 ( x > 0)
of the form
ψ1 ( x) ψ 2 ( x)
φ1 ( x) = 2
, φ 2 ( x) =
x x2
(b) Find all solutions of
x 2 y′′ + 4 xy′ + (2 + x 2 ) y = x 2 for x > 0 .
Ans. :
v
(a) Let φ = be a solution of the given homogeneous equation.
x2
4 2
L ( y ) = y ′′ + y ′ + 2 + 1 y = 0
x x
v′ 2v v′′ 4v′ 6v
Then, φ′ = 2
– 3 , φ ′′ = 2 – 3 + 4 and
x x x x x
v′′ 4v′ 6v 4 v′ 2v 2 v
L( y ) = 2 – 3 + 4 + 2 – 3 + 2 + 1 2 = 0.
x x x x x x x x
Therefore L(y) = 0 implies v′′ + v = 0.
ψ1 ( x ) = cos x and ψ 2 ( x ) = sin x are two linearly independent solutions of v′′ + v = 0.
ψ p = u1 ( x) φ1 ( x) + u2 ( x) φ2 ( x).
b ( x )W1 b ( x)W2
Then u1′ ( x ) = , u 2′ ( x ) =
W(φ1 , φ2 ) W(φ1 , φ2 )
where b(x) = 1
cos x sin x
x2 x2 1
W(φ1, φ 2 ) = = 4
sin x 2cos x cos x 2 sin x x
– 2 – –
x x3 x2 x3
sin x cos x
0 2
0
x sin x x2 cos x
W1 = = – 2 , W2 = =
1
cos x 2sin x
– x sin x 2 cos x
– 2 – 1 x2
2 3
x x x x3
– sin x
x2
u1′ ( x ) = = – x 2 sin x, u1 ( x ) = x 2 cos x – 2 x sin x – 2 cos x
1
x4
cos x
u2′ ( x ) = x2 = + x 2 cos x, u2 ( x ) = x 2 sin x + 2 x cos x – 2sin x
1
x4
ψ p = u1 ( x) φ1 ( x) + u2 ( x) φ2 ( x)
cos x sin x
= ( x 2 cos x – 2 x sin x – 2 cos x ) 2
+ ( x 2 sin x + 2 x cos x – 2sin x )
x x2
2
= 1–
x2
Therefore the general solution of non-homogeneous equation is
2 cos x sin x
ψ = ψ p + c1 φ1 + c2 φ2 = 1– 2
+ c1 2
+ c2 .
x x x2
b( x )W1 b( x )W2
u1′ ( x ) = , u2′ ( x ) =
W1 (φ1, φ2 ) W1 (φ1 ,φ 2 )
cos x sin x
W (φ1, φ2 ) = =1
– sin x cos x
0 sin x cos x 0
W1 = = – sin x, W2 = = cos x
1 cos x – sin x 1
x x
b(t )W1 (t ) – b(t ) sin t
u1 ( x ) = ∫ dt = ∫ dt
1 W( φ 1 , φ 2 ) (t ) 1 1
x x
b(t )W2 (t ) b(t ) cos t
u2 ( x ) = ∫ dt = ∫ dt
1 W(φ1, φ2 ) (t ) 1 1
x x
ψ p = – cos x ∫ b (t ) sin t dt + sin x ∫ b(t ) cos t dt
1 1
x
= ∫ b(t ) [sin x cos t – cos x sin t ] dt
1
x
= ∫ b(t ) sin ( x – t ) dt
1
EXERCISE
4. Use the method of variation of parameter and find the particular solution of the following
equations where the solutions for the related homogeneous equation are given.
Ans. :
ψ p = x3 log x – x3
2 2 1 3
(a) y ′′ – y ′ + 2 y = x log x, φ1 ( x ) = x, φ2 ( x) = x 2
x x 2 4
1 Ans. : ψ p = x3
(b) x y′′ + xy′ – 4 y = x , φ1 = x , φ2 =
2 3 2
x2 5
∞
If x0, x and an are real numbers the series ∑ an ( x – x0 ) is called a real power series. Its
n
n =0
disk of convergence intersects the real axis in an interval (x0 – r, x0 + r) called the interval of
convergence.
If g is a function defined on an interval I containing point x0 we say that g is analytic at x0
if g can be expanded in a power series about x0 which has a positive radius of convergence. Thus
g is analytic at x0 if it can be representd in the form
Differential Equations (82)
∞
g ( x ) = ∑ an ( x – x0 ) n
n =0
Where an are constants and the series converges for | x – x0 | < r , r > 0. If g has a power
series expansion then all the derivatives of g exist on | x – x0 |< r and they may be computed by
differentating the series term by term that is
∞ ∞
g ′( x ) = ∑ n an ( x – x0 ) n –1 , g ′′( x) = ∑ n (n – 1) an ( x – x0 ) n –2 etc.
n =0 n =0
The differentiated series converges on |x – x0 | < r.
In calculus there are certain tests by which one could determine an interval of converge of
a real power series. A simple one and one which is frequently used is known as ratio test.
∞
un+1
The series ∑ ui converges absolutely if
Lt
n→∞ = k < 1.
n =0 un
∞
xn
Example 1 : For the power series ∑ ,
n =0 n
un+1 x n +1
Lt Lt n+1 n
n→∞ = n→∞ = nLim
→∞ ⋅ x =| x |
un xn
n
n +1
Hence the series converges absoulately if | x | < |.
x 2 x 4 x6 (–1) n –1 2n –2
1– + – + .... + x + ....,
2! 4! 6! (2n – 2)!
(–1) n –1 2 n –2 (–1) n 2 n
un = x and un+1 = x
(2n – 2)! (2 n)!
Therefore
un+1 (–1) n x 2 n (2n – 2)! x2
Lim
n→∞ = Lim
n→∞ × = Lim
n→∞ = 0 for each x
un (2n )! (–1) n –1 x 2 n –2 2n(2n – 1)
Hence the series converges absoulately for all x. Its interval of convergence is the entire
real axis.
have convergent power series expansions in powers of (x – x0) on an interval | x – x0 |< r , r > 0.
If a1, a2, a3,.....an are any n constants, there exists a solution f of the problem
L( y ) = 0, y ( x0 ) = α1 , y′( x0 ) = α 2 ,...., y (n –1) ( x0 ) = α n
Differential Equations (83)
with a power series expansion
∞
φ ( x) = ∑ ck ( x – x0 ) k
k =0
is convergent for | x – x0 |< r . We have
Above equation is true only if all the coefficients of the power series of x are zero. Thus,
2a2 = 0, (n + 2) (n + 1) an+ 2 – an –1 = 0, n = 1, 2, 3....
This gives an infinite set of equations, and can be solved for an. Thus, for n = 1 we have
a0
(3) ⋅ (2) ⋅ a3 = a0 or a3 =
(3) ⋅ (2)
For n = 2 we find
a1
(4) ⋅ (3) a4 = a1 or a4 =
(4) ⋅ (3)
Differential Equations (84)
Continuing in this way we see that
a0; a1 a2 =0
a a1 a2
a3 = 0 ; a4 = ; a5 = =0
(3)⋅(2) (4)⋅(3) (5)⋅(4)
a a0 a a1
a6 = 3 = ; a7 = 4 = ; a8 =0
(6)⋅(5) (6)⋅(5)⋅(3)⋅(2) (7)⋅(6) (7)⋅(6)⋅(4)⋅(3)
a a0 a a1
a9 = 6 = ; a10 = 7 = ; a11 =0.
(9)⋅(8) (9)⋅(8)⋅(6)⋅(5)⋅(3)⋅ (2) (10)⋅(9) (10)⋅(9)⋅(7)⋅ (6)⋅ (4)⋅ (3)
In general
a0
a3m = ;
(2) ⋅ (3) ⋅ (5) ⋅ (6) ⋅ (8) ⋅ (9) ⋅⋅⋅ (3m – 1) (3m )
a1
a3m +1 = ; a3m + 2 = 0
(3) ⋅ (4) ⋅ (6) ⋅ (7) ⋅ (9) ⋅ (10) ⋅⋅⋅ (3m) (3m + 1)
Thus all the constants are determined in terms of a0 and a1. Collecting together terms
containing a0 and a1 as a factor we have
x3 x6 x4 x7
φ ( x ) = a0 1 + + + ... + a1 x + + + ...
(3) ⋅ (2) (2) ⋅ (3) ⋅ (5) ⋅ (6) (4) ⋅ (3) (3) ⋅ (4) ⋅ (6) ⋅ (7)
Let f1 and f2 represent the two series in the brackets.
∞ x 3m
Thus, φ1 ( x ) = 1 + ∑ ,
m =1 2 ⋅ 3 ⋅ 5 ⋅ 6 ⋅ 8 ⋅ 9....(3m –1) (3m)
∞ x 3m+1
φ2 ( x ) = x + ∑ .
m =1 3 ⋅ 4 ⋅ 6 ⋅ 7 ⋅ 9 ⋅ 10....(3m) (3m + 1)
We have shown, in a formal way that f satisfies y ′′ – xy = 0 for any two contants a0 and a1
In particular the choice a0 = 0 and a1 = 1 implies f2 (x) satisfies the equation and a0 = 1, a1
= 0 implies f1 (x) satisfies the equation.
The only question that remains is about the convergence of the series, defining f1 (x) and
f2 (x).
∞ x 3m
φ1 ( x ) = 1 + ∑ d m ( x ) = 1 + ∑
m =1 2 ⋅ 3 ⋅ 5 ⋅ 6 ⋅ 8 ⋅ 9....(3m – 1) (3m)
1. Find two linearly independent power series solutions (in powers of x) of the following
equations.
(a) y ′′ − xy ′ + y = 0 (b) y ′′ + 3x 2 y ′ – xy = 0
(c) y′′ – x 2 y = 0 (d) y′′ + 3x 3 y′ + x 2 y = 0
∞
Ans. (a) : Let φ ( x ) = a0 + a1x + a2 x 2 + a3 x 3 + ... + an x n + ... = ∑ an x n be a solution of
n =0
L ( y ) = y′′ – xy′ + y = 0 . Since it is a solution it satisfies the equation L(f ) = 0.
∞
φ ( x) = a0 + a1x + a2 x 2 + a3 x3 + ... + an x n + ... = ∑ an x n
n=0
∞
Then φ ′ ( x) = a1 + 2a2 x + 3a3 x 2 + ... + n an x n –1 + ... = ∑ n an x n –1
n =1
∞
φ ′′ ( x ) = 2a2 + 3 ⋅ 2a3 x + ... + n (n – 1) an x n –2 + ... = ∑ n (n – 1) an x n –2
n=2
∞ ∞ ∞
Thus, L(φ ) = ∑ n (n – 1) an x n –2 – x ∑ n an x n –1 + ∑ an x n = 0.
n=2 n=1 n =0
∞ ∞ ∞
= ∑ (n + 2) (n + 1) an+ 2 x n – ∑ n an x n + ∑ an x n = 0
n =0 n=1 n =0
∞
= (2 a2 + a0 ) + ∑ {(n + 2) (n + 1) an+2 – nan + an }x n = 0
n=1
∞
L(φ ) = (2a2 + a0 ) + ∑ {(n + 2) (n + 1) an+ 2 – (n – 1)an }x n = 0
n=1
We see that L(f ) = 0 if and only if 2a2 + a0 = 0 and (n + 2) (n + 1) an+ 2 – (n –1) an = 0 for
1 (n – 1) an
n = 1, 2,3,.... a2 = – a0 ; an + 2 = is called recurrence relation.
2 (n + 2) (n + 1)
a0 ; a1
– a0
a2 = ; ; a3 = 0.a1
2 ⋅1
a a 2
a4 = 2 = – 0 ; a5 = .0 = 0
4⋅3 2⋅3⋅ 4 5⋅ 4
3a 3a0
a6 = 4 = – ; a7 = 0
6⋅5 2 ⋅ 3 ⋅ 4 ⋅ 5⋅ 6
3 ⋅ 5 a0
a8 = – ; a9 = 0
2 ⋅ 3⋅ 4 ⋅ 5⋅ 6⋅ 7 ⋅ 8
3 ⋅ 5 ⋅ 7 a0
a10 = – ; a11 = 0
10!
Differential Equations (86)
In general a2n+1 = 0 n = 1, 2,3,....
3 ⋅ 5 ⋅ 7 ⋅ 9 ⋅⋅⋅ (2n – 3) a0
a2 n = –
(2n)!
2 ⋅ 3 ⋅ 4 ⋅ 5 ⋅ 6 ⋅ 7 ⋅ 8 ⋅ 9 ⋅⋅⋅ (2n – 3) (2n – 2) (2n – 1) (2n)
=– a0
2 ⋅ 4 ⋅ 6 ⋅ 8 ⋅⋅⋅⋅ 2n (2n –1) (2n)!
(2n) !
=– n
a0
2 n !(2n – 1) (2n)!
a0
=– n
2 n !(2n –1)
∞ ∞
φ ( x ) = ∑ a2n x 2n + ∑ a2n+1 x 2n+1
n =0 n =0
∞ a0
= a0 – ∑ n
x 2 n + a1x
n =1 2 n !(2n – 1)
∞
x 2n
= a0 1 – ∑ n + a1 x
n =1 2 n !(2n –1)
∞ x 2n
φ1 ( x) = 1– ∑ and φ2 ( x) = x are two solutions of the equation
n=1 2 n n !(2n –1)
∞
Let φ1 ( x ) = ∑ d m ( x )
m =0
x 2(n+1)
d n +1 2n+1 (n + 1)!(2n + 1) x 2 (2n – 1)
= =
dn x 2n 2(n + 1) (2n + 1)
2 n n !(2n –1)
Lt d n +1 Lt (2 n – 1)
n→∞ = n→∞ =0
dn 2(n + 1) (2n + 1)
Redius of convergence = ¥
The series converges if | x | < ¥ i.e. all values of x. Both the solutions are convergent for
all values of x.
∞
Ans. (b) : Let φ ( x ) = ∑ an x n be a solution.
n =0
∞ ∞
φ ′ ( x) = ∑ nan x n –1 , φ ′′( x) = ∑ n (n – 1)an x n –2
n=1 n=2
L(φ ) = φ ′′ + 3 x 2φ ′ – xφ
∞
Ans. (c) : Let φ ( x) = ∑ an x n be a solution of y′′ – x 2 y = 0. Since it is a solution f (x) satisfies
n= 0
L(φ ) = φ ′′ – x 2φ = 0.
∞ ∞ ∞
φ ( x) = ∑ an x n , φ ′ ( x) = ∑ n an x n –1 , φ ′′ ( x) = ∑ n (n –1) an x n –2
n= 0 n =1 n= 2
∞ ∞
L(φ ) = ∑ n(n –1) an x n –2 – x 2 ∑ an x n = 0.
n= 2 n= 0
∞ ∞
2 ⋅1a2 + 3 ⋅ 2 a3 x + ∑ n (n –1)an x n –2 – ∑ an x n+2 = 0
n= 4 n= 0
an
Thus, a2 = 0, a3 = 0 and an + 4 = .
(n + 4) (n + 3)
a0 a1 a2 = 0 a3 = 0
a0 a1
a4 = ; a5 = ; ; a6 = 0 ; a7 = 0
3⋅ 4 5⋅4
a a0 a a1
a8 = 4 = ; a9 = 5 = ; a10 = 0 ; a11 = 0
8 ⋅ 7 3⋅ 4 ⋅7 ⋅8 9 ⋅8 4 ⋅5 ⋅8 ⋅9
a0 a1
a12 = ; a13 = ; a14 = 0 ; a15 = 0
3 ⋅ 4 ⋅ 7 ⋅ 8 ⋅ 11⋅ 12 4 ⋅ 5 ⋅ 8 ⋅ 9 ⋅ 12 ⋅ 13
M M M M
Thus all the coefficients an’s are determined in terms of a0 and a1 since a2 = a3 = 0 implies
a4m+2 and a4m+3 = 0 for m = 0, 1, 2, 3,..... Therefore
∞ ∞
φ ( x) = ∑ a4m x m + ∑ a4m +1 x 4 m+1
m=0 m=0
∞
x 4m
= a0 1 + ∑
m = 1 3 ⋅ 4 ⋅ 7 ⋅11 ⋅12 ⋅⋅⋅ (4m – 1) (4m)
∞
x 4m +1
+ a1 x + ∑
m = 1 4 ⋅ 5 ⋅ 8 ⋅ 9 ⋅12 ⋅13 ⋅⋅⋅ (4m) (4m + 1)
Therefore two linearly independent solutions are
∞
x 4m
φ1 ( x ) = 1 + ∑ and
m = 1 3 ⋅ 4 ⋅ 7 ⋅11⋅12 ⋅13 ⋅⋅⋅ (4m –1) (4m)
∞ x 4m+1
φ2 ( x ) = x + ∑ .
m = 1 4 ⋅ 5 ⋅ 8 ⋅ 9 ⋅12 ⋅13 ⋅⋅⋅ (4m) (4m + 1)
∞
Ans. (d) : Let φ2 ( x) = ∑ an xn be a solution of L( y ) = y′′ + x3 y′ + x 2 y = 0. Therefore
n =0
∞ ∞ ∞
L(φ ) = ∑ n (n – 1) an x n –2 + x 3 ∑ n an x n –1 + x 2 ∑ an xn = 0
n=2 n =1 n=0
∞ ∞ ∞
that is ∑ n (n – 1) an xn –2 + ∑ n an xn + 2 + ∑ an x n+ 2 = 0.
n=2 n =1 n=0
a0 a1 a2 = 0 a3 = 0
– a0 –2 a1
a4 = ; ; a5 = ; a6 = 0 ; a7 = 0
3⋅ 4 5⋅ 4
5 a0 6 ⋅ 2 a1
a8 = + ; a9 = ; a10 = 0 ; a11 = 0
3⋅ 4 ⋅ 7 ⋅8 4 ⋅5 ⋅8 ⋅9
9 ⋅ 5 a0 10 ⋅ 6 ⋅ 2 a1
a12 = – ; a13 = – ; a14 = 0 ; a15 = 0
3 ⋅ 4 ⋅ 7 ⋅ 8 ⋅ 11⋅ 12 4 ⋅ 5 ⋅ 8 ⋅ 9 ⋅ 12 ⋅ 13
M M M M
(–1) m 5 ⋅ 9 ⋅ 13 ⋅⋅⋅ (4m – 3)
a4m = ;
3 ⋅ 4 ⋅ 7 ⋅ 8 ⋅ 11⋅ 12 ⋅⋅⋅ (4m – 1) (4m)
(–1) m 2 ⋅ 6 ⋅ 10 ⋅⋅⋅ (4m – 2)
a4m +1 = ; a4 m + 2 = a4m +3 = 0.
4 ⋅ 5 ⋅ 8 ⋅ 9 ⋅⋅⋅⋅(4 m) (4m + 1)
Therefore two linearly independent solutions are
∞ (–1) m 5 ⋅ 9 ⋅13 ⋅⋅⋅ (4m – 3)
φ1 ( x ) = 1 + ∑ x 4m
m =1 3 ⋅ 4 ⋅ 7 ⋅ 8 ⋅ 11⋅ 12 ⋅⋅⋅ (4m –1) (4m)
∞ (–1) m 2 ⋅ 6 ⋅ 10 ⋅⋅⋅ (4m – 2) 4 m+1
φ2 ( x ) = x + ∑ x
m =1 4 ⋅ 5 ⋅ 8 ⋅ 9 ⋅⋅⋅ (4m) (4m + 1)
∞
2. Find the solution f of y ′′ + ( x – 1) 2 y ′ – ( x – 1) y = 0 in the form φ ( x) = ∑ ak ( x – 1) k which
k =0
satsfies φ (1) = 1, φ ′(1) = 0.
∞
Ans. : Let φ ( x) = ∑ ak ( x –1) be a solution of L( y ) = y′′ + ( x –1) 2 y′ – ( x –1) y = 0.
k
k =0
φ ′(1) = 0 give a1 = 0.
∞
Ans. : Let φ ( x) = ∑ an x n be a solution of L( y ) = y′′′ – xy = 0. Then
k =0
∞ ∞
L(φ ) = ∑ n (n –1) (n – 2) an x n –3 – x ∑ an xn = 0
n= 3 n= 0
∞ ∞
Then 3 ⋅ 2 ⋅1 a3 + ∑ n ( n –1) (n – 2) an x n –3 – ∑ an xn +1 = 0
n= 4 n= 0
In the first sum replace n by n + 4, then
∞
3 ⋅ 2 ⋅1 a3 + ∑ [ ( n + 4) ( n + 3) ( n + 2) an+4 – an ] x n +1 = 0
n= 0
an
Thus, a3 = 0 and an + 4 = .
(n + 4) (n + 3) (n + 2)
a0 a1 ; a2 ; a3 = 0
a0 a1 a2
a4 = ; a5 = ; a6 = ; a7 = 0
4 ⋅ 3⋅ 2 5⋅ 4⋅ 3 6⋅ 5⋅ 4
a0 a1 a2
a8 = ; a9 = ; a10 = ; a11 = 0
8⋅ 7 ⋅ 6 ⋅ 4 ⋅3⋅ 2 9 ⋅8 ⋅7 ⋅5 ⋅ 4 ⋅3 10 ⋅ 9 ⋅ 8 ⋅ 6 ⋅ 5 ⋅ 4
a0 a1 a2
a12 = ; a13 = ;a =
2 ⋅ 3 ⋅ 4 ⋅ 6 ⋅ 7 ⋅ 8 ⋅10 ⋅11 ⋅12 3 ⋅ 4 ⋅ 5 ⋅ 7 ⋅ 8 ⋅ 9 ⋅11 ⋅12 ⋅13 14 4 ⋅ 5 ⋅ 6 ⋅ 8 ⋅ 9 ⋅ 10 ⋅ 12 ⋅13 ⋅ 14
;a15 = 0
M M M M
a0 a1
a4 m = ; a4m +1 =
2 ⋅ 3 ⋅ 4 ⋅ 6 ⋅ 7 ⋅ 8 ⋅⋅⋅ (4m – 2) (4m –1) (4m) 3 ⋅ 4 ⋅ 5 ⋅⋅⋅ (4m – 1) (4m) (4m + 1)
a2
; a4 m + 2 = ; a4 m +3 = 0
4 ⋅ 5 ⋅ 6 ⋅ 8 ⋅ 9 ⋅ 10 ⋅⋅⋅ (4m) (4m + 1) (4m + 2)
The general solution f (x) of the given equation contains three parameters a0, a1, a2.
The solution f (x) becomes
∞ x 4m
φ ( x) = a0 1 + ∑
m = 1 2 ⋅ 3 ⋅ 4 ⋅ 6 ⋅ 7 ⋅ 8 ⋅⋅⋅ (4 m – 2) (4 m – 1) (4 m)
∞
x 4 m+1
+ a1 x + ∑
m = 1 3 ⋅ 4 ⋅ 5 ⋅ 7 ⋅ 8 ⋅ 9 ⋅⋅⋅ (4m – 1) (4m) (4m + 1)
For f to satisfy L(f ) = 0 we must have all the coefficients of the powers of x equal to zero.
Since Wronkian W(φ1, φ2 ) ≠ (0), φ1, φ2 are linearly independent and therefore forms a basis.
If a is a non-negativ even integer a = 2n, then f1 has only a finite number of non-zero
terms. In this case f1 is a polynomial of degree 2n containing only even powers of x. for example,
α = 0, φ1 ( x ) = 1 = p0 ( x )
(–1) (2 + 1)
α = 2, φ1 ( x) = 1 + α x 2 = 1– 3 x 2 = p2 ( x)
2
or the recurrance relation
n( n + 1) – α (α + 1)
an + 2 = an implies
(n + 1) (n + 2)
0(0) – 2 (3)
a2 = a0 = –3a0
1⋅ 2
2(3) – 2 (3)
a4 = a2 = 0
3⋅ 4
with a0 = 1 we get φ1 ( x ) = 1 – 3 x 2 = p2 ( x )
for a = 4,
0(0) – 4(5)
a2 = a0 = –10a0
1⋅ 2
2(3) – 4(5)
a4 = a2
3⋅ 4
6 – 20
= (–10 a0 )
12
140
=+ a0
12
35
= a0
3
4(5) – 4 (5)
a6 =
5⋅ 6
=0
3 5 4
φ1 ( x) = a0 1 – 10 x 2 + x with a0 = 1
3
35 4
φ1 ( x) = 1– 10 x 2 + x = p4 ( x)
3
The solution f2 is not a polynomial in this case since none of the coefficients in the series
of f2 vanish.
A similar situation occurs when a is a positive odd integer n. Then f2 is a polynomial of
degree n having only odd powers of x and f1 is not a polynomial.
Satisfying Pn (1) = 1 is called the nth Legendre polynomial and the differential equation is
called Legendre equation.
Let f be a polynomial of degree n defined by
dn
φ ( x) = n
( x 2 – 1) n
dx
u ( x ) = ( x 2 –1)n implies u ′( x) = n ( x 2 –1) n –1 2 x gives
( x 2 –1) u′( x) – 2n x u ( x) = 0
(1 – x 2 )u (n+ 2) – 2 x u (n+1) + n (n + 1) u (n ) = 0.
dn 2 n dn
Since φ ( x) = n
( x – 1) = n
u ( x) = u ( n) ( x),
dx dx
(1– x 2 ) φ ′′ – 2 x φ ′( x) + n (n + 1) φ ( x) = 0
dn 2
Thus the function φ ( x) = n
( x – 1) n is a solution of Legendre equation.
dx
dn 2 dn
( x –1) n
= n
( x + 1) n ( x – 1) n
dx n
dx
d n d n –1 d
= n ( x – 1)n ( x + 1)n + n –1 ( x – 1)n ( x + 1)n + ...
dx
dx dx
= n (n –1) (n – 2)....2 ⋅1( x + 1) n + terms containing (x – 1) as factor.
= n ! (x + 1)n + terms containing (x – 1) as factor.
Thus, at x = 1,
dn
n
( x 2 – 1) n = 2 n ⋅ n !
dx
1 1 dn 2
Define Pn ( x) = n
φ ( x ) = n n
( x –1) n then Pn (x) is a solution of Legendre
2 n! 2 n ! dx
1
equation witha = n Pn (1) = n 2 n ⋅ n ! = 1. Thus, Pn(x) is a Legendre polynomial of degree n.
2 n!
Supposey is a polynomial solution of Legendre equation with a = n. Since f1 and f2 are
basic solutions of Legendre equation ψ = c1 φ1 + c2 φ2 on | x | < | for some constants c1 and c2 is
a solution. If n is even f1 is polynomial solution and f2 is not a polynomial ψ – c1 φ1 is polynomial
EXERCISES
(v) y ′′ + x 2 y = 0
Answers :
∞ (–3) k x 2k
3. (i) y x a0 ∑
( ) = 1 + k
k = 1 2 k !
∞ (–3)k x 2k +1
+ a1 x + ∑
m =1 3 ⋅ 5 ⋅ 7, ⋅⋅⋅(2k + 1)
∞ (–1) k +1 2 2k x 2 k +1
(ii) y = a0 (1 + 4 x 2 ) + a1 ∑
k =0 4k 2 – 1
x3
(iii) y ( x ) = a0 (1– 3 x 2 ) + a1 x –
3
1 ∞ 3(–1)k x 2 k +1
(iv) y ( x) = a0 1 + x 2 + x 4 + a1 ∑ 2 k
12 k = 0 2 k !(2k – 3) (2k – 1) (2k + 1)
∞ (–1) k x 4k
(v) y x a0 ∑ 2 k
( ) = 1 +
k = 1 2 k ! 3 ⋅ 7 ⋅11⋅⋅⋅ (4k –1)
∞ (–1) k x 4 k +1
+ a1 x + ∑ 2k
k = 1 2 k ! 5 ⋅ 9 ⋅13 ⋅⋅⋅ (4k + 1)
Contents :
Unit 1 : Euler equation
Unit 2 : Second order equations with regular singular points
Unit 3 : The Bessel equation
Unit 4 : Regular singular points at infinity
Introduction
The simplest example of a second order equation that follows defination 3.1.1(a) is the
Euler equation
L( y) = x 2 y′′ + a x y′ + b y = 0
where a, b are constants.
Theorem 3.1.1
Consider the second order Euler equation
L( y) = x 2 y′′ + a x y′ + b y = 0 (a, b constants),
and the polynomial q given by
q ( r ) = r ( r –1) + ar + b
A basis for the solutions of the Euler equation on any interval not containing x = 0 is given
by
φ1 ( x ) = | x |r1 , φ2 ( x ) = | x |r2 ,
in case r1, r2 are distinct roots of q and by
φ1 ( x) = | x |r1 , φ 2 ( x) = | x |r1 log | x |,
if r1 is a root of equation q of multiplicity two.
Proof :
Case 1 : r1 ¹ r2
(a) We first consider the equation for x > 0. Let xr be a solution of Euler equation
L( y) = x 2 y′′ + a x y′ + b y = 0
L( x r ) = x 2 r (r – 1) x r –2 + a x r x r –1 + bx r = 0
(– x )r ′ = – r (– x )r –1 , (– x) r ′′ = r (r – 1) (– x )r –2
′
x (– x )r = r (– x ) (– x )r –1 = r (– x )r and
L( y) = r (r –1)(– x) r + a r (– x)r + b (– x) r = q (r )(– x) r
if r1 ¹ r2 then
x r = e r log x ( x > 0)
r r
then ( x r )′ = (r ⋅ log x )′ ⋅ e r log x = ⋅ x = r x r –1 and the result follows on the same lines
x
for complex roots also.
Thus, we have proved that if (x > 0) r1 ¹ r2 φ1 ( x) = x r1 and φ2 ( x) = x r2 are solutions of
L(y) = 0 and for x < 0, r1 ¹ r2 we have φ1 ( x ) = (– x) r1 and φ2 ( x ) = (– x )r2 are solution of L(y) = 0.
Since | x | = x for x > 0 and | x | = – x for x < 0 φ1 ( x) =| x |r1 and φ2 ( x ) =| x |r2 are solutions of
L(y) = 0 if r1, r2 are distinct roots of q(r) = 0.
We prove that f1 and f2 are linearly independent.
Let c1φ1 + c2 φ 2 = 0 i.e. c1 | x |r1 + c2 | x |r2 = 0 then c1 + c2 | x |r2 – r1 = 0 for every x ÎR.
Differentiating above equation w.r.t. x for x > 0 or x < 0 we get,
c2 (r2 – r1 )| x |r2 – r1 = 0
But r1 ¹ r2 and x ¹ 0 therefore c2 = 0 and c1φ1 + c2 φ 2 = 0 for all x implies c1 = 0 since
c2 = 0 and φ1 ( x) ≠ 0 .
Thus, f1 and f2 are linearly independent solutions.
Therefore if r1 and r2 are distinct roots of q(x) = 0, then φ1 ( x ) =| x |r1 and φ2 ( x ) =| x |r2
forms a basis for the solutions of L(y) = 0.
Differential Equations (102)
Case 2 : r1 = r2
(a) x > 0 : If r1 = r2 then q(r1) = 0 and q¢ (r1) = 0. We have proved that if r1 is a root of q(x) = 0
then φ1 ( x) = x r is a solution. To construct second solution consider
∂ ∂
L( x r ) = q (r ) x r
∂r ∂r
= [q′(r ) + q (r ) log x ] x r
∂ r
Since, x = x r log x
∂r
But if r1 = r2 = r then q (r) = 0 and q¢ (r) = 0 and we have
∂
L( x r ) = 0 .
∂r
∂ ∂
L( x r ) = L x r = L( x r log x)
∂r ∂r
Thus, L( x r log x) = 0 implies xr log x is a solution of L(y) = 0.
If r1 is a root of q (r) = 0 of multiplicity two then φ1 ( x) = x r1 and φ2 ( x) = x r1 log x are two
solutions of L(y) = 0.
(b) x < 0 : If x < 0, then – x > 0 and φ1 ( x ) = (– x ) r1 and φ2 ( x ) = (– x )r1 log (– x ) are solution of
L(y) = 0.
Thus φ1 ( x ) = | x |r1 and φ2 ( x ) = | x |r1 log | – x | are two solution of L(y) = 0.
c1φ1 + c2φ 2 = 0 implies c1 + c2 log | x | = 0 for all x and therefore c1 = c2 = 0 and f1, f2 are
linearly independent.
Thus if r1 is a repeated root of q(r) = 0 then φ1 ( x ) =| x |r1 and φ2 ( x) =| x |r1 log | x | is a
basis for solutions of the Euler equation L( y) = x 2 y′′ + a xy′ + by = 0.
Illustration :
x 2 y ′′ + xy′ + y = 0 for x ¹ 0 is Euler equation with a = b = 1.
The polynomial q(r ) = r ( r – 1) + r + 1 = r 2 + 1 and r = +i, – i are roots of q(r). A basis for
the solutions by theorem 3.1.1 are
φ1 ( x) =| x |i and φ2 ( x) =| x |– i ( x ≠ 0)
Theorem 3.1.2
Consider the Euler equation of order n.
L( y ) = x n y (n ) + a1x n –1 y (n –1) + a2 x n –2 y (n –2) + ... + an y = 0,
Differential Equations (103)
where a1, a2, a3,..., an are constants. Let r1, r2,..., rs be distinct roots of the indicial polynomial
q( r ) = r ( r – 1) (r – 2)...( r – n + 1) + a1 r ( r – 1)...( r – n + 2) + ... + an and suppose
ri has multiplicity mi. Then the n functions
| x |r1 , | x |r1 log | x |,....,| x |r1 (log | x |) m1 –1; | x |r2 ,| x |r2 log | x |,...,| x |r2 (log | x |) m2 –1;...;
| x |rs , | x |rs log | x |,...., | x |rs (log | x |) ms –1
form a basis for the solution of L(y) = 0 on any interval not containing zero.
Proof : Let | x |r be a solution of L(y) = 0.
(| x |r )′ = r | x |r –1 , (| x |r )′′ = r (r –1) | x |r –2 ,....
(| x |r )( n) = r (r – 1) ( r – 2)...(r – n + 1)| x |r – n
Hence, L(| x |r ) = r (r – 1) (r – 2)...(r – n + 1)| x |r + a1r (r – 1) (r – 2)
...(r – n + 2) | x |r +... + an | x |r
= q ( r ) | x |r
where q (r ) = (r ) (r – 1) (r – 2)...(r – n + 1) + r (r – 1) (r – 2)...(r – n + 2)a1 + ... + an .
The polynomial q(r) is called indicial polynomial. Thus, | x |r is a solution of L(y) = 0 if
q (r) = 0 i.e. if r is a root of indicial polynomial then | x |r is a solution of L(y) = 0.
Differentiating L(| x |r ) = q(r ) | x |r with respect to ‘r’ we get
∂ ∂
L (| x |r ) = L | x |r
∂r ∂r
= ( q ′( r ) + q ( r ) log | x | ) | x |r
In general k times differentiation gives
∂k r ∂k r
k
L (| x | ) = L k | x|
∂r ∂r
= q( k ) (r ) + kq ( r –1) (r ) log | x | +k (k –1)q( k –2) ( x) (log | x |)2 + ... + q(r ) (log | x |)k | x |r .
If r is a root of q (r) with multiplicity (k + 1) then q ( r ) = 0, q′( r ) = 0, q ′′( r ) = 0...,
∂i r
q (k )
(r ) = 0 and therefore L i | x | = 0 for i = 1, 2,3,...., k .
∂r
i
∂ r
Thus φ ( x) = i | x | , i = 1, 2,3,...., k. are solution of L(y) = 0.
∂r
If r1 is a root of q(r) of multiplicity m1 then
∂ ∂2 ∂ m1 –1
| x |r1 , | x |r1 =| x |r1 log | x |, 2 | x |r1 =| x |r1 (log | x |) 2 ,..., m –1 | x |r1
∂r1 ∂r1 ∂r1 1
x3 x 3 log x
W (φ1, φ2 ) = = x5 ,
2 2 2
3x x + 3x log x
0 x 3 log x
W1 = = – x 3 log x,
2 2
1 x + 3x log x
– x3 log x ⋅ x log x 1
u1 ( x ) = ∫ 5
dx = – ∫ dx = – (log x ) 2
x x 2
x3 ⋅ x dx dx
u2 ( x ) = ∫ 5
= ∫ = log x
x x
ψ ( x ) = u1 ( x ) φ1 ( x ) + u2 ( x ) φ2 ( x )
1 1
= – (log x) 2 x 3 + (log x) x 3 log x = x 3 (log x) 2
2 2
The general solution
1
φ = c1φ1 + c2φ2 + ψ = c1 x 3 + c2 x 3 log x + x 3 (log x) 2 .
2
(e) The indicial equation
q( r ) = r (r –1) (r – 2) + 2r ( r – 1) – r + 1
Since both the roots are distinct, φ1 ( x) =| x |2i and φ2 ( x) =| x | –2i . The general
solution of homogeneous equation is
φ ( x) = c1 | x |2i + c2 | x |–2i
The particular solution will be calculated by variation of constant method.
Case 1 : x > 0,
If x > 0 then | x | = x
φ1 ( x ) = x 2i and φ2 ( x ) = x –2i
Let ψ ( x ) = u1 ( x ) φ1 ( x ) + u 2 ( x) φ2 ( x ) be a solution of
x 2i x –2i 4i
W(φ1, φ2 ) = = –2i x –1 – 2i x –1 = –
2i –1 –2i –1 x
2i x –2ix
0 x –2i –2i x 2i 0
W1 ( x) = =–x ; W2 = = x 2i
–2i –1 2i –1
1 –2ix 2i x 1
1
– x –2i 2 –2i –2+1
x dx = x x –2 i x –2i
u1 ( x ) = ∫ ∫ 4i dx = =
4i 4i (–2i) 8
–
x
1
x 2i 2
x x 2i –2+1 x2i x 2i
=
u2 ( x ) ∫ =
dx ∫ dx = = .
4i – 4i – 4i (+2i ) 8
–
x
x –2i 2i x 2i –2i 1
Thus, ψ ( x ) = u1 ( x ) φ1 ( x ) + u2 ( x )φ 2 ( x ) = ⋅x + ⋅ x = . For x > 0,
8 8 4
–2i 1
φ ( x ) = c1x 2i + c2 x + is a solution of given equation.
4
Case 2 : x < 0
If x < 0 then | x |= – x and φ1 ( x) = (– x) 2i and φ2 ( x) = (– x) –2i
Let ψ ( x ) = u1 ( x ) φ1 ( x ) + u 2 ( x) φ2 ( x ) be a solution of the given differential equation.
(– x )2i (– x ) –2i 2i 2i 4i
W(φ1, φ2 ) = = + =
2i –1 –2i –1 (– x ) (– x ) – x
–2i(– x) +2i(– x)
0 (– x) –2i (– x) 2i 0
W1 = = –(– x) –2 i , W2 = = (– x) +2 i ,
–2i –1 2i –1
1 2i (– x ) –2i (– x ) 1
1
b ( x) =
x2
1
–(– x ) –2 i ⋅
W1 ( x) b ( x) (+ x ) 2 1 (– x ) –2i (– x ) –2i
u1 ( x ) = ∫ dx = ∫ dx = – = .
W(φ1, φ2 ) 4i 4i 2i 8
(– x )
distinct φ1 ( x ) =| x |2+ i and φ2 ( x ) =| x |2– i are two independent solutions. The general
solution
(
φ = c1φ1 ( x ) + c2 φ2 ( x) = c1 | x |2+i + c2 | x |2–i = x 2 c1 | x |i +c2 | x |– i . )
(c) The indical equation q ( r ) = r ( r – 1) + r – 4π has roots 2 π and –2 π . Since both the
π π
roots are distinct φ1 ( x) =| x |2 and φ2 ( x) =| x |–2 are two solutions. The general
π π
solution of corresponding homogeneous equation is φ ( x) = c1 | x |2 + c2 | x |–2 .
We solve the non-homogeneous equation using the variation of constants method.
Case 1 : x > 0
π π π π
If x > 0 then | x |2 = x2 , | x | –2 = x –2
Let ψ ( x ) = u1 ( x ) φ1 ( x ) + u 2 ( x) φ2 ( x ) be a solution of given equation then
π π
x2 x –2 –2 π 2 π – 4 π
W(φ1 , φ2 ) = = – =
π –1 π –1 x x x
2 π x2 –2 π x –2
π π
0 x –2 π x2 0 π
W1 ( x) = = – x –2 , W2 ( x ) = = x2
–2 π –1 +2 π –1
1 –2 π x +2 π x 1
EXAMPLES
Q. 2. Classify the singular points in the finite plane for the equation
x 4 ( x 2 + 1) ( x –1) 2 y′′ + 4 x3 ( x –1) y′ + ( x + 1) y = 0
1. For each equation, locate and classify all its singular points in the finite plane
(a) x 2 y′′ + y = 0 (Ans.: x = 0 is regular, no irregular)
A second order equation with a regular singular point at x0 has the form
L( y ) = ( x – x0 ) 2 y ′′ + ( x – x0 ) a ( x) y′ + b ( x ) y = 0,
where a(x), b(x) are analytic functions at x0 i.e. they have power series expansions
∞ ∞
a( x) = ∑ α k ( x – x0 )k and b( x) = ∑ β k ( x – x0 )k
k=0 k=0
which are convergent on some interval |x – x0| < r0 for some r0 > 0.
Without loss of generality we assume x0 = 0. Then
L( y ) = x 2 y ′′ + x a ( x) y′ + b( x ) y = 0 and
∞ ∞
a( x) = ∑ α k xk , b( x) = ∑ β k xk which are convergent on an interval
k=0 k=0
| x | < r0 , r0 > 0. The Euler equation is a particular case of L(y) = 0 with a, b constants.
Differential Equations (112)
A second order equation with regular singular point has a power series solution. If functions
a (x), b (x) have power series expansion on some interval | x | < r0 then the power series solution
converges on the interval | x | < r0.
Theorem 3.2.1
Consider the equation
x 2 y′′ + a ( x ) xy′ + b ( x ) y = 0,
where a and b have convergent power series expansions for | x | < r0, r0 > 0. Let r1, r2
(Re r1 ³ Re r2 ) be the roots of the indicial polynomial
q ( r ) = r ( r –1) + a (0) r + b (0)
for 0 < | x | < r0 there is a solution f1 of the form
∞
φ1 ( x ) = | x |r1 ∑ ck x k (c0 = 1),
k =0
where the series converges for | x | < r0. If r1 – r2 is not zero or a positiove integer, there is a
second solution f2 for 0 < | x | < r0 of the form
∞ : :
φ2 ( x) = | x |r2 ∑ ck xk (c0 = 1),
k =0
where the series converge for | x | < r0.
:
The coefficients ck , ck can be obtained by substitution of the solution into the differential
equation.
Proof :
Suppose we have a solution f of the form
∞
φ ( x ) = x r ∑ ck x k (c0 ≠ 0, x > 0)
k =0
∞ ∞
b ( x ) φ ( x ) = ∑ β k x k x r ∑ ck x k
k =0 k =0
Differential Equations (113)
∞ : : k
= x r ∑ β k x k where β k = ∑ c j β k – j
k =0 j =0
∞ ∞
xa ( x ) φ ′( x) = x ∑ α k x k x r –1 ∑ (k + r ) ck x k
k =0 k =0
∞ ∞
= x r ∑ α k x k ∑ ( k + r ) ck x k
k =0 k =0
∞ : : k
= xr ∑ α k x k where α k = ∑ ( j + r ) c j α k – j
k =0 j =0
∞
x 2φ ′′( x) = x r ∑ (k + r ) (k + r –1) ck x k .
k =0
∞ ∞ : ∞ :
Thus, L(φ ) ( x) = x
r
∑ (k + r ) (k + r –1) ck xk + x r ∑ α k x k + x r ∑ β k x k
k =0 k =0 k =0
∞ : :
= xr ∑ (k + r ) ( k + r –1) ck + α k + β k x k
k =0
: :
L(f ) = 0 implies [ ]k = (k + r ) (k + r – 1) ck + α k + β k ) = 0
: : k = 0, 1, 2, 3,....
Using the definitions of α k , β k we can write [ ]k as
k k
[ ]k = (k + r ) (k + r – 1) ck + ∑ ( j + r) c j αk– j + ∑ c j βk– j
j =0 j =0
k –1
= [(k + r )(k + r –1) + (k + r ) α 0 + β 0 ]ck + ∑ ( j + r ) α k – j + β k – j c j
j =0
for k = 0 we must have
r ( r – 1) + r α 0 + β 0 = 0.
Since c0 ¹ 0 the second degree polynomial q given by
q (r ) = r (r – 1) + r α 0 + β 0
is called the indicial polynomial and the only admissible values of r are the roots of q.
[ ]k = q ( r + k )c k + d k = 0 (k = 1, 2,3,....) ......(3.2.1)
k –1
where dk = ∑ ( j + r )α k – j + β k – j c j (k = 1, 2,3,....) ......(3.2.2)
j =0
Note that dk is a linear combination of c0, c1, c2,..... ck–1 with coefficients involving the
known functions a, b and r. Leaving r and c0 indeterminant for the moment we solve equations
(3.2.1) and (3.2.2) successively in terms of c0 and r. The solutions we denote by Ck (r) and the
corresponding dk by Dk(r). Thus,
D (r)
D1 (r ) = (r α1 + β1 )c0 , C1 (r ) = – 1 ,
q (r + 1)
Differential Equations (114)
and in general
k –1 D (r )
Dk (r ) = ∑ ( j + r ) α k – j + β k – j C j (r ), Ck (r ) = – k (k = 1, 2, 3,...)
j =0 q (r + k )
The Ck thus, determined are rational functions of r, and the only points where they cease
to exist are the points r for which the denominator q(r + k) = 0 for some k = 1, 2, 3,.... Only two
such possible points exist.
Define F by
∞
Φ ( x, r ) = c0 x r + x r ∑ Ck ( r ) x k ....(3.2.3)
k =0
If the series converges for 0 < x < r0, then clearly
L(Φ ) ( x, r ) = c0 q (r ) x r ,
since Ck (r ) satisfies equation 3.2.1 for every k = 1, 2, 3....
∞
Thus if the function φ = x
r
∑ Ck x k is a solution of L(y) = 0 then r must be a root of the
k =0
indicial polynomial
q ( r ) = r (r – 1) + r α 0 + β 0
and ck (k ³ 1) are determined uniquely in terms of r and c0 given by equation (3.2.2), provided
q(r + k) ¹ 0 k = 1, 2, 3,... Conversely if r is a root of q and if Ck(r) can be determined then the
function f given by equation (3.2.3) is a solution of L(y) = 0 for any choice of c0, provided the
series in equation (3.2.3) is convergent.
Let r1, r2 be two roots of q and suppose Re r1 ³ Re r2. Then q(r1 + k) ¹ 0 for all k = 1, 2, 3,...
Thus, Ck(r1) exists for all k = 1, 2, 3,... and for c0 = C0 (r) = 1 we get a solution.
∞
φ1 ( x ) = Φ ( x, r1 ) = x r1 ∑ C k (r1 ) x k (C0 (r ) = 1),
k =0
is a solution of L(y) = 0, provided the series converges.
If r2 is a root of q distinct from r1 and q(r2 + k) ¹ 0 for k = 1, 2, 3,..., then clearly Ck(r2) is
defined for k = 1, 2, 3,.... and the function F2 defined by
∞
Φ 2 ( x) = Φ ( x, r2 ) = x r2 ∑ Ck ( r2 ) x k (C0 (r2 ) = 1)
k =0
is another solution of L(y) = 0, provided the series is convergent. The condition q(r2 + k) ¹ 0 for
k = 1, 2, ... is same as r2 + k ¹ r1 for any k = 1, 2, 3,.... or r1 – r2 ¹ k i.e. r1 – r2 is not a positive
integer and the result follows.
Illustration :
Consider the equation
3
L( y ) = x 2 y ′′ +
xy ′ + xy = 0
2
As per theorem 3.2.1 we assume the solution f of the equation L(y) = 0 as
3 3
= r (r – 1) + r c0 x r + (r + 1) (r ) + (r + 1) c1 + c0 x r +1
2 2
+ (r + 2) (r + 1) + (r + 2) c2 + c1 x r + 2 + ....
3
2
3
q( r ) = r ( r – 1) + r is the indicial polynomial
2
L(φ ) = q (r )c0 x r + [q (r + 1) c1 + c0 ] x r +1 + [q (r + 2) c2 + c1 ] x r + 2 + .....
∞
= q (r )c0 x r + x r ∑ [q (r + k ) ck + ck –1 ] x k
k =1
Theorem 3.2.2
Consider the equation
L( y ) = x 2 y′′ + a( x) xy′ + b( x ) y = 0,
where a, b have power series expansions which are convergent for | x | < r0, r0 > 0 . Let r1, r2
(Re r1 ³ Re r2) be the roots of the indicial polynomial
q ( r ) = r (r –1) + a (0) r + b(0).
If r1 = r2 there are two linearly independent solutions f1, f2 for 0 < | x | < r0 of the form
φ1 ( x ) = | x |r1 σ 1 ( x ), φ 2 ( x ) =| x |r1+1 σ 2 ( x ) + (log | x |) φ1 ( x ),
Differential Equations (117)
where s1, s2 have power series expansions which are convergent for | x | < r0 and s1(0) ¹ 0.
If r1 – r2 is a positive integer there are two linearly independent solutions f1, f2 for
0 < | x | < r0 of the form
φ1 ( x ) = | x |r1 σ 1 ( x),
φ2 ( x) = | x |r2 σ 2 ( x) + c(log | x |) φ1 ( x),
where s1, s2 have power series expansions which are convergent for | x |< r0 , σ 1(0) ≠ 0,
σ 2 (0) ≠ 0, and c is a constant. It may happen that c = 0.
Proof :
For x > 0, suppose we have a solution f of the form
∞
φ ( x ) = x r ∑ ck x k .
k =0
∞ : :
L(φ ) ( x ) = x r ∑ (k + r ) (k + r –1) ck + α k + β k x k
k =0
: k : k
where α k = ∑ ( j + r ) c j α k – j and β k = ∑ c j β k – j
j=0 j=0
∞ ∞
a( x) = ∑ α k x k , b( x) = ∑ β k x k .
k =0 k =0
L (φ ) ( x ) = 0 implies
: :
[ ]k = (k + r ) (k + r –1) ck + α k + β k = 0, k = 0,1, 2,3,...
k k
= (k + r ) (k + r – 1) ck + ∑ ( j + r ) c j α k – j + ∑ c j β k – j
j =0 j =0
= [(k + r ) (k + r – 1) ck + (k + r ) α 0 + β 0 ] ck
k –1
+ ∑ ( j + r ) α k – j + β k – j c j
j =0
For k = 0 we must have
q(r ) = r ( r – 1) + r α 0 + β 0 = 0
Then
[ ]k = q (r + k ) ck + d k = 0 .... 3.2.4
k –1
where dk = ∑ ( j + r ) α k – j + β k – j c j .... 3.2.5
j =0
Here, we are going to consider two cases according as the roots r1, r2 (Re r1 ³ Re r2) of the
indicial polynomial q(r) satisfy.
q(r ) = (r – r1 ) ( r – r2 )
and hence,
q ( r + m) = ( r + m – r1 ) ( r + m – r2 )
= (r + m – r2 – m) (r + m – r2 )
= (r – r2 ) (r + m – r2 )
If Dm(r) also has (r1 – r2) as a factor (i.e. Dm(r2) = 0 ), then it will get cancel from both the
sides of equation q(r + m) Cm(r) = – Dm(r) and would give Cm(r2) as a finite number. Then
Cm +1 (r2 ), Cm + 2 (r2 ),.....
all exist. In this special situation we will have a solution f2 of the form
∞
φ2 ( x) = x r2 ∑ Ck (r2 ) x k (C0 (r2 ) = 1)
k =0
L(ψ ) ( x, r ) = (r – r2 ) q (r ) x r
Therefore L(ψ ) ( x, r2 ) = 0 and
ψ ( x ) = ψ ( x, r2 )
is the second solution of L(y) = 0
Since C0 (r2 ) = C1 (r2 ) = ..... = Cm –1 (r2 ) = 0, the series y actually starts with the m-th power
in x.
To get a solution associated with r2 differentiate
L(ψ ) ( x, r ) = (r – r2 ) q (r ) x r
with respect to r then
∂ ∂ψ
[L(ψ ) ( x, r ) ] = L
( x, r )
∂r ∂r
= q(r ) x r + (r – r2 ) [q′( x) + (log x) q( r )] x r
∂ψ
and L = 0 at r = r2
∂r
∂ψ
and φ2 ( x) = ( x, r2 )
∂r
∞ ∞
= x r2 ∑ Ck′ (r2 ) x k + (log x) x r2 ∑ Ck +m (r2 ) x k +m
k =0 m=0
∞ ∞
= x r2 ∑ Ck′ (r2 ) x k + (log x) x r2 + k ∑ Ck +m (r – m) x m
k =0 m=0
∞ ∞
= x r2 ∑ Ck′ (r2 ) x k + (log x) x r1 C ∑ Cm (r1 ) x m
k =0 m=0
∞
φ2 ( x) = x r2 ∑ Ck′ (r2 ) x k + (log x) ⋅ c ⋅ φ1 ( x)
k =0
Where c is constant.
For x < 0, we replace x r1 , x r2 , log x everywhere by | x |r1 ,| x |r2 , log |x| respectively and the result
follows.
The method used in the theorem 3.2.2 is called the Frobenius method. The solutions f1, f2
are linearly independent. Thus, if the roots are equal or they differ by an integer then theorem
3.2.2 gives two linearly independent solutions of the differential equation
L( y ) = x 2 y′′ + x a( x) y′ + b( x) y = 0.
EXAMPLES
∞ ∞
then φ ′( x) = x r –1 ∑ (k + r ) ck x k and φ ′′( x ) = x
r –2
∑ (k + r ) (k + r –1) ck x k
k =0 k =0
d k +1 ( x ) –3 x 3| x|
= = → 0
d k ( x) + 1 + (3k + 1) ( k + 1)
3 k (k 1)
3
as k ® ¥ provided | x | < ¥. Thus, series defining f2 is convergent for all finite x.
Thus f1, f2 are solutions of the given equation.
Differential Equations (123)
(b) Suppose for x > 0 we have a solution f of the form
∞
φ ( x) = x r ∑ Ck x k , c0 ≠ 0
k =0
Q. 2. Obtain two linearly independent solutions of the following equations which are
valid near x = 0.
(a) x 2 y′′ + 3 x y′ + (1 + x ) y = 0
(b) x 2 y′′ + 2 x 2 y ′ – 2 y = 0
Ans. :
(a): For x > 0 suppose we have a solution f of the form
L(φ ) ( x, r ) = x r ∑ (k + r ) (k + r – 1) Ck x k + 3 x r ∑ (k + r ) Ck x k + ∑ Ck x k + r + x r ∑ Ck x k +1
∞ ∞
= ∑ [(k + r )(k + r –1) + 3(k + r ) + 1) ]Ck x k +r + ∑ ck x k +r +1
k =0 k =0
L (φ ) ( x, r ) = 0 implies
∞ ∞
∑ [(k + r )(k + r –1) + 3(k + r ) + 1]Ck x k + r + ∑ Ck x k + r +1 = 0
k =0 k =0
Since, q( r ) = ( r + 1) 2 , q( r + k ) = (r + k + 1) 2 and
– Ck –1
Ck =
(r + k + 1) 2
–1 –1 –1 –1 –1
=
2
2 2
2 ⋅⋅⋅ 2
(r + k + 1) (r + k ) (r + k – 1) (r + k – 2) (r + 2)
(–1) k C0
=
[(r + 2) (r + 3) ⋅⋅⋅ (r + k – 2) (r + k – 1) (r + k ) (r + k + 1)] 2
The first solution will be constructed by substituting Ck’s at r = – 1 in the series. Ck at
r = – 1 is
(–1)k C0
Ck =
k !2
(–1) k C0
Ck = .
[(r + 2) (r + 3) (r + 4) ⋅⋅⋅ (r + k + 1) ] 2
Define D = (r + 2) (r + 3) (r + 4) ..... (r + k + 1) then
(–2) 1 1 1 1
Ck' = (–1) k C0 + +
2 r + 2 r +3 r + 4
+ ⋅⋅⋅ +
D r + k + 1
(–2) (–1) k 1 1 1 1
Ck′ (–1) = + + + ⋅⋅⋅ +
k! 2
1 2 3 k
∞ (–2) (–1) k 1 1
1 + + + ⋅⋅⋅ + x k + (log x) φ1 ( x)
1
φ2 ( x) = x –1 ∑
k =0 k! 2
2 3 k
To obtain solution for x < 0 we replace x by | x |.
Thus, the two solutions are
∞ (–1) k x k
φ1 ( x ) = | x | –1 ∑
k =0 k !2
∞ (–2) (–1) k 1 1 1 k
φ2 ( x ) = | x |–1 ∑ 1 +
2 3+ + ⋅⋅⋅ + x
k =0 k !2 k
∞ (–1) k x k
+ (log | x |) | x |–1 ∑ .
k =0 k !2
Check that series in both the solutions converge.
(b): For x > 0 suppose we have a solution f of the form
∞
φ ( x) = x r ∑ Ck xk (C0 ≠ 0)
k =0
Let L( y ) = x 2 y ′′ + 2 x 2 y ′ – 2 y then
∞
L(φ ) ( x, r ) = ∑ (k + r ) (k + r –1) Ck x k + r + 2 (k + r )Ck x k + r +1 – 2 Ck x k + r
k =0
∞ ∞
= ∑ [(k + r ) (k + r – 1) – 2]Ck x k + r + ∑ 2 (k + r ) Ck x k + r +1
k =0 k =0
Differential Equations (127)
= [r (r –1) – 2]C0 x r + {[(r + 1)(r ) – 2]C1 + 2 r C0 } x r +1
(k 2 – 3 k ) ck + 2 (k – 2) ck –1 = 0
k (k – 3) ck + 2(k – 2) ck –1 = 0
2 (k – 2) ck –1
ck = – , k = 4,5, 6, 7,....
k (k – 3)
–2 (k – 2) –2 (k – 3) –2 (k – 4) –2 (2)
ck = ⋅⋅⋅ c3
k ( k – 3) ( k – 1) (k – 4) ( k – 2) ( k – 5) 4 ⋅1
2 ∞ (–2) k –3 ( k – 2) 6 k –1
= c0 ( x –1
– 1) + c3 x + ∑ x
k=4 k!
∞ (–2) k ( k + 1) 6
= c0 x –1 (1 – x ) + c3 x 2 + ∑ x k +2
k =1 (k + 3)!
∞ (–2) k ( k + 1) 6
= c0 x –1 (1– x) + c3 x 2 1 + ∑ xk
k = 1 (k + 3)!
Thus, we get two solutions
∞ (–2) k (k + 1) 6
φ1 ( x ) = x –1 (1– x) and φ2 ( x ) = x 2 1 + ∑ xk
k = 1 (k + 3)!
These are two solutions for x > 0 for x < 0 replace x by | x | we get,
∞ (–2) k ( k + 1) 6
φ1 ( x ) = | x |–1 (1– x) and φ 2 ( x ) =| x |2 1 + ∑ xk
k = 1 (k + 3)!
Check that series appearing in f2 is convergent series.
EXERCISE
1. Compute indicial polynomials and their roots for the following equations.
(a) x 2 y′′ + ( x + x 2 ) y ′ – y = 0
1
(b) x 2 y ′′ + x y′ + x 2 – y = 0
4
(c) 4 x 2 y′′ + (4 x 4 – 5 x) y′ + ( x 2 + 2) y = 0
(d) x 2 y′′ + ( x – 3x 2 ) y′ + e x y = 0
(b) x 2 y′′ + ( x – x 2 ) y′ + y = 0
3. For each equation obtain two linearly independent solutions valid near origin
(a) 2 x ( x – 1) y′′ + 3( x –1) y′ – y = 0
(b) 2 x y′′ + 5(1 + 2 x ) y′ + 5 y = 0
(c) 3 x y′′ + (2 – x) y′ – 2 y = 0
(d) 2 xy′′ + (1 – 2 x 2 ) y′ – 4 xy = 0
5. Obtain two linearly independent solutions of the following equations which are valid
near x = 0.
(a) x 2 y′′ – 2 x ( x + 1) y ′ + 2( x + 1) y = 0
Answers :
1. (a) q( r ) = r 2 – 1 ; r1 = 1, r2 = –1
2 1 1 1
(b) q( r ) = r – ; r1 = , r2 = –
4 2 2
2 9 1 1
(c) q( r ) = r – r+ ; r1 = 2, r2 =
4 2 4
(d) q(r ) = r 2 + 1 ; r1 = i, r2 = – i
∞ (–1) k x k 1 ∞ (–1)k x k 1 –x
2. (a) y1( x) = x ∑ , y 2 ( x) = x 2 ∑ = x 2e 2
k =0 1⋅ 3 ⋅ 5 ⋅ 7 ⋅⋅⋅⋅(2k + 1) k =0 2k k !
3
3(–5) x n – – 12
(b) y1 = 1 + ∑ ; y2 = x 2 – 10 x
n !(2n + 1) (2n + 3)
∞ n + 13 ∞
(3 n + 4) x ( n + 1) x n
(c) y1 = ∑ ; y2 = 1 + ∑
n =0 4 ⋅ 3n n ! n =1 2 ⋅ 5 ⋅ 8 ⋅ ⋅ ⋅ ⋅(3 n – 1)
∞
x 2 k +1 ∞ 2k x 2k
(d) y1 = ∑ ; y2 = 1 + ∑
n =0 2k k ! n =1 3 ⋅ 7 ⋅11 ⋅⋅⋅⋅(4k – 1)
1
4. (a) r1 = ; r2 = –2
2
3 1
(b) r1 = ; r2 =
2 2
5. (a) y1 ( x ) = x ; y2 ( x ) = x (e 2 x – 1)
∞
2k xk
(b) y1(x) = x y2(x) = x–1 1+ 3x + 6x2 –3 ∑ − 4x log| x |
2 2
;
k =4 (k –3)k !
1 1 ∞ 60 x n
(c) y1 = 1 + x + x2 ; y2 = x 5 + ∑
2 12 n = 6 ( n – 5)! n (n – 1) ( n – 2)
∞ 6(–2) n –3 x n+1
(d) y1 = x – 2 x + 2 x y2 = x 4 + ∑
2 3
;
n=4 n!
Let L( y ) = x 2 y ′′ + xy′ + ( x 2 – α 2 ) y.
L(φ ) ( x, r ) = ∑ (k + r ) (k + r –1) C k x k +r + ∑ (k + r ) C k x k +r + ( x 2 – α 2 ) ∑ C k x k +r
∞ ∞
= ∑ (k + r )(k + r –1) + (k + r ) – α 2 Ck x k + r + ∑ Ck x k + r + 2
k =0
k =0
= r ( r – 1) + r – α 2 C0 x r + ( r + 1) r + ( r + 1) – α 2 C1x r +1
{ }
∞
+ ∑ (k + r ) 2 – α 2 Ck + Ck –2 x k +r
k =2
The indicial equation is
q( r ) = r ( r –1) + r – α 2 = 0
q( r ) = r 2 – α 2 = 0.
The indicial polynomial q(r) has two roots r1 = a and r2 = –a . We shall construct solutions
for x > 0. We consider three cases namely a = 0, 2a is not a positive integer and 2a is a positive
integer.
Case 1 : a = 0
Since the roots are both equal to zero by theorem 3.2.2, there are two solutions f1, f2 of
the form
φ1 ( x) = σ 1 ( x) and φ2 ( x) = x σ 2 ( x) + (log x) φ1 ( x),
Where s1(x), s2(x) have power series expansions which converge for all finite x.
Since a = 0,
L( y ) = x 2 y′′ + x y′ + x 2 y
∞ ∞
Suppose σ 1 ( x ) = xα ∑ Ck x k = ∑ Ck x k (C0 ≠ 0) be a solution of L(y) = 0. Then
k =0 k =0
∞ ∞ ∞
L(σ 1 ) ( x ) = ∑ k (k – 1) Ck x k + ∑ k Ck x k + ∑ Ck x k +2
k =2 k =1 k =0
∞ ∞ ∞
= ∑ k (k –1) Ck x k + C1x + ∑ k Ck x k + ∑ Ck –2 x k
k =2 k =2 k =2
∞
= C1x + ∑ {[k (k – 1) + k ]Ck + Ck –2 } x k
k=2
L(φ 2 ) ( x) = x 2φ 2′′ + x φ 2′ + x 2φ 2 = 0
∞ φ ( x)
φ2′ ( x) = ∑ k ck x k –1 + 1 + (log x) φ1′ ( x)
k =1 x
(–1) m 2m
[2m (2m –1) + 2m]c2m + c2m–2 = –2 2m
2 (m!) 2
(–1) m+1 m
[4m 2 ]c2m + c2m –2 = , m = 1, 2,3, 4,.....
22m –2 (m!) 2
1 (–1) m+1 m
c2 m = 2m –2 – c2 m –2 , m = 1, 2,3, 4....
(2 m) 2 2 (m !) 2
c0 = 0
∞ ∞
+ ∑ ck x k +α + 2 – α 2 ∑ ck xk +α
k =0 k =0
∞ ∞
= ∑ (k + α ) (k + α – 1) + (k + α ) – α 2 ck x k +α + ∑ ck x k +α + 2
k =0 k =0
∞ ∞
= ∑ (k + α ) 2 – α 2 ck x k +α + ∑ ck x k +α +2
k =0 k =0
∞
= 0 ⋅ c0 xα + (α + 1) 2 – α 2 c1 xα +1 + ∑ (k + α ) 2 – α 2 ck x k +α
k =2
∞
+ ∑ ck –2 x k +α
k =2
– c4 – c4 – c4 – c0
c6 = = = = 6
6(2α + 6) 12(α + 3) 2 ⋅ 3!(α + 3) 2 ⋅ 3!(α + 1) (α + 2) (α + 3)
– c6 – c0
c8 = = 8
8(2α + 8) 2 4!(α + 1) (α + 2) (α + 3)(α + 4)
In general,
(–1) m c0
c2m = .
22m m!(α + 1)(α + 2) (α + 3) ⋅⋅⋅ (α + m)
Thus the solution f1 becomes
∞
φ1 ( x) = xα ∑ ck xk (c0 ≠ 0)
k =0
∞ (–1) m x 2 m
= c0 xα + xα ∑
m =1 22 m m !(α + 1) (α + 2) ⋅⋅⋅ (α + m)
for a = 0 and c0 = 1, f1(x) becomes J0 (x). Before going for the second solution let us define
gamma function and study some properties of gamma function.
z e– x T T – x
= TLim
→∞ x | – ∫ – e ⋅ ze z –1dx
–1 x = 0 0
T
= TLim – x z –1
→∞ z ∫ e x dx
0
∞
= z ∫ e – x x z –1dx
0
= z Γ( z )
T
Observe that Lim
T →∞ x ze – x | = Lim
T →∞ T z ⋅ e –T – 0 = 0
0
By definition 3.3.1
∞
Γ(1) = ∫ e – x dx = 1
0
Thus, if z is a positive integer n,
Γ(n + 1) = n Γ(n) = n (n – 1) (n – 2)
= ( n ) ( n –1) ( n – 2) ( n – 3) ⋅⋅⋅⋅Γ(1)
=n!
Thus, gamma function is an extension of the factorial function to numbers which are not
integers.
Suppose Re z < 0 and z is not a negative integer then there is a natural number N such that
– N < Re z < – N + 1
But then Re (z + N) > 0 and therefore we can define
Γ( z + N ) = ( z + N –1) ( z + N – 2) ⋅⋅⋅⋅ ( z + 1) z Γ( z ). Then
Γ( z + N )
Γ( z ) = , (Re z < 0)
( z + N –1) ( z + N – 2) ( z + N – 3) ⋅⋅⋅⋅ ( z + 1) z
The gamma function is not defined at 0, –1, –2, –3,.......
We have a solution f1 (x) as
Observe that this formula for Ja reduces to J0 when a = 0 as G(m + 1) = m!. Ja (x) is one
solution of Bessel equation with a ¹ 0 and Re a ³ 0.
To determine second solution we have to consider two situations. Either 2a is not a positive
integer or 2a is a positive integer. We determine second solution for both the situations.
On repeating the same calculations we have carried out for the root a, (replace a by – a
everywhere)
We get the second solution
–α ∞ 2m
x (–1) m x
J –α ( x ) = ∑ .
2 m = 0 m ! Γ ( m – α + 1) 2
Observe that G (m – a + 1) exists for all m = 0, 1, 2, 3,..... since a is not a positive integer.
∞ ∞
or ∑ (k – n) 2 – n 2 ck x k – n + ∑ ck x k – n +2 + 2 c x J n′ ( x) = 0
k =0 k =0
∞ ∞
that is 0 ⋅ c0 x – n + (1– n)2 – n 2 c x1– n + ∑ (k – n)2 – n 2 ck x k – n + ∑ ck x k – n +2
k =2 k =0
+2 c x J n′ ( x) = 0.
∞
Since ∑ ck x k – n+ 2 = c0 x – n+ 2 + c1 x1– n+ 2 + c2 x 2– n+ 2 +c3 x3– n+ 2 + ⋅⋅⋅⋅
k =0
∞
= ∑ ck –2 x k – n , we get
k =2
{ }
∞
0 ⋅ c0 x – n + (1– n)2 – n2 c1 x1– n + ∑ (k – n)2 – n 2 ck + ck –2 x k – n + c x J n′ = 0.
k =2
On multiplying by xn we have
∞
0 ⋅ c0 + (1– 2n) c1 x + ∑ [k (k – 2n) ck + ck –2 ] xk = –2 c x J n′ ( x) ⋅ x n .............. (3.3.1)
k =2
∞ ∞ (2m + n ) (–1)m x 2m + 2 n
0 ⋅ c0 + (1– 2n) c1 x + ∑ [k (k – 2n) cn + ck –2 ] x k = – c ∑
k =2 m=0 m ! (m + n)! 22m + n –1
.............. (3.3.2)
2n
The series on the right side begin with x and since n is positive integer, the right side do
not contain any odd terms. Therefore c1 = 0 c2k+1 = 0 for k = 1, 2, 3,.... and if n > 1 then
k ( k – 2n) ck + ck –2 = 0, for k = 2,3, 4,...., 2n – 1.
c0
in general c2 j = 2j
, j = 1, 2,3,....., n –1 ...(3.3.3)
2 j !(n – 1) (n – 2) (n – 3) ⋅⋅⋅⋅ (n – j )
In particular
c0
c2n –2 =
22n –2 (n – 1)!(n –1)!
On comparing the coefficients of x2n in equation (3.3.2) we get
n c
c2n –2 = – c n –1
=–
n !2 (n –1)!2n –1
c0 c
Thus c2n –2 = =–
22n –2 (n – 1)!2 2n –1 (n –1)!
c0
and therefore c=– n –1
. .....( 3.3.4 )
2 (n –1)!
Since c2 n–2 is used to find c, c2n remains undetermined, but the remaining coefficients
c2n+2, c2n + 4, c2n + 6,..... can be obtained by comparing the coefficients of x2(n + j) in equation
(3.3.2).
(–1) j (2 j + n) 1
(2n + 2 j ) (2n + 2 j – 2n) c2n+ 2 j + c2n+ 2 j –2 = – c j = 1, 2,3.....
j !(n + j )! 2 +n –1
2 j
for j = 1 we have
(–1) (n + 2)
4(n + 1) c2 n+2 + c2 n = – c
(n + 1)! 2n+1
c 1 c2n
c2n+ 2 = n+1 1 + –
4⋅2 ⋅ (n + 1)! n + 1 4(n + 1)
c2n
=
–c + 1 + 1 + 1 + ⋅⋅⋅⋅ + 1
Choose 1
4(n + 1) 4 ⋅ 2 n+1 ⋅ (n + 1)! 2 3 4 n
–c 1 1 1 1
i.e. c2 n = n +1 1 + + + + ⋅⋅⋅⋅ +
2 n! 2 3 4 n
+c + + 1 + 1 + 1 + ⋅⋅⋅⋅ + 1 + 1
Then c2n+ 2 = n+1 1 1
4 ⋅ 2 (n + 1)! 2 3 4 n n +1
for j = 2 we have
c (–1) 2 ( n + 4)
4 ⋅ 2 (n + 2) c2 n + 4 = – – c2 n+2
2!(n + 2)!2n +3
–c 2 1 c 1 1 1
c2 n + 4 = n+ 3
1+ – + 1 + 1 + + + ⋅⋅⋅⋅ +
4 ⋅ 2 ⋅ 2!(n + 2)! 2 n + 2 4 ⋅ 2( n + 2) n 1
4 ⋅ 2 (n + 1)! 2 3 n +1
=
–c 1 + 1 + + + 1 + 1 + ⋅⋅⋅⋅ + 1 + 1
n+2 2 n + 2
1 1
2
4 (n + 2)! 2 2 3 n n +1
–c 1 1 1 1 1
= 1 + 2 + 1 + 2 + 3 + ⋅⋅⋅⋅ + n + n + 1
4 (n + 2)! 2n + 2
2
– (–1) 2 c 1 1 1 1 1
= 1 + 2 + 1 + 2 + 3 + ⋅⋅⋅⋅ + n + n + 1
2 ⋅ 2!(n + 2)! 2n+4
It can be shown by induction that
– (–1) m c 1 1 1 1 1 1
c2n + 2m = n+2 m 1 + 2 + 3 + ⋅⋅⋅⋅ + m + 1 + 2 + 3 + ⋅⋅⋅⋅ + n + m ,
2 ⋅ m!(n + m)! 2
m = 1, 2,3,....
Finally we get a solution f2
∞
φ2 ( x ) = x – n ∑ ck x k + (c log x ) J n ( x )
k =0
2n –1 ∞
= x – n ∑ ck xk + x – nc2n –1x2n−1 + x – nc2n x2 n + x – n ∑ ck xk + (c log x) J n ( x)
k =0 k = 2n–1
Since all odd terms c2k+1 = 0, k = 1, 2, 3,....., we get,
EXAMPLES
1
Ex. 1. Suppose f is any solution of x y′′ + xy ′ + x y = 0 for x > 0 and let ψ ( x ) = x 2 φ ( x ) .
2 2
1
y ′′ + 1 + 2 y = 0, ( x > 0)
4x
The function satisfies f ( x ) = sin x satisfies y ′′ + y = 0
1
Since 1 + > 1 and sin x = 0 has infinitely many zeroes x = n π , n = 0,1, 2,3,....
4x 2
1 1
∫ φλ ( x) φ µ ( x) d x = ∫ x J 0 (λ x) J 0 ( µ x) dx = 0
0 0
1
Ans. (b) : φ µ′′ ( x) + φ ( x) = – µ 2φµ
2 µ
4x
1
φλ ′′ ( x ) + 2 φλ ( x ) = – λ 2φλ .
4x
Differential Equations (145)
Multiply first equation by fl and second equation by fm and subtract these equations.
φ µ ′′ φλ – φλ ′′ φ µ = – µ 2φ µ φλ + λ 2φ µ φλ = (λ 2 – µ 2 )φλ φ µ .
Thus, (
(λ 2 – µ 2 )φλ φ µ = φ µ ′ φλ – φλ ′ φ µ )′
Integrate above equation between 0 to 1. Since
1 1
φλ ( x) = x2J 0 (λ x), φλ (0) = 0 and φ µ ( x) = x2J 0 (µ x), φ µ (0) = 0.
Therefore
1
(λ 2 – µ 2 ) ∫ φλ ( x) φ µ ( x) d x = φ µ ′ (1) φ λ (1) – φλ ′ (1) φ µ (1)
0
1
Ans. (c) : (λ 2 – µ 2 ) ∫ φλ ( x) φ µ ( x) d x = φ µ ′ (1) φ λ (1) – φλ ′ (1) φ µ (1)
0
1 1 1 1
i.e. ∫ φλ ( x) φ µ ( x) dx = ∫ x2J 0 (λ x) x 2 J 0 (µ x)d x = 0
0 0
1
i.e. ∫ x J 0 (λ x) J 0 ( µ x) dx = 0.
0
∞ 2m
(–1) m x
Ans. : J 0 ( x) = ∑ 2
m = 0 m! 2
∞ (–1) m 2m x 2 m –1
J 0′ ( x ) = ∑ 2
m = 1 m! 2 2m
∞ (–1) m m x 2m –1
= ∑ 2
m = 1 m! 2 2 m –1
∞ (–1) m +1 (m + 1) x 2 m +1
= ∑ (Replace m by m + 1)
m = 0 ( m + 1)!
2
2 2m +1
2m
x ∞ (–1) m x
= – ∑
2 m = 0 m !( m + 1)! 2
= – J1 ( x)
–n ∞ 2m
x (–1)m x
Ans. : J – n ( x) = ∑
2 m = 0 m! (m – n + 1) 2
–n ∞ 2m
(–1)m x
=
x 1
∑ (As = 0 for k ≤ 0)
2 m = n m! (m – n)! 2 (k )
2 m+ 2 m
(–1) m+ n x
–n ∞
=
x
∑ (Replace m by m + n)
2 m = 0 ( m + n)! m ! 2
∞ n 2m
(–1) m
= (–1) n ∑ x
x
2 m = 0 m ! ( m – n)! 2
= (–1) n J n ( x ).
(b) ( x –α Jα )′ ( x ) = – x –α Jα +1 ( x )
α ∞ 2m
x (–1)m x
Ans. (a) : Jα ( x) = ∑
2 m = 0 m! ( m + α + 1) 2
x 2α ∞
2m
(–1)m
\ xα Jα ( x) =
x
α ∑
2 m = 0 m! (m + α + 1) 2
1 ∞ (–1)m x 2 m+ 2 α
= ∑
2α m = 0 m! (m + α + 1) 22m
= xα Jα –1 ( x )
α 2m
x (–1) m x
Ans. (b) : Jα ( x) = ∑
2 m! (m + α + 1) 2
x ∞ (–1)m x2m
=– ∑
2α +1 m = 0 m! (m + α + 2) 2m
α +1 ∞ 2m
x (–1)m x
But Jα +1 ( x) = ∑
2 m=0 m! (m + α + 2) 2
∞ 2m
–α x (–1)m x
\ –x Jα +1 = – α +1 ∑
2 m = 0 m! ( m + α + 2) 2
(
Thus x –α Jα )′ = – x –α Jα +1
8. Show that
(a) Jα –1 ( x ) – Jα +1 ( x) = 2 Jα ′ ( x)
(b) Jα –1 ( x ) + Jα +1( x) = 2 α x –1Jα ( x)
(a) Multiply equation (i) by x–a and equation (ii) by xa and add.
Jα –1 ( x ) – Jα +1 ( x ) = (α x –1 – α x –1 ) Jα + Jα ′ + Jα ′ = 2 Jα ′ ( x )
Thus, Jα –1 ( x) – Jα +1 ( x) = 2 Jα ′ ( x)
(b) Multiply equation (i) by x–a and equation (ii) by xa and subtract.
Jα –1 ( x ) + Jα +1( x) = 2 α x –1Jα ( x)
Thus, using results of example 7 we have proved the required result.
∞ 2m
(–1) m + 1 + 1 + ⋅⋅⋅⋅ + 1 x
Ans. : K 0 ( x ) = – ∑ 1 + log x J 0 ( x )
m =1 m! 2 2 3 m 2
1 m x 2m –1 1 ∞ (–1) m x
∞ 2m
(–1) m 1 1
=– ∑ 1 + + + ⋅⋅⋅⋅ + + ∑ + log x J 0′
m =1 m! 2 2 3 m 2 2 m –1 x m = 0 m !2 2
x ∞ (–1) m
+ 1 + 1 + ⋅⋅⋅⋅ + 1 m x
2 m –1
1 x 2 ∞ (–1) m x 2 m –1
= – ∑ 1 + 1 – + ∑
2 m = 2 m! 2 2 3 m 2 2 m –1 x 4 m = 2 m !2 2 2 m
+ log x J 0′ ( x)
∞ (–1) m x 2 m –1 1 1
1 x 1 1
= + + ∑ – m 1 + + + ⋅⋅⋅⋅ + + log x J 0′ ( x)
2
x 4 m = 2 m! 2 2 m –1
2 2 3 m
1 x x ∞ (–1) m x 2m –2 1 1 1 1
= + + ∑ – – m 1 + + + ⋅⋅⋅⋅ + + log x J 0′
x 4 2 m = 2 m! 2
2 2 m –2
2 2 3 m – 1
1 x x ∞ (–1) m+1 x
2m
1 (m + 1) 1 1 1
K 0′ ( x) = + + ∑ – – 1 + + + ⋅⋅⋅⋅ +
x 4 2 m = 1 (m + 1)!m ! 2 2 (m + 1) (m + 1) 2 3 m
– log x ⋅ J1 ( x )
2m
1 x 1 x (–1) m x 1 1 1 1
= + + ∑ m + 1 + 2 1 + 2 + 3 + ⋅⋅⋅⋅ + m
x 4 2 2 m !( m + 1)! 2
– log x ⋅ J1 ( x )
= – K1 ( x)
Thus, K 0′ ( x ) = – K1 ( x )
EXAMPLES
1
4. For a fixed a > 0 and l > 0 let φλ ( x) = x 2 Jα (λ x) show that
1 α2
4 –
′′
φλ + 2 φλ = – λ 2 φλ
x
5. If l, m are positive show that
1
(λ 2 – µ 2 ) ∫ φλ ( x ) φ µ ( x )d x = φ λ (1) φ µ′ (1) – φ µ (1) φλ′ (1)
0
At the beginning of chapter 3 we have defined singular points of linear differential equation
of order n on the domain | x | < ¥. In unit 2 of chapter 3 we have discussed the power series
solutions of second order differential equation with regular singular points. These singular points
lie in a finite plane | x | < ¥. Often it is necessary to investigate solution of the differential
equation for large values of | x |. A simple way of doing this is to change the independent variable
1
by its reciprocal x = and study the solution of the resulting equation near t = 0. If the resulting
t
equation possesses the regular singular point of t = 0. We say that the original equation has a
regular singular point at infinity. The results on analytic solution and equations with regular
singular point at t = 0 can be applied to the transformed equation. Analysis of equation at t = 0
gives the analysis of given equation for infinite x.
Let us consider the second order differential equation
L( y ) = y ′′ + a1 ( x) y′ + a2 ( x) y = 0
for large values of | x |.
Suppose f is a solution of L(y) = 0 for | x | > r0 for some r0 > 0.
1
and let φ% (t ) = φ , a%1 (t ) = a1 , a% 2 (t ) = a2 . These functions will
1 1 1
Define t =
x t t t
1
exist for | t |< .
x0
dφ% (t ) d dx –1 d φ%
= φ ( x) ⋅ = φ ′( x) 2 ⇒ φ ′( x) = – t 2 .
dt dx dt t dt
d 2φ% (t ) 1 d dx d –1
2
= – 2 (φ ′ ( x) ) + φ ′( x) 2
dt t dx dt dt t
1 ′′ 2
= 4
φ ( x) + 3 φ ′ ( x)
t t
Therefore
EXAMPLES
t r , r >/ 1 and denominator of q contain a factor t r , r >/ 2, t = 0 is a regular singular point. Thus,
infinity is a regular singular point of the given differential equation.
Ex. 2. Show that infinity is not a regular singular point for the equation
y ′′ + a y ′ + b y = 0
where a, b are constants, not both zero.
Ans. : y ′( x) = – t 2 y& (t )
y ′′( x ) = t 4 &&
y (t ) + 2 t 3 y& (t )
and y ′′( x ) + a y ′( x ) + b y ( x ) = 0 gives
t 4 &&
y (t ) + 2 t 3 y& (t ) – a t 2 y& (t ) + b y (t ) = 0.
Therefore
2t – a b
y (t ) + 2 y& (t ) + 4 y (t ) = 0
&&
t t
2t – a b
Here p (t ) = 2 and q (t ) = 4
t t
2
[If a = b = 0 then p (t ) =and q (t) = 0. Since denominator of p(t) contain a factor t r , r >/ 1
t
and denominator of q(t) contains a factor t r , r = 0 >/ 2 t = 0 is a regular singular point and
infinity is regular singular point of the equation.]
2t – a
Since either a or b is non-zero, p (t ) = 2
contains the determinator tr with r = 2 > 1 or
t
b
q (t ) = 4
contains the denominator tr with r = 4 > 2. Therefore t = 0 is an irregular singular
t
point of the transformed equation and infinity is an irregular singular point of the given equation.
1
Ans. : y ′( x ) = – t 2 y& (t ), y ′′( x) = t 4 &&
y (t ) + 2t 3 y& (t ), x =
t
\ Legendre equation becomes,
1 4
1 – 2 t &&
t
1
( )
y + 2 t 3 y& – 2 –t 2 y& + α (α + 1) y = 0
t
(t 2 – 1) t 2 &&
y + 2(t 2 – 1) t y& + 2 t y& + α (α + 1) y = 0
2t α (α + 1)
y+
&&
2
y& + 2 2 y=0
(t –1) t (t – 1)
2t
Here p (t ) = 2
contains a factor t r in the denominator with r = 0 >/ 1 and
(t – 1)
α (α + 1)
q (t ) = 2 2
contains a factor tr in the denominator with r = 2 >/ 2.
t (t –1)
Therefore by definition 3.1.3(b) t = 0 is a regular singular point of the transformed equation
and infinity is a regular singular point of a given equation.
1
Ans. : Put x = then
t
1 4
1– 2 t &&
t
1
(
y + 2 t 3 y& – 2 –t 2 y& + 2 y = 0
t
)
or (t 2 –1) t 2 &&
y + 2 t y& + 2 t y& + 2 y = 0
L( y ) = t 2 (t 2 – 1) &&
y + 2 t 3 y& + 2 y = 0
From example 2 we observe that x = 0 is a regular singular point of L(y) = 0.
Let f be a solution of L(y) = 0 of the form
∞
φ (t ) = t r ∑ ck t k
k =0
∞ ∞
= ∑ (k + r ) (k + r – 1) ck t k +r + 2 – ∑ (k + r ) (k + r – 1) ck t k +r
k =0 k =0
∞ ∞
+2 ∑ (k + r ) t k + r + 2 + 2 ∑ ck t k + r
k =1 k =0
∞
= ∑ [(k + r )(k + r –1) + 2(k + r )] ck t k + r + 2
k =0
∞
– ∑ [(k + r ) (k + r – 1) – 2 ]ck t k +r
k =0
∞
= ∑ [(k + r ) (k + r + 1)] ck t k +r + 2 – (r (r –1) – 2) c0 t r – ((r + 1) (r ) – 2 ) c1t r +1
k =0
∞
– ∑ [(k + r ) (k + r – 1) – 2] ck t k +r
k =2
∞
= – r 2 – r – 2 c0 t r – r 2 + r – 2 c1 t r +1 – ∑ [(k + r ) (k + r – 1) – 2] ck t k +r
k =2
∞
+ ∑ [(k + r – 2) (k + r –1) ck –2 ]t k + r
k =2
( ) (
= – r 2 – r – 2 c0 t r – r 2 + r – 2 c1 t r +1)
∞
–∑ {[(k + r) (k + r –1) – 2] ck – [(k + r – 2) (k + r –1)]ck –2}t k +r
k =2
In particular
(2k – 2) (2 k – 1)
c2k +1 = c2k –1 k = 2,3, 4,....
2k (2k –1) – 2
(2k – 2) (2k – 1)
= c2 k –1
(2k – 2) (2k + 1)
2k –1
= c2 k –1 , k = 2,3, 4,.....
2k + 1
2 k –1 2 k – 3 2 k – 5 2 k – 7 1
c2 k +1 = ⋅⋅⋅⋅⋅ c3
2k + 1 2k – 1 2k – 3 2k – 5 3
c3
=
2k + 1
Thus we get a solution
φ (t ) = c0t –1 + c3t 2 + c5t 4 + c7 t 6 + ....
∞
= c0t –1 + ∑ c2 k +1 t 2 k
k =1
∞ c3
= c0t –1 + ∑ t 2k
k =1 (2 k + 1)
∞
t 2k
= c0t –1 + c3 ∑
k =1 (2k + 1)
Thus, we get two solutions
∞ t 2k
φ1 (t ) = t –1 and φ2 (t ) = ∑
k =1 2k + 1
∞
Let φ2 (t ) = ∑ d k (t ).
k =1
By ratio test
d k +1 (t ) t 2 k + 2 2k + 1 (2k + 1) t 2
= × 2k =
d k (t ) 2k + 3 t (2k + 3)
2k + 1
Since Lt
k →∞ =1
2k + 3
The series converger for | t | < 1.
1
But x = , therefore
t
(b) x 2 ( x 2 – 4) y′′ + 2 x 3 y ′ + 3 y = 0
(c) y′′ + xy = 0
(d) x 2 ( x – 4) 2 y′′ + 3 x y′ – ( x – 4) y = 0
Ans. : In chapter 3 Unit I, example 3, we have classified all its singular points in a finite plane.
It remains to check whether infinity is a singular point and whether it is a regular singular
point.
(a) a0 ( x) = x 3 ( x –1) = 0 gives x = 0, x = 1 are singularities x = 0 is irregular singular point
whereas x = 1 is a regular singular point.
put x=
1
t
then
dy dt
y ′( x ) =
dt dx
⋅ = y& (t ) (–t 2 ) and y ′′( x ) =
dt
(
d &
)
y (t ) (–t 2 )
dt
dx
= –t &&
2
y (t ) – 2 t y& (t ) –t so y ′′( x ) = t 4 &&
2
y + 2 t 3 y&.
L( y ) = x 3 ( x – 1) y′′ + ( x – 1) y′ + 4 x y
1 1 4 1 1
= 3t
y + 2t 3 y& + –1 (–t 2 ) y& + 4 y
– 1 t &&
t t t
y + – 2 – t + t 2 y& + y
2 4
= (1 – t ) &&
t t
2 – (2 + t ) t + t 3 4
L ( y ) = (1 – t ) &&
y+ y& + y
t t
L(y) = 0 can be put in the form
(t – 1) (t 2 – 2) 4
y+
&& y& + y =0
(1– t ) t t (1 – t )
2
t –2 4
or y+
&& y& + y=0
t t (1– t )
t2 – 2 4
This equation is of the type y ′′ + p (t ) y ′ + Q (t ) y = 0 where p (t ) = and Q (t ) =
t t (1– t )
Since the denominator of p(t) contains a factor t r , for r = 1 >/ 1 and q(t) contains a
r
denominator t , for r = 1 >/ 2, t = 0 is a regular singular point.
Thus, x = 1 and infinity are regular singular whereas x = 0 is irregular singular point.
L( y) = x 2 ( x 2 – 4) y′′ + 2 x3 y′ + 3 y
11 2
= 2 2
– 4 t 4 &&
y + 2t 3 y& + 3 (– t 2 y& ) + 3 y
t t t
2 2
= (1 – 4t 2 ) &&
y + – 8t – y& + 3 y
t t
= (1– 4t 2 ) &&
y – 8 t y& + 3 y
2
y – 8t y& + 3 y = 0, infinity is not a singular
Since t = 0 is not a singular point of (1– 4 t ) &&
point of the given equation.
1
(c) y ′′ + x y = (t 4 &&
y + 2t 3 y& ) + y.
t
4 1
y + 2t 3 y& + y = 0 can be written in the form &&
Therefore L( y ) = t && y + p (t ) y& + q (t ) y = 0
t
2 1
where p (t ) = and q (t ) = 5 .
t t
Here t = 0 is a singular point but since the denominator of q(t) contains a factor
tr, r = 5 > 2, t = 0 is not a regular singular point.
Since t = 0 is irregular singular point infinity is irregular singular point of the equation
y ′′ + x y = 0.
(d) x = 0 is regular singular point and x = 4 is irregular singular point.
1
Put x = then y′ = –t 2 y& , y′′ = t 4 &&
y + 2 t 3 y&
t
L ( y ) = x 2 ( x – 4) 2 y′′ + x y ′ – ( x – 4) y
2
1 1 4 1 1
= 2 t
y + 2t 3 y& + [–t 2 y& ] – – 4 y
– 4 t &&
t t t
2(4t –1) 2
– t y& +
4t –1
= (4t –1) 2 &&
y+ y
t t
31t 2 – 16t + 1 4t – 1
= (4t – 1)2 &&
y+ y& + y
t t
31t 2 –16t + 1
y + p (t ) y& + q (t ) y = 0 where p (t ) =
L(y) = 0 can be written in the form && and
(4t –1) 2 ⋅ t
Contents :
Unit 1 : The method of successive approximations.
Introduction :
In the last three chapters we have seen the methods of finding a solution to the given linear
differential equations. For linear differential equation with constant coefficients there is a method
to find all the solutions whereas for linear equations with variable coefficients, there are very
few types of equations whose solutions can be expressed in terms of elementary functions and
therefore we go for power series solutions. All the equations considered so far were linear
differential equations.
In this chapter we consider the general first order equation y ′ = f ( x, y ) where f is some
continuous function (need not be linear in y) Only in special cases it is possible to find explicit
analytic expressions for the solutions of y ′ = f ( x, y ) .
Our main purpose in this chapter is to prove that a wide class of initial value problems
y ′ = f ( x, y ), y ( x0 ) = y0
has a solution. Though it may not be possible to find out the exact solution, it is feasible to
construct a sequence of approximate solutions that may converge to the exact solution.
Differentiate this equation with respect to x and use the fundamental theorem of integral
calculus. The integral equation becomes
φ ′( x ) = f ( x, φ ( x )) for all x ∈ I .
From (4.1.2) it is obvious that φ ( x0 ) = y0 .
Thus f is a solution of equation (4.1.1).
Successive approximate solutions
As a first approximation to a solution defined
φ0 ( x ) = y0 .
Then f 0 satisfies an initial condition but does not in general satisfy the differential equation.
Since f 0 is a first approximate solution, substitute y = f0 in equation (4.1.2) to generate
second approximate solution. Call this solution as f1 then
x
φ1 ( x) = y0 + ∫ f (t , φ0 (t )) dt.
x0
∞
We get a sequence of functions {φ k }k = 0 . If this sequence converges then it may happen
that the limit function will turn out to be the solution of differential equation (4.1.1).
We now show that there is an interval I containing x0 where all the functions φk , k = 1, 2,.....
exist. Since f is continuous on a compact set R, it is bounded on R, that is there exists a constant
M > 0 such that
| f ( x, y ) | ≤ M for all ( x, y ) ∈ R.
Theorem 4.1.2
The successive approximations defined by (4.1.3) exist and are continuous on
I = {x / | x – x0 | ≤ α where α = min {a, b / M }}
and for x ∈ I , ( x, φ k ( x)) ∈ R.
The function fk satisfy
| φ k ( x) – y0 |≤ M | x – x0 | for all x in I
Proof : We will prove this result by mathematical induction,
(i) Clearly φ0 ( x ) = y0 is continuous on I and
| φ0 ( x) – y0 | = 0
Thus the theorem is true for k = 0.
x
(ii) φ1 ( x) = y0 + ∫ f (t , φ0 (t )) d t
x0
x
= y0 + ∫ f (t , y0 ) dt
x0
Therefore | φ1 ( x) – y0 | ≤ M | x – x0 |
Since f is continuous on R the function F0 defined by
F0 (t ) = f (t , y0 )
is continuous on I. Therefore f1 defined by
(Since Fk (t ) = f (t , φ k ) ≤ M )
Thus fk+1 exist is continuous and satisfies the required inequality.
Definition : Let f be a function defined for (x, y) in a set S. We say f satisfies a Lipschitz
condition on S if there exists a constant K > 0 such that
| f ( x, y1 ) – f ( x, y2 ) | ≤ K | y1 – y2 |
for all ( x, y1 ), ( x, y2 ) ∈ S . The constant K is called Lipschitz constant.
Theorem 4.1.3 :
Suppose S is either a rectangle
| x – x0 | ≤ a, | y – y0 | ≤ b (a, b > 0);
or a strip
| x – x0 | ≤ α , | y | < ∞ (a > 0)
and that f is real valued function defined on S.
∂f
Such that exists, is continuous on S and
∂y
∂f
( x, y ) ≤ K , for ( x, y ) ∈ S and for some K > 0. Then f satisfies a
∂y
Lipschitz condition on S with Lipschitz constant K.
y1
∂f
Proof : f ( x, y1 ) – f ( x, y2 ) = ∫ ( x, t ) d t
y2 ∂ y
y1
∂f
≤ ∫ ∂ y ( x, t ) d t
y2
y1
≤ ∫ K dt
y2
≤ K y1 – y2
Thus, f ( x, y1 ) – f ( x, y2 ) ≤ K | y1 – y2 | for all (x, y1) (x, y2) in S.
EXAMPLES
Here f (t , y ) = 3 y + 1 . Therefore
x
φ1 ( x) = 2 + ∫ [3 φ0 + 1] d t
0
x
= 2 + ∫ 7dt = 2 + 7 x
0
x
φ2 ( x) = 2 + ∫ [3 φ1 + 1] d t
0
x
= 2 + ∫ [3 (2 + 7 t ) + 1] d t
0
x
= 2 + ∫ (21t + 7) d t
0
21x 2 21x 2
= 2+ + 7x = 2 + 7x + .
2 2
x
φ3 ( x) = 2 + ∫ [3 φ2 (t ) + 1] d t
0
x
21
= 2 + ∫ 3 2 + 7 t + t 2 + 1 d t
0 2
x
63
= 2 + ∫ 7 + 21t + t 2 d t
0 2
x 2 63 x 3
= 2 + 7 x + 21 +
2 2 3
21 21
= 2 + 7 x + x 2 + x3 .
2 2
(c) y′ – 3 y = 1
y = e3 x ∫ 1e –3 x d x + c1
e –3x
= e3 x + c1
–3
1
= – + c1 e3 x
3
2. For each of the following problems compute the first four successive approximations
f0, f1, f2 , f3 .
(c) y ′ = y 2 , y (0) = 1
Answers :
(a) φ0 ( x ) = y0 = 0, f ( x, y ) = x 2 + y 2
x
φ1 ( x) = y0 + ∫ f (t , φ0 (t ) ) d t
x0
x
(
= 0 + ∫ t 2 + 02 d t
0
)
x3
=
3
x
φ2 ( x) = y0 + ∫ f (t , φ1 (t ) ) d t
x0
x t3
= 0 + ∫ f t, d t
0 3
x t3
2
= ∫ t +
2 dt
0 3
x3 x7
= +
3 63
x t3 t 7
2
φ3 ( x ) = ∫ t + + d t
2
0 3 63
Differential Equations (165)
x
t 6 t14 2 t10
= ∫ t 2 + + 2 + dt
0 9 63 3 × 63
x3 x 7 x15 2 x11
= + + +
3 7.9 15 × 63 × 63 11× 3 × 63
x3 x 7 2 x11 x15
= + + +
3 7.9 11.3.68 15.63.63
(b) φ0 ( x ) = 1 f ( x, y ) = 1 + x y
x
φ1 ( x ) = y0 + ∫ (1 + t φ0 ) d t
x0
x
= 1 + ∫ [1 + t ] d t
0
x2
=1+ x +
2
x t 2
φ2 ( x ) = 1 + ∫ 1 + t 1 + t + d t
0 2
x
t3
= 1 + ∫ 1 + t + t 2 + d t
0 2
x 2 x3 x 4
=1+ x + + +
2 3 8
x t 2 t 3 t 4
φ3 ( x ) = 1 + ∫ 1 + t 1 + t + + + d t
0 2 3 8
x
t3 t 4 t5
= 1 + ∫ 1 + t + t 2 + + + d t
0 2 3 8
x 2 x3 x 4 x5 x6
=1+ x + + + + +
2 3 8 15 48
(c) φ0 ( x) = y0 = 1 f ( x, y) = y 2
x
φ1 ( x) = 1 + ∫ [(1) 2 ] d t
0
= 1+ x
x
φ2 ( x) = 1 + ∫ [1 + t ]2 d t
0
x3 x5 x7 x 2 x3 x 4 x 4 x5 x6
= 1+ x + + + +2 + + + + +
3 5 7 ⋅9 2 3 12 4 15 18
2 4 1 5 x6 x7
φ3 ( x ) = 1 + x + x + x + x + x + +
2 3
3 3 9 63
(b) On what interval containing x = 0 will all the sucdessive approximations exist and be
such that their graphs are in R.
Answers :
(a)
sup
M= R f ( x, y )
sup
= | x|≤1, | y|≤1 ( x
2
+ y2 )
=2
(b) By theorem 4.1.2
I = x / | x – x0 | ≤ α where α = min a,
b
M
Here, x0 = 0, y0 = 0, a = b = 1 and M = 2
1
\ I = x / | x | ≤ α where α = min 1,
2
\ {
I = x/| x| ≤
1
2 }
.
Answers :
(a) f ( x, y ) = 4 x 2 + y 2
f ( x, y1 ) – f ( x, y 2 ) = 4 x 2 + y12 – 4 x 2 + y 22
= y12 – y22
= y1 + y2 y1 – y2
But | y |≤1 ∴ | y1 | ≤ 1 and | y2 | ≤ 1
f ( x, y1 ) – ( f ( x, y2 ) ≤ (| y1 | + | y2 |) | y1 – y2 |
≤ 2 | y1 – y 2 |
Therefore Lipschitz constant K = 2.
( )
= x 2 cos2 y1 – cos 2 y2 + sin 2 x ( y1 – y2 )
≤ cos 2 y1 – cos 2 y2 + y1 – y2
By mean value theorem f (b) – f (a) = f ′(c ) (b – a )
cos2 y1 – cos 2 y2 = –2cos y sin y ( y1 – y 2 )
Therefore
f ( x, y1 ) – f ( x, y2 ) ≤ 2 cos y sin y y1 – y2 + y1 – y2
≤ 3 y1 – y2
Therefore k = 3 is a Lipschitz constant.
5. (a) Show that the function f given by
f ( x, y ) = x 2 | y |
Satisfies Lipschitz condition on R = {( x, y ) | x | ≤ 1, | y | ≤ 1}
∂f
(b) Show that does not exist at (x, 0) if x ¹ 0.
∂y
Differential Equations (168)
Answer :
(a) f ( x, y1 ) – f ( x, y2 ) = x 2 | y1 | – x 2 | y 2 |
≤ | x 2 | | y1 – y2 |
≤ 1| y1 – y2 |
Thus, function satisfies Lipschitz condition with Lipschitz constant k = 1.
∂f
(b) Since | y | is not differentiable at y = 0, do not exist at (x, 0) unless x = 0 if x = 0 then
∂y
the function itself is zero.
EXERCISE
{ 1
on S = ( x, y ) / | x | ≤ , | y | ≤ 1 .
2 }
Differential Equations (169)
Unit 2 : Convergence of successive approximations
In the last unit we have found the successive approximate solutions to a differential equation
(4.1.1). In this unit let us prove that this sequence of successive approximate solutions actually
converges to the exact solution of differential equation (4.1.1).
y ′ = f ( x, y ), y ( x0 ) = y0 on I
By theorem 4.1.2 the functions fp all exist, each is continuous on I and ( x,φ p ( x)) ∈ R for
x in I.
Moreover φ1 ( x ) – φ0 ( x ) ≤ M | x – x0 | for x in I
x x
φ2 ( x ) – φ1 ( x ) = y0 + ∫ f (t , φ1 (t )) d t – y0 + ∫ f (t , φ0 (t )) d t
x0 x0
x
= ∫ [ f (t , φ1 (t )) – f (t, φ 0 (t )) ] dt
x0
But φ1 ( x) – φ0 ( x) ≤ M x – x0 for x in I
x
Therefore φ2 ( x) – φ1 ( x) ≤ K ∫ M t – x0 d t
x0
( x – x0 ) 2
and φ2 ( x) – φ1 ( x) ≤ K M .
2
By mathematical induction we will prove that
| x – x0 | p
φ p ( x ) – φ p –1 ( x ) ≤ M K p –1 for every x in I............ (4.2.1)
p!
We have seen that this inequality is true for p = 1 and p = 2. Let us assume the result for
p = m and we will prove it for p = m + 1.
Without loss of generality assume that x ³ x0.
By definition of fm+1 and fm we get
x x
φm +1 ( x ) – φm ( x ) = y0 + ∫ f (t , φm (t )) d t – y0 + ∫ f (t , φm –1 ) d t
x0 x0
x
= ∫ [ f (t , φ m (t )) – f (t , φ m –1 (t )) ] d t
x0
x
Thus, φm +1 ( x) – φm ( x) ≤ ∫ f (t, φm (t )) – f (t, φm –1 (t )) d t
x0
Since f satisfies Lipschitz condition we get
x
φm +1 ( x) – φm ( x) ≤ K ∫ φm (t ) – φm –1 (t ) d t
x0
m
t – x0
But φm (t ) – φm –1 (t ) ≤ M K m –1
.
m!
MK m x m
Therefore φm +1 ( x ) – φm ( x ) ≤ ∫ t – x0 dt
m! x
0
MK m
( x – x0 )m +1 M ( K | x – x0 |)m +1
≤ =
m! m +1 K ( m + 1)!
= f (x)
Thus the sequence {fk} of successive approximations is a convergent sequence.
x2
= ∫ f (t , φk (t ) d t
x1
Therefore | φ ( x) – φ ( x0 ) | ≤ M | x – x0 | ≤ b
k
and φk ( x ) = φ0 ( x ) + ∑ φ p ( x ) – φ p –1 ( x)
p =1
Therefore
∞
φ ( x) – φk ( x) = ∑ φ p ( x) – φ p –1 ( x)
p = k +1
∞
φ ( x ) – φ k ( x ) ≤ ∑ | φ p ( x ) – φ p –1 ( x ) |
p = k +1
p p
M K | x – x0 |
But φ p ( x) – φ p –1 ( x) ≤ and | x – x0 | < α
K p!
M K pα p
Therefore φ p ( x) – φ p –1 ( x) ≤
K p!
∞ M K pα p
Thus, φ ( x) – φk ( x) ≤ ∑
p = k +1 K p!
M ( K α )k +1 ∞ ( K α ) p
≤ ∑
k (k + 1)! p = 0 p!
k +1
M (K α )
≤ ek α
K (k + 1)!
for every k we have
k +1
M (K α )
φ ( x) – φ k ( x) ≤ eK α .
K ( k + 1)!
(d) The limit f is a solution
We must show that
x
φ ( x) = y0 + ∫ f (t , φ (t )) d t , for all x in I .
x0
x x
∫ f (t , φk (t )) d t – ∫ f (t , φ (t )) d t
x0 x0
x
≤ ∫ f (t , φ k (t )) – f (t , φ (t )) d t
x0
x
≤K ∫ φ k (t ) – φ (t ) d t
x0
But by (c)
k +1
M (K α )
φ k (t ) – φ (t ) ≤ eK α
K (k + 1)!
x x
( K α ) k +1 K α
Therefore ∫ f (t , φk (t )) d t – ∫ f (t , φ (t )) d t ≤ M e | x – x0 |
x0 x0 (k + 1)!
( K α ) k +1
Since → 0 as k → ∞,
(k + 1)!
x x
∫ f (t , φk (t )) d t → ∫ f (t , φ (t )) d t that is
x0 x0
x x
∫ f (t , φ k (t )) d t = ∫ f (t , φ (t )) d t
Lt
k →∞
x0 x0
x
And φ ( x) = y0 + ∫ f (t, φ (t )) d t
x0
Differential Equations