Ode PDF
Ode PDF
INSTITUTE OF MATHEMATICS
189/35 BEHIND RAILWAY STATION, VAISH COLLEGE ROAD, ROHTAK PIN-124001 (HARYANA)
E-mail us on - [email protected], [email protected] Mob. 8607383607, 9802177766
Page 1
Chapter - 4
Method of Variation of Parameters :
Here we shall explain the method of finding the complete solution of a linear differential equation of
second order whose complementary function is known.
d2y dy
Let 2
P Qy R ......(1)
dx dx
be a linear differential equation of the second order where P, Q and R are functions of x.
Working Rule :
d2y dy d2y
Step I : Put the equation in the standard form P Qy R in which the co-efficient of
dx 2 dx dx 2
must be unity.
d2y dy
Step II : Taking R = 0 in standard form , we get 2
P Qy 0 . Then find two independent
dx dx
solutions of this equation and call them y1 and y2. So, we have
y1 y2 y2 R yR
Step IV : Find W , u(x) = dx c1 and v x 1 dx c2
y1 y2 W W
Step V : Putting the values of u(x) and v(x) in y = u(x)y1 + v(x) y2 , we get the required solution.
d2y
Example 1 : Apply the method of variation of parameters to solve 4 y tan 2 x
dx 2
d2y
Solution : The given differential equation is 4 y tan 2 x ......(1)
dx 2
d2y dy
Comparing (1) with 2
P Qy R, we get P 0, Q 4, R tan 2 x
dx dx
The symbolic form of the equation (1) is
d
( D 2 4) y tan 2 x where D =
dx
2
A.E. is m2 4 0
or m2 4 m = 2i
cos 2 x sin 2 x
= 2cos 2 2 x 2sin 2 2 x = 2
2sin 2 x 2cos 2 x
y2 R
Now u(x) = dx c1
W
(sin 2 x ) tan 2 x
= dx c1
2
1 sin 2 2 x
= dx c1
2 cos 2 x
1 1 cos 2 2 x
=
2 cos 2 x
dx c1
1
= (sec 2 x cos 2 x )dx c1
2
1 log(sec 2 x tan 2 x) sin 2 x
= c1
2 2 2
1
= log(sec 2 x tan 2 x) sin 2 x c1
4
y1 R
and v(x) = dx c2
W
cos 2 x tan 2 x
= dx c2
2
1
= sin 2 x dx c2
2
cos 2 x
= c2
4
From equation (2), we have
y = u ( x) cos 2 x v( x)sin 2 x
J.R. INSTITUTE OF MATHEMATICS
189/35 BEHIND RAILWAY STATION, VAISH COLLEGE ROAD, ROHTAK PIN-124001 (HARYANA)
E-mail us on - [email protected], [email protected] Mob. 8607383607, 9802177766
Page 3
1 cos 2 x
= log(sec 2 x tan 2 x) sin 2 x cos 2 x c1 cos 2 x sin 2 x c2 sin 2 x
4 4
cos 2 x
or y = log(sec 2 x tan 2 x) sin 2 x sin 2 x c1 cos 2 x c2 sin 2 x
4
So the complete solution of (1) is
cos 2 x
y = c1 cos 2 x c2 sin 2 x [log(sec 2 x tan 2 x)]
4
The Method Of Undetermined Coefficients :
Def. UC functions : A function is called UC function if it is either
(I) A function defined by one of the following
(i) x n , n being a positive integer or zero.
(ii) e ax , a being a non-zero constant.
(iii) sin(bx + c), where b and c are constants and b 0.
(iv) cos(bx +c), where b and c are constants and b 0.
or
(II) A function defined as a finite product of two or more function of these four types. For example
x 2 , e 3 x , sin(3x / 4), x 2e3 x , x3 cos 2 x , xe 4 x sin 3 x are all UC function.
The method of undetermined coefficients (briefly called UC method) is applicable only when the non-
homogeneous function is a finite linear combination of UC functions.
Def. UC set : Consider a UC function f. The set of functions consisting of f itself and all linearly
independent UC functions of which the successive derivatives of f are either constant multiples or
linear combinations is called the UC set of f. For example
(i) consider the function f(x) = x3
then
f ( x ) 3 x 2 , f ( x ) 6 x , f ( x) 6 ,......, f n ( x) 0 , n 3
We see that all the successive derivatives of f are either multiples or linear combinations of the
functions x2, x , 1. Thus the UC set of x3 is the set
S {x 3 , x 2 , x ,1}
(ii) As another example, consider the function f(x) = x2sin x
then f ( x ) 2 x sin x x 2 cos x
Working Rule : We now outline the method of undetermined co-efficient for finding a particular
d2y dy
integral yp of the equation a0 2
a1 a2 y f ( x)
dx dx
where f is a finite linear combination.
i.e., f = 1 1 A2u2 .......... Amu m
Au of UC functions u1 , u2,......., um where
A1 , A2 ......... Am are arbitrary constant. Assuming that the complementary function yc has already been
obtained, we proceed as follows :
Step I. For each of the UC functions u1 , u2 ,........., um , form the corresponding UC set, thus obtaining
Step II. In case one of the UC sets , so formed, say Sj is identical with or completely contained in
another, say Sk ; we omit the set Sj from further consideration, retaining the set Sk.
Step III. Now we consider the UC sets which remain after step (II).
In case one of these UC sets, say ‘Si’ includes one or more members which are solutions of
the corresponding homogeneous equation (i.e. appear in the complementary function), we
multiply each member of Si by the lowest positive integral power of x so that the resulting
revised set contains no member which is a part of complementary function. We then replace
Si by this new revised set. This step is repeated for all UC sets one by one.
J.R. INSTITUTE OF MATHEMATICS
189/35 BEHIND RAILWAY STATION, VAISH COLLEGE ROAD, ROHTAK PIN-124001 (HARYANA)
E-mail us on - [email protected], [email protected] Mob. 8607383607, 9802177766
Page 5
Step IV. Now, we are left with (i) Certain of the original UC sets, which are neither omitted in step
(II) nor revised in step (III) and (ii) Certain revised set resulting form step (III).
Now form a linear combination of all the sets of these two categories, with unknown constant
co-efficients (called undetermined coefficients).
Step V. Assuming that this linear combination is a particular solution of given equation we can find
the unknown contstant by substituting the linear combination in the differential equation.
This completes the method.
The following examples will illustrate the method more clearly.
Example 2 : Find the general solution using the method of undetermined co-efficent,
d2y dy
2
2 3 y 2e x 10sin x
dx dx
d2y dy
Solution : The corresponding homogeneous equation is 2
2 3 y 0 and the complementary
dx dx
function is yc c1e3 x c2e x
The non-homogeneous terms is the linear combination of the two UC functions ex and sinx.
1. Form the UC set for each of these two functions. i.e. S1 {e x }, S2 {sin x, cos x}
2. Note here, that neither of S1 and S2 is identical with nor included in the other, hence both are
retained.
3. By examining the complementary function, we see that no element of S1 and S2 is a part of C.F. ,
so no revision is required.
4. Thus the original sets S1 and S2 remain as such and we form the linear combination,
Ae x B sin x C cos x of the elements of S1 and S2 with the undetermined coefficients A , B , C
5. We take y p Ae x B sin x C cos x as a particular solution, then
Exercise 4.1
Apply the method of variation of parameters to solve the following differential equations :
d2y 2 d2y
1. 2
y 2. a 2 y sec ax
dx 1 ex dx 2
Using the method of undetermined coefficients find the general solution of the following differential
equation.
d 2 y dy
3. y " 3 y ' 2 y 14sin 2 x 18cos 2 x 4. 6 y 10 e 2 x 18 e3 x 6 x 11
dx 2 dx
Answers
ex 1
1. y c1e x c2 e x e x log x 1 e x log(e x 1)
e
cos ax log cos ax x
2. y c1 cos ax c2 sin ax sin ax
a2 a
3. y c1e x c2e 2 x 2sin 2 x 3cos 2 x
4. y c1e2 x c2e 3 x 2 x e 2 x 3 e3 x x 2
J.R. INSTITUTE OF MATHEMATICS
189/35 BEHIND RAILWAY STATION, VAISH COLLEGE ROAD, ROHTAK PIN-124001 (HARYANA)
E-mail us on - [email protected], [email protected] Mob. 8607383607, 9802177766
Page 1
Assignment-4 1. x5 2. x
--------------------- S C Q --------------------- 3. tan x 4. sin x
1. Solving by variation of parameter 6. A UC set of e2x sin 3x is
y " 2 y ' y e x log x , then the value of
1. e2 x
wronskian is
1. e2x 2. 2
2. e2 x sin 3x
3. e2x 4. –2 3. e2 x cos3x
d2y
2. For
dx 2
4 y tan 2 x solving by variation 4. e2 x sin 3x, e2 x cos3x
of parameters. The value of wronskian W is 7. Using the method of undetermined
1. 1 2. 2 3. 3 4. 4 coefficients, find the general solution of
3. Solving by variation of parameter for the y " 3 y ' 2 y 14sin 2 x 18cos 2 x
equation y " y sec x , the value of 1. y c1e x c2e x 2sin 2 x 3cos 2 x
wronskian is
2. y c1e x c2e 2 x 2sin 2 x 3cos 2 x
1. 1 2. 2 3. 3 4. 4
4. Using the method of variation of 3. y c1e2 x c2e 2 x 2sin 2 x 3cos 2 x
parameters for the particular solution to the 4. y c1e x c2e 2 x sin 2 x 3cos 2 x
3
differential equation y " 4 y ,
sin 2 x
0 x Answers
2
3 3 ----------------------- S C Q --------------------
1. sin 2 x log cos 2 x cos 2 x
4 4 1. 1 2. 2 3. 1
3 3 4. 4 5. 3 6. 4
2. sin 2 x log cos 2 x cos 2 x
2 4 7. 2
3 3
3. sin 2 x log sin 2 x x cos 2 x
2 2
3 3
4. sin 2 x log sin 2 x x cos 2 x
4 2
5. Which of the following are not UC
function.
J.R. INSTITUTE OF MATHEMATICS
189/35 BEHIND RAILWAY STATION, VAISH COLLEGE ROAD, ROHTAK PIN-124001 (HARYANA)
E-mail us on - [email protected], [email protected] Mob. 8607383607, 9802177766
Page 1
Chapter - 5
Second order differential equations
Def. Linearly dependent function : Let f1 , f 2 ,...., f n are real function defined on an interval I, then
f1 , f 2 ,...., f n are said to be L.D. over I if, there exist real no. c1 , c2 ,...., cn (not all zero) such that
c1 f1 x c2 f 2 x .... cn f n x 0, x I
Def. Linearly Independent function : Let f1 , f 2 ,..., f n are real function defined on an interval I, then
have c1 c2 ... cn 0
Result 1 : Let f, g are any function defined on an interval I, then f and g are L.D. over I iff there exist a
constant c such that either f x c g x or g(x) c f x , x I
Result 2 : Suppose f1 , f 2 ,..., f n are defined on an interval I, and J be any interval s.t J I , then we
have
(i) If f1 , f 2 ,..., f n are L.D. on I then f1 , f 2 ,..., f n are L.D. on J.
(iii) If f1 , f 2 ,..., f n are L.I. on I then f1 , f 2 ,..., f n may or may not L.I. on J.
(iv) If f1 , f 2 ,..., f n are L.D. on J then f1 , f 2 ,..., f n may or may not L.D. on I.
Exercise 5.1
1. Which of the following functions are L.D. / L.I. on :
(i) {1, x} (ii) 1, x 2 (iii) x, x 2
(iv) x m
, xn , m n
(ix) {x, cos x} (x) x 2 ,cos x (xi){sin x, cos x} (xii) {sin x, sin 2x}
(xiii) {sin x, 2sin x} (xiv) sin x, sin x (xv) cos x, cos3 x (xvi) cos x, cos 2 x
(xxix) e x ,sinh x (xxx) e x , cosh x (xxxi) sin x,sinh x (xxxii) cos x, cosh x
2. Find the interval on which the following functions are L.D. and L.I.
2
2 x 2 , x 0 sin x, x 0
(i) f x x , g x 2 (ii) f x sin x , g x
3 x , 0 x sin x, 0 x
x
x e , x 0
(iii) f x e , g x x
e , 0 x
3. Find the interval on which the following function are L.D. / L.I.
(i) f x x , g x x (ii) f x x 2 , g x x 2 (iii) f x x3 , g x x 3
(iv) e x ,sin x,cos x (v) sinh x,cosh x,cos x (vi) sin x,cos x,sin x
3
(x) e x , e x , cosh x (xi) 2sin x,3cos x,cos x (xii) 1, x, x
4
Answers
1. (i) L.I. (ii) L.I. (iii) L.I. (iv) L.I. (v) L.I. (vi) L.I. (vii) L.I.
(viii) L.I. (ix) L.I. (x) L.I. (xi) L.I. (xii) L.I. (xiii) L.D. (xiv) L.D.
(xv) L.I. (xvi) L.I. (xvii) L.I. (xviii) L.I. (xix) L.I. (xx) L.I. (xxi) L.I.
(xxii) L.I. (xxiii) L.I. (xxiv) L.I. (xxv) L.I. (xxvi) L.I. (xxvii) L.I. (xxviii) L.I.
(xxix) L.I. (xxx) L.I. (xxxi) L.I. (xxxii) L.I. (xxxiii) L.I. (xxxiv) L.I.
2. (i) L.I. on and L.D. on (,0) and [0, )
(ii) L.I. on and L.D. on (,0] and [0, )
(iii) L.I. on and L.D. on (,0] and [0, )
3. (i) L.I. on and L.D. on (,0] and [0, )
(ii) L.D. on
(iii) L.I. on and L.D. on (,0] and [0, )
(iv) L.I. on and L.D. on (,0] and [0, )
J.R. INSTITUTE OF MATHEMATICS
189/35 BEHIND RAILWAY STATION, VAISH COLLEGE ROAD, ROHTAK PIN-124001 (HARYANA)
E-mail us on - [email protected], [email protected] Mob. 8607383607, 9802177766
Page 3
---------------------------------------------------------------------------------------------------------------------------
Wronskian theory
dny d n 1 y dy
Def. Linear O.D.E. of order n : a0 x n
a1 x n 1
... an1 x an x y F x
dx dx dx
If F x 0 then equation is homogeneous and if F x 0 then equation is nonhomogeneous.
Linear : In any term there should not be two or more y’s multiply together. In other words, in any
dy dny
term there should not be product of two or more terms of the set y , ,........, n .
dx dx
2
dy dy 2 3 d 2 y
e.g, y. , , y , y , y. 2 etc. should not be present .
dx dx dx
Result 3 : FEUT : Consider the linear O.D.E. of order n,
dny d n 1 y dy
a0 x n
a1 x n 1
... an1 x an x y F x ......(1)
dx dx dx
Requirements :
(i) Continuity requirement : The functions a0 x , a1 x ,..., an x and F x are continuous on an
interval I.
(ii) Non zero requirement : a0 x 0, x I
Conclusion : Let c0 , c1 , c2 ,..., cn1 are any real constant and x0 I be any point then exist a unique
solution f of (1) s.t. f x0 c0 , f x0 c1..., f n 1 x0 cn 1 . Moreover this solution f is defined over the
entire interval I.
4
d2y dy
Illustration : If we take x 2 2
2 x x 3 y e x on 1,10 , then both requirement are fulfilled so
dx dx
FEUT is applicable. But if we take the interval 1,1 , then nonzero requirement is not fulfilled and
Def. Wronskian : Let f1 , f 2 ,..., f n be any n function defined on an interval I. Then the determinant
f1 x f 2 x ... fn x
f1 x f 2 x ... f n x
W f1 , f 2 ,..., f n x is called the Wronskian of f1 , f 2 ,..., f n at the
... ... ... ...
f1 n 1 x f 2 n1 ... f n n1 x
point x.
f1 x f2 x
In particular, W f1 , f 2 x f1 x f 2 x f 2 x f1 x
f x
1 f 2 x
Result 4 : In all the following results we shall consider the linear homogenous ODE
d2y dy
a0 x 2
a1 x a2 x y 0 .......(1)
dx dx
where both the requirements are fulfilled ( continuity and non zero ) on an interval I.
Results :
5. Abel’s Liouville’s formula : Suppose f1 , f 2 be any two solutions of (1) and x0 I be any point
x a t
then W x W f1 , f 2 x W f1 , f 2 x0 exp 1 dt , x I
x a0 t
0
Remark : By this formula, the value of wronskian can be calculated at the general point if it is given at
a particular point x0 .
6. Wronskian of any two solutions of (1) is either identically zero on I or never zero on I.
Note : This result is only for solutions of (1).
7. Two solutions f1 and f 2 of (1) are L.D. over I iff W f1 , f 2 x 0, x I
9. Let f1 , f 2 be any two function defined on I such that W f1 , f 2 x 0 , for some x I then, f1 , f 2
are L.I
10. If W f1 , f 2 x 0, x I , then f1 , f 2 are may or may not be L.D.
x 2x
e.g, Consider the functions x, 2 x on . Then W x, 2 x 0 and function are clearly L.D.
1 2
J.R. INSTITUTE OF MATHEMATICS
189/35 BEHIND RAILWAY STATION, VAISH COLLEGE ROAD, ROHTAK PIN-124001 (HARYANA)
E-mail us on - [email protected], [email protected] Mob. 8607383607, 9802177766
Page 5
11. Let f1 , f 2 are any two solution of (1) then W f1 , f 2 x keeps a constant sign on I, if it is non zero.
12. Let f1 , f 2 be any two solution of (1) and they have a common zero in interval I then f1 and f 2 are
L.D.
13. Contrapositive of above result : If two solution of (1) are L.I. then they can not have a common
zero in I.
14. Let f1 , f 2 are any two solution of (1) and x0 I be any point. If f1 and f 2 has a local maxima or
f1 x0 f 2 x0
W f1 , f 2 x0 0
0 0
15. Contrapositive : If the solutions f1 and f 2 of (1) are L.I. then both of them can not have a local
d f1 W f1 , f 2 x
16. Let f1 and f 2 be any two solutions of (1) s.t. f 2 x 0, x I then x
dx f 2 f 22 x
d f1 f 2 x f1' x f1 x f 2' x W f1 , f 2 x
Proof : x 2
=
dx f 2 f2 x f 22 x
f1
17. Let f1 , f 2 be any two solutions of (1) s.t. f 2 x 0, x I , then is a monotonic function.
f2
d f1 W f1 , f 2 x
Proof : We know that x
dx f 2 f 22 x
6
d f1 f1
Case (i): If W f1 , f 2 0 on I then x 0 for all x I is decreasing function on I
dx f 2 f2
d f1 f1
Case (ii): If W f1 , f 2 0 on I then x 0 for all x I is increasing on I
dx f 2 f2
18. The collection of all solutions of a linear homogeneous ODE of order n forms a vector space of
dim n.
d2y dy
Illustration : 2
5 6 y 0 , General solution = c1e2 x c2e3 x , V is a vector space of dim = 2,
dx dx
Basis = e 2 x , e3 x
24. The collection of all solutions of non homogeneous ODE does not form a vector space.
25. Let f1 , f 2 be a set of L.I. solutions of (1) and g1 , g 2 be another set of L.I. solution of (1) then ,
Proof : If f1 , f 2 be one set of L.I solutions of (1) and g1 , g 2 be another set of L.I. solutions of (1).
Then g1 and g 2 are also the linear combination of f1 and f 2 , so result is proved.
26. Let f1 and f 2 are two L.I. solution of (1) and x0 I be any point such that f1 x0 f 2 x0 0 then
a1 x0 a2 x0 0 .
a b
0
c d
or
a b
Let f and g are L.I solution of (1) , then af bg , cf dg are also L.I. solution of (1) iff 0
c d
28. If f x is any non trivial solution of (1) then f x cannot have a multiple zero in I.
J.R. INSTITUTE OF MATHEMATICS
189/35 BEHIND RAILWAY STATION, VAISH COLLEGE ROAD, ROHTAK PIN-124001 (HARYANA)
E-mail us on - [email protected], [email protected] Mob. 8607383607, 9802177766
Page 7
or
2
If f x is a solution of (1) and there is a point x0 I s.t f x x x0 g x s.t. g x0 0 then
Remark : The above result is also valid if the power 2 is replaced by 3, 4, 5, .....
Exercise 5.2
1. Find the Wronskian of the following :
(i) 1, x, x 2 (ii) 1,sin x, cos x (iii) 1, x, e x ,sin x,cos x
(iv) 1, x, e x ,sin x,cos x,sinh x (v) 1, x, e x ,sin x,cos x,cosh x (vi) e x ,sinh x,cosh x
2. Let y1 and y2 be two solutions of the problem y "(t ) ay '(t ) by (t ) 0, t y (0) 0, where
a and b are real constants. Let W be the wronskian of y1 and y2 . Then find W (t ), for all t .
2
y " 2 xy ' sin e 2 x y 0, x [0,1] with y1 (0) 0, y1' (0) 1, y2 (0) 1, y2' (0) 1 . Then find the
wronskian W ( x) of y1 ( x ) and y2 ( x ) at x 1 .
d 2u du
ODE : 2
P( x ) Q( x )u 0, x 1,1 satisfying u1 0, u2 0 and u1 (0) 0, u2 (0) 0 .
dx dx
Let W denote the wronskian of u1 and u2 . Then find W. Also is u1 and u2 are linearly dependent
or not.
5. Let P be continuous function on and W the wronskian of two linearly independent solutions y1
d2y dy
and y2 of the ODE : 2
(1 x 2 ) P( x ) y 0, x . Let W (1) a, W (2) b, W (3) c . Then
dx dx
find a, b, c .
Answers
1. (i) 2 (ii) –1 (iii) 2e x (iv) –4 (v) 4 (vi) 0
2. W (t ) 0, t 3. e 4. W 0, L.D. 5. a b c or a b c
8
Strum Theory
d2 d
2 0
a x y a1 x y a2 x y 0. ......(2)
dx dx
Remark : After taking the indicated derivatives in equation (2), we obtain
d2y dy
a0 x 2a0 x a1 x a0 x a1 x a2 x y 0 ......(3)
dx
2 dx
where the primes denote differentiation with respect to x.
d2y dy
Def. The second order homogeneous linear differential equation a0 x 2
a1 x a2 x y 0 is
dx dx
called self adjoint if it is identical with its adjoint equation.
Theorem 1 : A necessary and sufficient condition for equation (1) to be self adjoint is that
d
a0 x a1 x on a x b.
dx
d2y dy
Corollary : If the equation a0 x 2
a1 x a2 x y 0 is self adjoint, then it can be written in
dx dx
d dy
the form a0 x a2 x y 0 .
dx dx
Proof : Since the given equation is self adjoint, so by Theorem 1, we have a0 x a1 x . Thus the
d2y dy
given equation may be written as a0 x 2
a0 x a2 x y 0 or
dx dx
d dy
a0 x a2 x y 0.
dx dx
Theorem 2 : Let the coefficients a0 , a1 and a2 in the differential equation
d2y dy
a0 x 2
a1 x a2 x y 0 ......(1)
dx dx
are continuous on a x b and a0 x 0 on a x b .
J.R. INSTITUTE OF MATHEMATICS
189/35 BEHIND RAILWAY STATION, VAISH COLLEGE ROAD, ROHTAK PIN-124001 (HARYANA)
E-mail us on - [email protected], [email protected] Mob. 8607383607, 9802177766
Page 9
d dy
Then equation (1) can be transformed into the equivalent self adjoint equation P x Q( x) y 0
dx dx
1 a x
on a x b , by multiplication throughout by the factor exp 1 dx
a0 x a0 x
a x a x a x
Here, P x exp 1 dx , Q(x)= 2 exp 1 dx .
a0 x a0 x a0 x
d2y dy
Example : Consider the equation x 2 2
2x 2 y 0 .
dx dx
1 a x
adjoint. Let us form the factor exp 1 dx for this equation.
a0 x a0 x
1 a x 1 2 x 1
We have exp 1 dx = 2 exp 2 dx 4 .
a0 x a0 x x x x
d2y dy 1
Multiplying equation x 2 2
2t 2 x 0 by 4 on any interval a x b which does not include
dx dx x
1 d 2 y 2 dy 2
x=0, we obtain y 0.
x 2 dx 2 x 3 dx x 4
d 1 2 d 1 dy 2
Since 2 3 , this equation is self adjoint and may be written in the form y 0.
dx x x dx x 2 dx x 4
d dy
Theorem 3 : Let f be a solution of P x +Q x y 0 having first derivative f on a x b .
dx dx
d dy
Theorem 4 : Abel’s Formula :Let f and g be any two solutions of P x +Q x y 0 on the
dx dx
d dy
Theorem 5 : Let f and g be two solutions of P x +Q x y 0 s.t. f and g have a common
dx dx
zero on a x b . Then f and g are linearly dependent on a x b .
10
d dy
Theorem 6 : Let f and g be nontrivial linearly dependent solutions of equation P x Q( x) y 0
dx dx
Theorem 7 : Sturm Separation Theorem : Let f and g be real linearly independent solutions of
d dy
P x +Q x y 0 on the interval a x b .
dx dx
Between any two consecutive zeros of f there is exactly one zero of g.
Theorem 8 : Sturm’s Fundamental Comparison Theorem : On the interval a x b , Let 1 be a
d dy d dy
real solution of P x +Q 1 x y 0 . Let 2 be a real solution of P x +Q 2 x y 0 .
dx dx dx dx
Let P have a continuous derivative and be s.t. P x 0 and let Q 1 and Q 2 be continuous and s.t.
Q 2 x Q 1 x . If x1 and x2 are successive zeros of 1 on a, b , then 2 has atleast one zero at some
Exercise 5.3
1. Find the adjoint equation to each of the following equations :
d2y dy d2y dy
(i) x 2 2
3x 3 y 0 (ii) 2 x 1 2
x3 y 0
dx dx dx dx
d2y dy d2y dy
(iii) x 2 2
2 x 3 7 x 8 x 2 8 y 0 (iv) x 3 2
x 3 2 x 2 x x 2 x 1 y 0
dx dx dx dx
d2y dy
2. Show that the adjoint equation of the adjoint equation of the equation a0 x 2
a1 x a2 x y 0
dx dx
is the original equation itself.
3. Show that each of the following equations is self adjoint and write each in the form
d dy
a0 x a2 x y 0 .
dx dx
d2y dy d2y dy 2
x 1 d y 1 dy 1
(i) x 3 3x2 y0 (ii) sin x cos x 2 y 0 (iii) 2
2 y 0
dx 2
dx dx 2 dx x dx x dx x 3
4. Transform each of the following equations into an equivalent self adjoint equation :
d2y dy d2y dy
(i) x 2
dx 2
x y 0
dx
(ii) x 4 x2 dx 2
2 x3 3 x 0
dx
d2y dy d2y
(iii) 2
tan x y 0 (iv) f x g x y 0
dx dx dx 2
J.R. INSTITUTE OF MATHEMATICS
189/35 BEHIND RAILWAY STATION, VAISH COLLEGE ROAD, ROHTAK PIN-124001 (HARYANA)
E-mail us on - [email protected], [email protected] Mob. 8607383607, 9802177766
Page 11
2
5. Let x :[0,3 ] be a nonzero solution of the ODE x "(t ) et x (t ) 0, for t [0,3 ] . Then find
the cardinality of the set {t [0,3 ] : x (t ) 0} .
Answers
1. (i) x 2 y " xy ' 2 y 0 (ii) (2 x 1) y " (4 x 3 ) y ' (1 3x 2 ) y 0
d 3 dy d dy
3. (i) x y0 (ii) sin x 2 y 0
dx dx
dx dx
d ( x 1) dy 1
(iii) y0
dx x dx x3
d dy 1 d 2 dy 3
4. (i) x y0 (ii) ( x 1) 2 y 0
dx dx x
dx dx x
d dy d g ( x ) dy
(iii) cos x cos x y 0 (iv) exp dx 0
dx dx dx f ( x ) dx
5. 3
---------------------------------------------------------------------------------------------------------------------------
Sturm Liouville Boundary Value Problems
These problems represents a class of linear boundary value problems. The importance of these
problems lies in the fact that they generates set of orthogonal functions. The sets of orthogonal
functions are useful in the expansion of a certain class of functions.
S L equation : A classical “strum-lioulle equation”, is a real second-order linear differential equation
of the form
d dy
P ( x ) q( x) y r ( x ) y ……(1)
dx dx
In the simplest of cases all coefficients are continuous on the finite closed interval a, b , and P( x )
required to satisfy boundary conditions. The function r ( x) , is called the “weight” or “density”
function. The number of famous differential equations could be represented in the SL form :
Bessel’s equation : x 2 y " xy ' ( x 2 v 2 ) y 0 can be written in strum-liouville form as
v2
( xy ') ' x y 0
x
The Legendre equation : (1 x 2 ) y " 2 xy ' v(v 1) y 0 can easily be put into SL form, since
The general way to convert the 2nd order linear ODE to the SL form is to use an integrating factor
u ( x) such that the P( x) y " Q ( x) y ' R( x ) y 0 multiplied by u ( x) would has SL form, one can
1 Q( x)
easily show that u ( x) exp dx does the job.
P( x) P( x)
SL boundary value problem (SL-BVP) : We introduce the SL-operator as
d dy
L[ y ] P( x ) q( x) y and consider the SL equation L[ y ] r ( x) y 0 ……(2)
dx dx
where P( x ) 0, r ( x) 0 and P, q and r are continuous functions on the interval [a, b]; along with
the BC
1 y (a ) 2 P (a) y '(a) 0
……(3)
1 y (b) 2 P(b) y '(b) 0
Page 13
where p , q , r are real functions s.t. p(t) has a continuous derivative , q(t) and r(t) are continuous ,
and p(t) > 0 and r(t) > 0 for all t on a real interval a t b , and is a parameter independent
of t , and the conditions
A1u(a) + A2u'(a) = 0
.......(2)
B1u(b) + B2u'(b) = 0
where A1 , A2 , B1 , B2 are real constants s.t. A1 , A2 are not both zero and B1 and B2 are not
both zero.
Let m and n be any two characteristic values of the problem. Let m be a characteristic function
corresponding to m and let n be a characteristic function corresponding to n. Then
characteristic functions m and n are orthogonal w.r.t. the weight function r(t) on the interval
b
a t b i.e. n (t ) m (t ) r(t) = 0
a
r ( x) u( x) v( x) dx 0
a
4. If A : differential equation has non-trivial solution i.e. A is the set of eigen values. Then A
is always infinite set. Also elements of A can be arranged into strictly increasing sequence
1 2 3 .... n .... such that n as n
A is always bounded below and unbounded above set. Again, every element of A is an isolated
point.
A is no where dense set. Since set of isolated points is countable. So A is always countably
infinite set.
14
5. There does not exist any 0 such that the differential equation has non-trivial solution for
every 0 .
6. There does not exist any 0 such that differential equation has only trivial solution for all
0 .
7. There exists a 0 such that the differential equation has only trivial solution for all 0 .
8. The number of negative eigen values are always finite if exists.
Properties of regular BVP :
1. The set A of eigen values is of the form 1 2 3 .... n .... and n as n
2. The eigen values of the regular BVP are simple. Thus an eigen function that corresponds to an
eigen value is unique upto a constant multiple dim E 1
zeros in (a, b) corresponding of first eigen value, eigen function does not have any zeros and
corresponding to 2 every eigen function has one zero and corresponding to 3 every eigen
( a, b) .
6. If 0 is an eigen value and corresponding eigen function does not change its sign. This implies
there does not exists any negative eigen value.
7. Similarly if 0 is an eigen value and corresponding eigen function change its sign, then there
must exists negative eigen values.
Properties of periodic BVP : Let A : differential equation has non-trivial solution
1 2 3 .... n ....
1. dim E0 1
2. dim Ek 2, for all k 1
3. If 0 is an eigen value and dim E0 1 , then there does not exists any negative eigen value.
4. If 0 is an eigen value and dim E0 2 then there must exists negative eigen value.
Page 15
d du
and hence (1) is of the form
dt p(t) dt q(t) + r(t)u = 0
where p(t) = 1 , q(t) = 0 and r(t) = 1. The supplementary conditions are of the special form
u(a) = 0 , u(b) = 0.
d du
Example 2 : Show that the boundary value problem t 2t 2 t 3 u = 0 ......(1)
dt dt
with conditions 3 u(1) + 4 u' (1) = 0 ..... (2)
5u(2) – 3 u' (2) = 0
is a Sturm – Liouville problem.
d du
Solution : The differential equation is of the form
dt p(t) dt q(t) + r(t)u = 0
where p(t) = t , q(t) = 2t2 and r(t) = t3. The conditions (2) are of the form
1u(a) + 2u'(a) = 0
1u(b) + 2u'(b) = 0
where a = 1 , b = 2 , 1 = 3 , 2 = 4 , 1 = 5 , 2 = – 3.
Example 3 : Find non-trivial solutions of Sturm Liouville boundary value problem
d 2u
u = 0 ......(1)
dt 2
u(0) = 0 , u( ) = 0 ......(2)
Solution : We shall consider separately the three cases = 0 , < 0 and > 0. In each case , we
shall first find the general solution of the differential equation (1) and then attempt to determine two
arbitrary constants in this general solution so that the supplementary conditions (2) will also be
satisfied.
d 2u
Case(I) : 0 In this case (1) reduces to = 0 and so the general solution is
dt 2
u(t) = c1 + c2t .......(3)
We now apply conditions (2) to solution (3). Condition u(0) = 0 implies
0 = c1 + c2.0 c1 = 0.
and condition u( ) = 0 implies 0 = c1 + c2 c2 = 0 [Since c1 = 0].
16
Thus in order that solution (3) to satisfy conditions (2) , we must have c1 = c2 = 0.
But then the solution (3) becomes u(t) = 0 for all t. Thus , in case when the parameter = 0 , the
only solution of the given problem is the trivial solution.
d 2u
Case(II) : 0 Differential equation (1) is u = 0
dt 2
Its auxiliary equation is m2+ = 0 and m = . Since is negative , so these roots are real
Since = , we must have = 0. But < 0 in this case. Thus there are no non-trivial
solution of the given problem in this case < 0.
Case(III) : 0 In this case A.E. is m2+ = 0. So its roots are m =
These roots are conjugate complex numbers since > 0. Roots can be written as i . Thus in
this case the general solution is of the form
u(t) = c1sin t + c2cos t .......(7)
We now apply conditions (2) to this general solution. Condition u(0) = 0 implies
c1sin 0 + c2cos 0 = 0 c2 = 0
Condition u( ) = 0 implies
c1sin + c2cos = 0
c1sin = 0 .......(8)
[Since c2 = 0].
We must therefore satisfy (8).
J.R. INSTITUTE OF MATHEMATICS
189/35 BEHIND RAILWAY STATION, VAISH COLLEGE ROAD, ROHTAK PIN-124001 (HARYANA)
E-mail us on - [email protected], [email protected] Mob. 8607383607, 9802177766
Page 17
The values of parameter in (1) for which these exist non-trivial solutions of the Sturm Liouville
problem are called the characteristic values (or eigen values of the problem). The corresponding non-
trivial solutions are called the characteristic functions or the eigen functions of the problem.
Example 4 : Find the characteristic functions of the strum-liouville problem
y " y 0, y (0) y ( ) 0, y '(0) y '( ) 0 , A
Solution : y " y 0, y (0) y ( ) 0, y '(0) y '( ) 0
These boundary conditions are called periodic boundary condition. We consider three cases
corresponding to the value of :
Case I : u 2 0
The general solution of the ODE is given as yy Ae ux Be ux
By substituting BC we obtain the following system :
18
B 1 ).
Note that this eigen value is simple. The eigen value is called simple, if its eigenspace is of dimension
one ; otherwise the eigenvalue is called multiple.
Case III : u 2 0 The general solution of the ODE is given as y A cos(ux) B sin(ux )
By substituting BC we obtain the following system :
A(1 cos(u ) B sin(u ) 0
A sin(u ) B sin(1 cos(u )) 0
This problem has a non-trivial solution only when the determinant of the matrix of coefficients
D (u ) 2 2cos 0 . This corresponds to u 2n, n 1, 2,.... and hence n 4n 2 .
n cos
n x , n sin n x
Note that the eigen values n are positive and there are two linearly independent eigen functions
Page 19
d2y
y = 0 .......(4)
dt 2
We consider the separately two cases.
d2y
Case (i) : 0 Differential equation (4) reduces to = 0
dt 2
The general solution of this differential equation is y = c0t + c
or y = c0 logx + c .......(5)
Now we apply the conditions (2) to solution (5). Condition y' (1) = 0 implies
c0 c0
0 = c0 = 0 [Since y' (x) = ]
1 x
Condition y' e 2 = 0 implies , c0 = 0
Both of these conditions does not impose any restriction upon c. Thus for = 0 , we obtain the
solutions y = c where c is an arbitrary constant.
These are non-trivial solutions for all c 0. Thus = 0 is a characteristic value and the
corresponding characteristic functions are given by y = c , where c is an arbitrary non-zero
constant.
Case (ii) : 0 In this case A.E. is m2+ = 0 m = = i .
Condition y' (1) = 0 implies c1 cos ( log 1)– c2 sin ( log 1) = 0 or c1 = 0 [Since log 1= 0]
Thus we must have c1 = 0 .......(8)
Condition y' e 2 = 0 implies
Since c1 = 0, the choice c2 = 0 would give us the trivial solution. So we must have sin 2 = 0 and
hence
2 = n , n = 1 , 2 ,……………..
n2
Thus for the non-trivial solution , we must have = , n = 1 , 2 , 3 ,………… ….....(9)
4
Corresponding to these values of , we obtain the non-trivial solutions
n logx
y = cn cos (n = 1 , 2 , 3,……………..) .......(10)
2
where cn ’s are arbitrary non-zero constants.
A R : y1 ( ) 0 ,
B R : y2 ( ) 0
Then
(i) A and B must be countable.
(ii) A B , i.e., zeros of two L.I. solution never common zero or if two solution have common
zero they must be L.D.
(iii) Between any two successive zeros of y1 ( x), exactly one zero of y2 ( x ) and conversely.
Page 21
y "( x ) q( x ) y 0
Case II. q( x ) is continuous q( x ) 0 x R assuming that q( x ) dx
1
So 0 m q ( x ) M for all x,
m and M are infimum and supremum respectively.
(iv) If y ( x) be a non trivial solution and if x0 , x1 are successive zeros of y ( x) , let x0 x1 then
x1 x0
M m
(i) Distance between two successive zeros of any non-trivial solution are lying between ,
M m
K 2 2
(ii) If m 2
then y ( x) has atleast K zeros in a, b .
b a
22
Exercise 5.4
Find the characteristic values and characteristic functions of each of the following Sturm Liouville
problems.
d2y d2y
1. y 0 ; y(0) = 0, y 0. 2. y 0 ; y(0) = 0, y 0.
dx 2 dx 2 2
d2y d2y
3. y 0 ; y(0) = 0, y 2 0. 4. y 0 ; y(0) = 0, y l 0.
dx 2 dx 2
d2y d2y
5. y 0 ; y(0) = 0, y 0. 6. y 0 ; y(0) = 0, y 0.
dx 2 dx 2 2
d2y d2y
7. y 0 ; y(0) = 0, y 2 0. 8. y 0 ; y(0) = 0, y l 0.
dx 2 dx 2
d2y d2y
9. y 0 ; y (0) = 0, y 0. 10. y 0 ; y (0) = 0, y 0.
dx 2 dx 2
2
d2y d2y
11. y 0 ; y (0) = 0, y 2 0. 12. y 0 ; y (0) = 0, y l 0.
dx 2 dx 2
d2y d2y
13. y 0 ; y (0) = 0, y 0. 14. y 0 ; y (0) = 0, y 0.
dx 2 dx 2
2
d2y d2y
15. y 0 ; y (0) = 0, y 2 0. 16. y 0 ; y (0) = 0, y l 0.
dx 2 dx 2
d2y d dy
17. y 0 ; y ( ) = 0, y 0. 18. x y 0, y(1) = 0, y e 0.
dx 2 dx dx x
d dy d2y
19. x y 0 ; y(1) = 0, y e 0. 20. y 0, y(0) = 0, y y 0 .
dx dx x dx 2
d2y d2y
21. y 0, y 0 y 0 0, y y 0 . 22. y 0, y 0 y 0 0, y 1 y 1 0.
dx 2 dx 2
d2y d dy
23. y 0, y 0 y 2c , y 0 y 2c . 24. x y 0, y 1 y e 2 0
dx 2 dx dx x
d 3 dy
25. x xy 0 ; y(1) = 0, y(e) = 0. 26. y 2 y 1 y 0 ; y(0) = 0, y l 0.
dx dx
d 2 dy
27.
dx
x 1 2
dx x 1
y 0, y(0) = 0, y(1) = 0. [ Hint : Let x = tan t ]
d 1 dy
28. 3 x 2 1 y 0, y(0) = 0, y 0 . [ Hint : Let t x3 x ]
dx 3 x 2 1 dx
J.R. INSTITUTE OF MATHEMATICS
189/35 BEHIND RAILWAY STATION, VAISH COLLEGE ROAD, ROHTAK PIN-124001 (HARYANA)
E-mail us on - [email protected], [email protected] Mob. 8607383607, 9802177766
Page 23
Answers
n2 nx
3. , n 1,2,3,..... y ( x) c sin , n 1, 2,3,.....
4 2
n 2 2 n x
4. 2
, n 1, 2,3,..... y ( x) c sin , n 1, 2,3,.....
l l
2
1 1
5. n , n 0,1, 2,3,..... y ( x) c sin n x, n 0,1, 2,3,.....
2 2
2
6. 2n 1 , n 0,1, 2,..... y ( x) c sin(2n 1) x, n 0,1, 2,3,.....
2
2n 1 2n 1
7. , n 0,1, 2,3,..... y ( x) c sin x, n 0,1, 2,3,.....
4 4
2
8. (2n 1)2 , n 0,1, 2,3,..... y ( x) c sin (2n 1) x, n 0,1, 2,3,.....
4l 2 2l
2
1 1
9. n , n 0,1, 2,3,..... y ( x) c cos n x, n 0,1, 2,3,.....
2 2
2
2n 1 2n 1
11. , n 0,1, 2,3,..... y ( x) c cos x, n 0,1, 2,3,.....
4 4
2 2
12. 2n 1 , n 0,1, 2,3,..... y ( x) c cos (2n 1) x, n 0,1, 2,3,.....
4l 2 2l
n2 nx
15. , n 0,1, 2,3,..... y ( x) c cos , n 0,1, 2,3,.....
4 2
n 2 2 n x
16. 2
, n 0,1, 2,3,..... y ( x) c cos , n 0,1, 2,3,.....
l l
24
21. 1 y ( x) ce x
22. 1 y ( x) ce x
Page 1
d4y d2y
4. D3 D 2 4 D 4 y 0 the differential equation
dx 4
dx 2
0 is
u1 and u2 is 1
4. W ( x ) 0 for x 0,
2
1. W (t ) 2t 1, 0 t 1
14. Let y1 and y2 be two solutions of the
2. W (t ) sin 2 t , 0 t 1
y "(t ) ay '(t ) by (t ) 0, t
3. W (t ) cos 2 t , 0 t 1 problem
y (0) 0
4. W (t ) 1, 0 t 1
Where a and b are real constants. Let W be
12. Let y1( x ) and y2 ( x ) form a fundamental
the wronskian of y1 and y2 . Then
set of solutions of
1. W (t ) 0, t
J.R. INSTITUTE OF MATHEMATICS
189/35 BEHIND RAILWAY STATION, VAISH COLLEGE ROAD, ROHTAK PIN-124001 (HARYANA)
E-mail us on - [email protected], [email protected] Mob. 8607383607, 9802177766
Page 3
2. W (t ) c, t for some positive 18. The set of all eigen values of the strum
constant c Liouville problem y " y 0, y '(0) 0,
3. W is a non-constant positive function
y ' 0 is given by
4. There exists t1, t2 such that 2
2 d2y
3. 2, 2 4. 3, y 0, with boundary conditions
3 dx 2
16. The strum Liouville problem dy dy
y (0) y (2 ), (0) (2 ) are
dx dx
y " 2 y 0, y (0) 0 , y '( ) 0 has its
1. n 2. n 2 2
eigen vectors given by y is equal to
1 3. n 4. n 2
1. sin n x
2 20. The set of real numbers for which the
boundary value
2. sin nx
d2y
problem 2 y 0, y 0 0, y 0
1 dx
3. cos n x
2 has nontrivial solutions is
1. ,0
4. cos nx, where n 0,1, 2,...
17. The eigen values of the strum Liouville
2. n | n is a positive integer
system y " y 0, 0 x , 3. n 2
| n is a positive integer
4.
y (0) 0, y '( ) 0 are (CSIR NET Dec 2017)
n2 (2n 1)2 2 21. Consider the ordinary differential equation
1. 2.
4 4 y " P x y ' Q x y 0 where P and Q
2 2 2
(2n 1) n
3. 4. are smooth functions. Let y1 and y2 be any
4 4
two solutions of the ODE. Let W x be the
4
corresponding Wronskian. Then which of to each non-zero eigen value . There
the following is always true ? correspond
1. If y1 and y2 are linearly dependent then 1. only one eigenfunction
2. two eigen function
x1 , x2 such that W x1 0 and
3. two linearly independent eigen functions
W x2 0
4. two orthogonal eigen functions
2. If y1 and y2 are linearly independent
d2y
3. Let q( x ) y 0, 0 x
then W x 0 x dx 2
dy
3. If y1 and y2 are linearly dependent then y (0) 1, (0) 1, where q( x ) is a
dx
W x 0 x positive monotonically increasing
4. If y1 and y2 are linearly independent continuous function. Then
then W x 0 x 1. y ( x) as x
dy
(CSIR NET June 2018) 2. as x
dx
3. y ( x) has finitely many zeros in [0, )
-------------------------- M C Q ---------------------
4. y ( x) has infinitely many zeros in [0, )
1. Let y1( x ) and y2 ( x ) form a fundamental
4. For the boundary value problem
set of solutions to the differential equation
y " y 0 ; y (0) 0, y (1) 0 , there
y " p( x) y ' q( x) y 0 , a x b , where
exists an eigen value for which there
p( x) and q( x ) are continuous in a, b ,
corresponds an eigen function in (0,1) that
and x0 is a point in (a, b) . Then
1. does not change sign
1. both y1( x ) and y2 ( x ) cannot have a 2. change sign
local maximum at x0 . 3. is positive
Page 5
1 1 3. y '(a ) 0
2 2 2
2. u ' ( x)dx u ( x) dx 0 4. y has infinitely many zeros in a, b
0 0
3. y x x t f t dt
ODE :
0
t x xt f t dt
x
d2y dy
dx 2
1 x 2
dx
P( x) y 0, x . Let is a solution of the given BVP
x 1
4. y x
W (1) a, W (2) b and W (3) c , then x t xt f t dt xt f t dt
0 x
1. a 0 and b 0
(CSIR NET Dec 2017)
2. a b c or a b c
16. Consider the differential equation
3.
a b c
d2y dy
2
2 tan x y 0 defined on
a b c dx dx
4. 0 a b and b c 0 , . Which among the following are
2 2
13. The problem y " (1 x) y y , x (0,1) true ?
y (0) y (1) 0 has a non zero solution 1. there is exactly one solution y y x
with y 0 y ' 0 1 and
1. for all 0
2. for all 0,1 y 2 1
3 3
3. for some (2, ) 2. there is exactly one solution y y x
Page 7
y " 0 y 0 Answers
4. if y1 and y2 are any two solutions then ----------------------- S C Q --------------------
ax b y1 cx d y2 for some 1. 1 2. 3 3. 4 4. 1
a, b, c, d
5. 1 6. 4 7. 3 8. 1
(CSIR NET Dec 2017)
17. Consider the Sturm-Liouville problem 9. 3 10. 4 11. 4 12. 1
0,
(CSIR NET June 2018)
J.R. INSTITUTE OF MATHEMATICS
189/35 BEHIND RAILWAY STATION, VAISH COLLEGE ROAD, ROHTAK PIN-124001 (HARYANA)
E-mail us on - [email protected], [email protected] Mob. 8607383607, 9802177766
Page 1
Chapter - 6
System of linear ordinary differential equations
Consider the homogeneous linear system
dx dy
ax by and cx dy ......(1)
dt dt
where a, b, c, d are real constants.
x t a b dx
Suppose X and A , then above system can be expressed as AX or X t AX
y t c d dt
be v1 and v2 . Then the general solution of (1) is given by X t c1v1e1t c2v2e 2t . Here v1e1t
Case (b): If 1 and 2 are complex and distinct, then case (a) is still valid. But a shorter method is given
as : Suppose v1 is an eigen vector corresponding to the eigen value 1 , then one solution of
system (1) is v1e1t . Let w1 and w2 be real and imaginary parts of v1e1t then it can be proved
that w1 and w2 are two L.I. solutions of (1) and so the general solution of (1) is given by
Case (c) : If 1 2 (real) and dimension of eigen space of is 2 then suppose v1 and v2 are two L.I.
eigen vectors w.r.t. . Then v1et and v2e t are two L.I. solutions of (1) and so the general
Case (d) : If 1 2 (real) and dimension of eigen space of is one, then suppose one eigen vector
is v1 . Now, construct a new vector v2 which satisfy the equation A I v2 v1. Then two
L.I. solutions of (1) are given by v1et and tv1 v2 e t and so the general solution is given by
2
Autonomous Systems
Page 3
Def. Center : The isolated critical point (0, 0) of autonomous system is called a center if there exists a
neighbourhood of (0, 0) which contains a countably infinite number of closed path Cn (n = 1, 2 ,…..) ,
each of which contains (0 , 0) in its interior , and which are such that the diameters of the paths
approach 0 as n , but (0 , 0) is not approached by any path either as t or as t .
x
0
The critical point (0 , 0) of adjoining figure is called a center. Such a point is surrounded by an
infinite family of closed paths , members of which are arbitrarily close to (0 , 0) , but it is not
approached by any path either as t or as t .
4
Def. Saddle Point : The isolated critical point (0 , 0) is called a saddle point if there exist a
neighbourhood of (0 , 0) in which the following two conditions hold :
(i) There exist two paths which approach and enter (0, 0) from a pair of opposite directions as t ,
and there exist two paths which approach and enter (0, 0) from a different pair of opposite
directions as t .
(ii) In each of the four domains between any two of the four directions in (i) there are infinitely many
paths which are arbitrarily close to (0 , 0) but which do not approach (0 , 0) either as t or
as t .
The critical point (0 , 0) of adjoining figure is a saddle point which is such that
D R2 y
B
R3
x
O R1
C
A R4
(i) It is approached and entered by two half-line paths (A0 and B0) as t , these two paths
forming the geometric curve AB.
(ii) It is approached and entered by two half – line path (C0 and D0) as t , these two paths
forming the geometric curve CD.
(iii) Between the four half – line paths described in (i) and (ii) there are four domains R1 , R3 , R3 , R4 ,
each containing an infinite family of semi – hyperbolic like paths which do not approach (0 , 0)
as t or as t , but which become asymptotic to one or another of the four half – line
paths as t and as t .
Def. Focal point/Spiral point : The isolated critical point (0, 0) is called a spiral point (or focal point)
if there exists a neighbourhood of (0 , 0) such that every path C in this neighbourhood has the
following properties :
(i) C is defined for all t > t0 (or for all t < t0 ) for some number t0.
(ii) C approaches (0 , 0) as t (or as t ).
(iii) C approaches (0 , 0) in a spiral – like manner , winding around (0 , 0) an infinite number of times
as t (or as t ).
J.R. INSTITUTE OF MATHEMATICS
189/35 BEHIND RAILWAY STATION, VAISH COLLEGE ROAD, ROHTAK PIN-124001 (HARYANA)
E-mail us on - [email protected], [email protected] Mob. 8607383607, 9802177766
Page 5
The critical point (0 , 0) in the above figure is a spiral point (or focal point). This point is approached
in a spiral – like manner by an infinite family of paths as t (or as t ). Here , while the
paths approach (0, 0) , they do not enter it. That is a point R tracing such a path C approaches O (0 , 0)
as t (or as t ), but the line OR does not tend to a definite direction , since the path
constantly winds about O.
Def. Node : The isolated critical point (0 , 0) is called a node if there exist a neighbourhood of (0 , 0)
such that every path C in this neighbourhood has the following properties :
(i) C is defined for all t > t0 (or for all t < t0) for some number t0.
(ii) C approaches (0 , 0) as t (or as t ).
(iii) C enters (0 , 0) as t [or as t ].
denote the distance between the critical point (0 , 0) and the point R : (x(t) , y(t)) on C. The critical
point (0 , 0) is called asymptotically stable if
(i) It is stable and
(ii) There exist a number 0 > 0 such that if
J.R. INSTITUTE OF MATHEMATICS
189/35 BEHIND RAILWAY STATION, VAISH COLLEGE ROAD, ROHTAK PIN-124001 (HARYANA)
E-mail us on - [email protected], [email protected] Mob. 8607383607, 9802177766
Page 7
(a – )A + bB = 0
a A + (b – ) B = 0 .......(3)
This system obviously has the trivial solution A B 0 . But this would give only the trivial solution
x = 0 , y = 0 of the system (1). Thus we seek non – trivial solution of the system (3). A necessary and
sufficient condition that this system have a non-trivial solution is that the determinant
a b
0 ......(4)
a d
This gives the quadratic equation
2 – (a + d) + (ad – bc) = 0 ......(5)
8
This equation is called the characteristic equation associated with the system (1). Its roots 1 and 2
are called the characteristic roots.
If the pair (2) is to be a solution of the system (1) , then in (2) must be one of these roots.
Note : We summarize our results in the following table :
Page 9
unequal and of the same sign (both negative) , we conclude that the critical point (0 , 0) of (1) is a
node. Since the roots are real and negative , the point is asymptotically stable.
Example 2 : Determine the nature of the critical point (0, 0) of the system
dx
2x 4 y
dt
…...(1)
dy
2x 6 y
dt
Exercise 6.1
Find the general solution of each of the following linear systems:
dx dy dx dy dx dy
(1) 5 x 2 y, 4x y (2) x 2 y, 3x 2 y (3) 3x y, 4x 3 y
dt dt dt dt dt dt
dx dy dx dy dx dy
(4) 3 x 4 y, 2x 3 y (5) x 3 y, 3x y (6) x 4 y, x y
dt dt dt dt dt dt
dx dy dx dy dx dy
(7) x 3 y, 3x y (8) 4 x 2 y, 5x 2 y (9) 3x 2 y, 2x 3 y
dt dt dt dt dt dt
dx dy dx dy
(10) 3 x y, 4x y (11) 5 x 4 y, x y
dt dt dt dt
Find the particular solution of each of the following linear systems:
dx dy
(12) 2 x 7 y, 3x 2 y; x 0 9, y 0 1
dt dt
dx dy
(13) 2 x 8 y, x 6 y; x 0 4, y 0 1
dt dt
dx dy
(14) 6 x 4 y, x 2 y; x 0 2, y 0 3
dt dt
10
x t
Find the general solution of each of the following homogeneous linear systems, where X has
y t
been used.
dX 1 2 dX 1 4 dX 3 1 dX 6 4
(15) X (16) X (17) X (18) X
dt 3 2 dt 1 1 dt 4 1 dt 1 2
Answers
1. x1 c1et c2e3t , x2 2c1et c2e3t
8. x1 2c1 cos3t e3t c2e3t (cos 3t 3sin 3t ), x2 2c1 e3t sin 3t c2 e3t (sin 3t 3cos3t )
14. x1 2 8t e4t , x2 (3 4t ) e 4t
Page 11
Theorem (A) : Basic theorems on non-linear system : We now consider the non-linear real
autonomous system
dx
P ( x, y )
dt
……(*)
dy
Q ( x, y )
dt
We assume that the system (*) has an isolated critical point which we shall choose to be origin (0, 0).
We now further assume that the functions P and Q in the right members of (*) are such that P( x, y )
and Q( x, y ) can be written in the form
P( x, y ) ax by P1 ( x, y )
…..(**)
Q ( x, y ) cx dy Q1 ( x, y )
where
a b
I. a, b, c, d are real constants and 0 and
c d
II. P1 and Q1 have continuous first partial derivatives for all ( x, y ) and are such that
P1( x, y ) Q1 ( x, y )
lim lim 0 …..(1)
( x , y )(0,0) x 2 y 2 ( x , y )(0,0) x 2 y 2
obtained from (2) by neglecting the nonlinear terms P1 ( x, y ) and Q1 ( x, y ) . Both systems have an
isolated critical point at (0, 0). Let 1 and 2 be the roots of the characteristic equation
dx
8x y 2
dt
Example 1 : Determine the nature of the critical point (0, 0) of the system ……(1)
dy 2
6 y 6 x
dt
Page 13
dx
8x
dt
dy
6 y
dt
of the form (3). The characteristic equation of this system is 2 2 48 0 and thus the roots are
1 8 and 2 6 . Since roots are real, unequal and of opposite sign, so critical point (0, 0) of the
given nonlinear system is a saddle point. So it is unstable. By definition, the critical points of this
system must simultaneously satisfy the system of algebraic equations 8 x y 2 0 , 6 y 6 x 2 0
From the second equation of this pair, y x 2 . Then, substituting this into the first equation of the pair,
into the origin 0 in the plane. We now transforms the given system into , coordinates.
d
8 8 2 …….(2)
dt
d
24 6 6 2
dt
The system (2) is of the form (2) and hypotheses of theorem (A) and (B) are satisfied in these
coordinates. To determine the type of the critical point 0 of (2), we consider the system
d
8 8
dt
d
24 6
dt
The characteristic equation of this linear system is 2 2 144 0
The roots of this system are 1 143 i which are conjugate complex with real part not zero, so
0 of non linear system (2) is a spiral point and since real part is positive so it is a unstable
spiral point. So the given system (1) has two real critical point namely
1. critical point (0, 0) ; a saddle point ; unstable.
2. critical point (2, 4) ; a spiral point ; unstable.
14
Exercise 6.2
1. Determine the nature of the critical point (0, 0) of the system
dx
x 4 y x2
dt
dy
6 x y 2 xy
dt
and determine whether or not the point is stable.
2. Determine the nature of critical point (0, 0) of the system
dx
sin x 4 y
dt
dy
sin 2 x 5 y
dt
and determine whether or not the point is stable.
Answers
1. saddle point ; unstable 2. node ; asymptotically stable
J.R. INSTITUTE OF MATHEMATICS
189/35 BEHIND RAILWAY STATION, VAISH COLLEGE ROAD, ROHTAK PIN-124001 (HARYANA)
E-mail us on - [email protected], [email protected] Mob. 8607383607, 9802177766
Page 1
Assignment-6 y t dy
4. Let y t 1 satisfy Ay; t 0
--------------------- S C Q --------------------- y2 t dt
x (t ) 0
1. The general solution and y 0 where A is a 2 2
of the system 1
y (t )
x ' x 2 y constant matrix with real entries satisfying
is given by trace A=0 and detA>0. Then
y ' 4x y
y1 t and y2 t both are
c e3t c2e 3t 1. monotonically decreasing functions of
1. 1
2c e3t c e 3t t.
1 2
2. monotonically increasing functions of t.
ce 3t
3. oscillating functions of t.
2. 1 4. constant functions of t.
c e 3t
2 (CSIR NET Dec 2012)
c e3t c2e3t
3. 1 5. Consider the system of ODE in
2c e3t c e 3t
1 2 dy 0
c e3t c2e 3t 2 , AY , Y 0 , t 0 where
4. 1
dt 1
2c e3t c e 3t 1 1 y t
1 2
A and Y t 1 .
(GATE 2001) 0 1 y2 t
Then
2. The general solution of the system of
1. y1 t and y2 t are monotonically
differential equations
dz increasing for t 0
y 0 2. y1 t and y2 t are monotonically
dx
dy increasing for t 1
z0
dx 3. y1 t and y2 t are monotonically
is given by decreasing for t 0
x x x x
1. y e e , z e e 4. y1 t and y2 t are monotonically
2. y cos x sin x, z sin x cos x decreasing for t 1
3. y sin x cos x, z cos x sin x (CSIR NET Dec 2015)
4. y e x e x , z e x e x
(GATE 2005) 2 1 0 x1 (t )
6. Let A 0 2 1 , x(t ) x2 (t ) and
y1' y2 , y2' ay1 by2 . Then any solution of the first order system
Every solution y ( x) 0 as x , if of the ordinary differential equation
1. a 0, b 0 2. a 0, b 0 x(t ) Ax(t )
3. a 0, b 0 4. a 0, b 0 satisfies
x (0) x0
(GATE 2008) 1. lim | x(t ) | 0
t
2
2. lim | x(t ) | 1
t
dy
3. lim | x(t ) | 2 Ay; t 0; y 0 1 . Then
t dt 1
4. lim | x(t ) | 12
t 2
t
(CSIR NET June 2016) 1. y1 t 1 t , y2 t 1 t , y3 t 1.
2
t2
---------------------- M C Q -------------------- 2. y1 t 1 t , y2 t 1 t , y3 t 1.
2
1. Consider the first order system of linear t2
3. y1 t 1, y2 t 1 t , y3 t 1 t .
dX 3 2 2
equations AX ; A ; tA
dt 2 1 4. y1 t e y 0 .
x t (CSIR NET Dec 2013)
X 1 . Then
x2 t 4. The critical point of the system
1. the coefficient matrix A has a repeated dy dy
eigenvalue =1 4 x y, x 2 y is an
dx dt
2. there is only one linearly independent
1. asymptotically stable node
1 2. unstable node
eigenvector X 1
1 3. asymptotically stable spiral
3. the general solution of the ODE is 4. unstable spiral
aX 1 bX 2 e , where a, b are arbitrary (CSIR NET June 2015)
Page 3
t 1 t
et e e
4. and 2
0
t
e
7. Consider the system of differential equations
dx
2x 7 y
dt
dy
3x 8 y
dt
Then the critical point 0,0 of the system
is an
1. asymptotically stable node
2. unstable node
3. asymptotically stable spiral
4. unstable spiral
(CSIR NET June 2018)
Answers
----------------------- S C Q --------------------
1. 1 2. 3 3. 1
4. 3 5. 4 6. 1
---------------------- M C Q --------------------
1. 1,2,3,4 2. 1 3. 1
4. 1 5. 1,3 6. 3,4
7. 1
J.R. INSTITUTE OF MATHEMATICS
189/35 BEHIND RAILWAY STATION, VAISH COLLEGE ROAD, ROHTAK PIN-124001 (HARYANA)
E-mail us on - [email protected], [email protected] Mob. 8607383607, 9802177766
Page 1
Chapter - 7
Existence and uniqueness problems
Def. Lipschitz condition : A function f x, y defined in domain D is said to satisfy the Lipschitz
condition w.r.t. y in D if, these exists constant K such that f x, y1 f x, y2 K y1 y2 , for all
exists and bounded in D then f x, y satisfy the Lipschitz condition with the Lipschitz constant
f
max , x, y D .
y
f
Remark: Boundness of is only a sufficient condition for the Lip. Condition but not necessary i.e.
y
f
If is not bounded then nothing can be said about Lipschitz condition.
y
Theorem (2) : CauchyPeano Existence theorem :
Hypothesis : Consider the initial value problem
dy
f x, y ......(1)
dx
and y x0 y0 ......(2)
R x, y : x x0 a, y y0 b
b
(ii) This solution x is defined at least on the interval x0 h, x0 h , where h min a, and
M
f x, y M for all x R
2
b
(ii) This solution x is defined atleast on the interval x0 h, x0 h , where h min a, and
M
f x, y M for all x R
Remark : If the function f x, y is not C , Lip then nothing can be said about existence and
Exercise 7.1
1. Which of the following functions satisfies a Lipschitz condition in the rectangle R defined by
R x, y : x a, y b, a 0, b 0
2 2 3
(vii) f x, y y 3 (viii) f x, y y 5 (ix) f x, y y 2
5 2 3
(x) f x, y y 2 (xi) f x, y xy 3 (xii) f x, y xy 2
(xxiv) f x, y x y
J.R. INSTITUTE OF MATHEMATICS
189/35 BEHIND RAILWAY STATION, VAISH COLLEGE ROAD, ROHTAK PIN-124001 (HARYANA)
E-mail us on - [email protected], [email protected] Mob. 8607383607, 9802177766
Page 3
2. Which of the following functions satisfies a Lipschitz condition in the rectangle R defined by
R x, y : x 3 1, y 2.
(xxiv) f x, y x y
3. Which of the following functions satisfies a Lipschitz condition in the rectangle R defined by
R x, y : x 3 1, y 3 2.
(xxiv) f x, y x y
4. Discuss the existence and uniqueness of the following initial value problems. Also find the
solution, wherever possible:
4 4 4
dy dy dy
(i) y 3 , y 0 0 (ii) y 3 , y 0 1 (iii) y 3 , y 1 0
dx dx dx
4 2 2
dy dy dy
(iv) y 3 , y 1 1 (v) y 3 , y 0 0 (vi) y 3 , y 0 1
dx dx dx
2 2
dy dy dy
(vii) y 3 , y 1 0 (viii) y 3 , y 1 1 (ix) y , y 0 0
dx dx dx
dy dy dy
(x) y , y 0 1 (xi) y , y 1 0 (xii) y , y 1 1
dx dx dx
dy dy
(xiii) y 2 , y 0 1 [ NET June 2013 ] (xiv) 1 x 2 y, y 0 1
dx dx
dy
(xv) 1 f 2 x y , y 0 1, where f is continuous on R [ NET DEC 2011 ]
dx
dy dy 2 y
(xvi) 1 y 2 , y 0 0 (xvii) e , y 0 0
dx dx
1
dy dy
(xviii) xy 3 , y 0 0 [ NET June 2013 ] (xix) y , y 1 0
dx dx
2
dy dy
(xx) 60 y 5 , y 0 0 [ NET Dec. 2012 ] (xxi) x 3 y , y 1 1
dx dx
dy
(xxii) f x y x , y 0 1, where f is a continuous function on R [ NET June 2012 ]
dx
5. Find the largest possible interval for x in which Picard’s uniqueness theorem (PUT) guarantees the
existence of a unique solution. Also, wherever possible, find the interval for x in which the
solution actually exist.
dy dy 2 y dy 2
(i) 1 y 2 ; y 0 0 (ii) e , y 0 0 (iii) x y 2 , y 0 0
dx dx dx
dy dy
(iv) 16 y 2 , y 0 0 (v) y , y 0 1, R : x 1, y 1 1
dx dx
J.R. INSTITUTE OF MATHEMATICS
189/35 BEHIND RAILWAY STATION, VAISH COLLEGE ROAD, ROHTAK PIN-124001 (HARYANA)
E-mail us on - [email protected], [email protected] Mob. 8607383607, 9802177766
Page 5
dy 1 2 x 3 y
(vi) , y 0 0, R : x 2, y 1
dx 2 x 2 y 2
dy 1 x y
(vii) , y 1 1, R : x 1 1, y 1 1
dx 2 x 2 y 2
dy 2 xy
(viii) , y 1 1, R : x 1 1, y 1 1
dx 4 x 2 y 2
dy 3xy
(ix) , y 0 0, R : x 1, y 2
dx 2 cos xy
Answers
1. (i) Yes (ii) Yes (iii) Yes (iv) Yes (v) Yes
(vi) No (vii) No (viii) No (ix) Yes (x) Yes
(xi) No (xii) Yes (xiii) No (xiv) No (xv) No
(xvi) No (xvii) No (xviii) Yes (xix) Yes (xx) Yes
(xxi) Yes (xxii) Yes (xxiii) Yes (xxiv) Yes
2. (i) Yes (ii) Yes (iii) Yes (iv) Yes (v) Yes
(vi) No (vii) No (viii) No (ix) Yes (x) Yes
(xi) No (xii) Yes (xiii) No (xiv) No (xv) No
(xvi) No (xvii) No (xviii) Yes (xix) Yes (xx) Yes
(xxi) Yes (xxii) Yes (xxiii) Yes (xxiv) Yes
3. (i) Yes (ii) Yes (iii) Yes (iv) Yes (v) Yes
(vi) Yes (vii) Yes (viii) Yes (ix) Yes (x) Yes
(xi) Yes (xii) Yes (xiii) Yes (xiv) Yes (xv) Yes
(xvi) Yes (xvii) Yes (xviii) Yes (xix) Yes (xx) Yes
(xxi) Yes (xxii) Yes (xxiii) Yes (xxiv) Yes
4. (i) trivial solution , no actual solution, unique solution
1
3
(ii) y 3 , non trivial solution, unique solution
3 x
(iii) trivial solution, no actual solution, unique solution
1
1
3 y03
(iv) y3 1 2
, non trivial solution, unique solution
y03 x 3 x0 y03
6
1
x
(v) trivial solution, y3 , infinite solution
3
1
x3
(vi) y3 , non trivial solution
3
1
x 1
(vii) trivial solution, y 3 , infinite solution
3
1 1
x x0
(viii) y3 y03 , non trivial solution
3
1
x
(ix) trivial solution, y2 , infinite solution
2
1
x2
(x) y 2 , non trivial solution, unique solution
2
1
x 1
(xi) trivial solution, y2 , infinite solution
2
1 1
x x0
(xii) y2 y02 , non trivial solution
2
1
(xiii) y , non trivial solution
1 x
x3
x
(xiv) y e 3 , non trivial solution, unique solution
x f 2 ( x ) dx
(xv) y e , non trivial solution
(xvi) y tan x , non trivial solution, unique solution
1
(xvii) y log(1 2 x ) , non trivial solution, unique solution
2
2
x2
(xviii) y3 , trivial solution, infinite solution
3
1
x 1
(xix) y2 , trivial solution, infinite solution
2
3
(xx) y 5 36 x , trivial solution, infinite solution
Page 7
f ( x ) dx
(xxii) y e , non trivial solution, unique solution
J.R. INSTITUTE OF MATHEMATICS
189/35 BEHIND RAILWAY STATION, VAISH COLLEGE ROAD, ROHTAK PIN-124001 (HARYANA)
E-mail us on - [email protected], [email protected] Mob. 8607383607, 9802177766
Page 1
Page 3
0 for 0 t 10
1. If
a x dx , then y is bounded
4. u t 4
0
t 10 for t 10
(CSIR NET Dec 2015) 2. If
a x dx , then lim y x exists
0
x
Answers
----------------------- S C Q --------------------
1. 1 2. 2 3. 1
4. 4 5. 4 6. 4
7. 2 8. 1 9. 2
10. 3 11. 2 12. 3
13. 3
---------------------- M C Q --------------------
1. 4 2. 1,4 3. 2,3
4. 1,2,3,4 5. 1,2,3,4 6. 2,3
7. 2,4 8. 1,2,3,4
J.R. INSTITUTE OF MATHEMATICS
189/35 BEHIND RAILWAY STATION, VAISH COLLEGE ROAD, ROHTAK PIN-124001 (HARYANA)
E-mail us on - [email protected], [email protected] Mob. 8607383607, 9802177766
Page 1
Chapter - 8
Equations of First Order and Higher Degree and singular solution
Introduction : In this chapter, we shall study those differential equations, which are of first order but
not of first degree. The general form of such differential equation is :
n n1 n 2
dy dy dy dy
P0 P1 P2 .... Pn 1 Pn 0 ......(1)
dx dx dx dx
where P0 , P1 , P2 ,................, Pn1 , Pn are functions of x and y.
dy
Usually , we denote by p, so that equation (1) takes the form
dx
P0 p n P1 p n1 P2 p n 2 ..... Pn 1 p Pn 0
We shall solve such equations and find the general and singular solutions of these equations.
The differential equation of first order but not of first degree can also be expressed as
dy
f (x , y, p) = 0, where p = .
dx
To solve such differential equations, we shall consider the following five cases.
(i) Equations solvable for p
(ii) Equations solvable for x
(iii) Equations solvable for y
(iv) Lagrange’s equation and Clairaut’s equation.
(v) Equations reducible to Clairaut’s equation .
We shall also find the singular solution of such equation. Let us discuss these cases one by one.
1. Equations solvable for p :
A differential equation of first order and nth degree is of the form
P0 p n P1 p n 1 P2 p n 2 ........ Pn 1 p Pn 0 ......(1)
where P0 , P1 , P2 ,......, Pn 1 , Pn are functions of x and y. Let left hand side of (1) be solvable for p, then it
dy dy dy
or f1 ( x , y ) , f 2 ( x , y ) , ...., f n ( x , y)
dx dx dx
2
These are differential equations of first order and first degree and can be solved easily. Let the
solutions of these equations be
F1 ( x , y , c1 ) 0, F2 ( x , y , c2 ) 0,............, Fn ( x , y , cn ) 0
Since the given differential equation is of first order ; therefore , the general solution can not have
more than one arbitrary constant , so we take
c1 c2 c3 ..............cn c (say)
Working Rule :
1. Reduce the given differential equation in linear factors of p
2. Solve individually these linear factors of p
3. Multiply all the solutions obtained in step 2 and equate it to zero. This gives the general solution of
the given differential equation.
Example 1 : Solve p2 7p + 12 = 0
Solution : The given differential equation is
p2 7p + 12 = 0. ......(1)
or (p 3) (p 4) = 0
p=3 or p=4
dy dy
i.e. 3 or 4
dx dx
dy = 3 dx dy = 4 dx
Integrating , y = 3x + c Integrating , y = 4x + c
or y 3x c = 0 or y 4x c = 0
Therefore the general solution of given equation is y 3x c y 4 x c 0.
2. Equations Solvable for x :
A differential equation solvable for x can be put in the form
x f y, p .......(1)
Differentiating w. r. t. y , we get
dx dp
g y, p,
dy dy
J.R. INSTITUTE OF MATHEMATICS
189/35 BEHIND RAILWAY STATION, VAISH COLLEGE ROAD, ROHTAK PIN-124001 (HARYANA)
E-mail us on - [email protected], [email protected] Mob. 8607383607, 9802177766
Page 3
1 dp
g y , p, .......(2)
p dy
This equation involves two variables y and p. Let solution of this equation be
y , p, c 0 .......(3)
Elimination of p from (1) and (3) gives the required solution. If elimination of p is not possible, then
expressions of x and y in terms of p together form the solution of the given equation.
Working Rule :
1. Express x explicitly in terms of p and y.
dx 1
2. Differentiate the equation w.r.t. y and replace by . The equation obtained involves p and y.
dy p
3. Solve the equation obtained in step 2 and eliminate p from this solution and the given equation.
Example 2 : Solve y 2px = y2p3
Solution : The given differential equation is y = 2px + y2p3 ......(1)
y y 2 p2
x
2p 2
Differentiating w.r.t. y, we have
dx 1 y dp dp
2 yp 2 py 2
dy 2 p 2 p dy dy
1 1 y dp dp
or 2 yp 2 py 2
p 2 p 2 p dy dy
1 y dp
or yp 2 2 py 2
2p 2p dy
1 y 1 dp
or yp 2 p2 y
2p p 2p dy
1 y 1 dp
or yp 2 yp 2 0
2p p 2p dy
1 y dp
or yp 2 1 0
2p p dy
y dp
or 1 0
p dy
4
dy dp
or 0
y p
Integrating , we get log y log p log c
or log yp log c
c
i.e. yp c p
y
Thus from equation (1), we have
2c c3
y x y2 3
y y
or y 2 2cx c 3
which is the general solution of given equation.
3. Equations Solvable for y :
A differential equation solvable for y can be put in the form
y f x, p .......(1)
dp
p g x, p, .......(2)
dx
This equation involves two variables x and p. Let solution of this equation be
x, p, c 0 .......(3)
Elimination of p from (1) and (3) gives the required solution. If elimination of p is not possible, then
expressions of x and y in terms of p together form the solution of the given equation.
Working Rule :
1. Express y explicitly in terms of p and x.
dy
2. Differentiate the equation w.r.t. x and replace by p. The obtained equation involves p and x.
dx
3. Solve the equation obtained in step 2 and eliminate p from this solution and the given equation.
Example 3 : Solve y = 3x + log p
Solution : The given differential equation is
y = 3x + log p ......(1)
Differentiating w.r.t. x, we have
J.R. INSTITUTE OF MATHEMATICS
189/35 BEHIND RAILWAY STATION, VAISH COLLEGE ROAD, ROHTAK PIN-124001 (HARYANA)
E-mail us on - [email protected], [email protected] Mob. 8607383607, 9802177766
Page 5
dy 1 dp
= 3+
dx p dx
1 dp
or p = 3
p dx
1 dp
or p3 =
p dx
dp
or dx =
p( p 3)
Integrating both sides, we have
1 1 1
x= dp c [Making partial fractions]
3 p p 3
1
or x= log p log p 3 c
3
p
or 3x = log 3c
p3
p
or 3x = log log c 3c log c
p 3
p 3
or log e3x = log c x log e x
p
p 3
or e3x = c
p
3
or e3x = c 1
p
3 3 1
or 1 ce3x =
p
or p
1 ce3x c c
3
Then from equation (1), we have y = 3 x log which is the required solution.
1 ce3 x
6
Exercise 8.1
Solve the following differential equations :
2
dy dy
1. 2
p 5p 6 0 2. x xy 6 y 2 0
2
dx dx
3. y 2 px p 2 y 4. x y p2
5. y px p 2 x 4 6. y p sin p cos p
Answers
1. ( y 3x c)( y 2 x c) 0. 2. ( x 3 y c)( y cx 2 ) 0.
e.g. y p 2 x sin p
dx p p
x
dp p p p p
dx p p
x
dp p p p p
J.R. INSTITUTE OF MATHEMATICS
189/35 BEHIND RAILWAY STATION, VAISH COLLEGE ROAD, ROHTAK PIN-124001 (HARYANA)
E-mail us on - [email protected], [email protected] Mob. 8607383607, 9802177766
Page 7
This is a linear differential equation in terms of x and p. Let the solution of this equation be
f x, p, c 0 .......(2)
Then (1) and (2) taken together provide the required solution.
dy
Differentiate this equation w.r.t. x and replace by p. We get a linear differential equation in p and
dx
x. Solve it by the method of integrating factor. This solution taken together with the given equation
will provide the required solution.
Example 4 : Solve y xp 2 p.
( p 1)2 x p log p c
8
p log p c
or x ......(3)
( p 1)2
Putting the value of x in (1), we have
p3 p 2 log p cp 2
y p ......(4)
( p 1) 2
Equation (3) and equation (4) taken together provide the required solution.
5. Clairaut’s equation (a particular case of Lagrange’s equation ) :
In the Lagrange’s equation
y x p p
y xp p
y = p x sin1p
which is Clairaut’s equation.
The solution is given by y = cx sin1 c.
J.R. INSTITUTE OF MATHEMATICS
189/35 BEHIND RAILWAY STATION, VAISH COLLEGE ROAD, ROHTAK PIN-124001 (HARYANA)
E-mail us on - [email protected], [email protected] Mob. 8607383607, 9802177766
Page 9
or e y ce x c 5
Exercise 8.2
Solve the following differential equations :
a
1. y (1 p ) x p 2 2. y xp
p
Answers
a
1. y 2 p 2 ce p (1 p ) and together with given relation. 2. y = cx +
c
a 2c c2
3. y 2 cx 2 4. x2 + y2 = c(x + y)
c 1 4
---------------------------------------------------------------------------------------------------------------------------
Singular solution : A singular solution in this stronger sense is often given as tangent to every
solution from a family of solutions. A singular solution is the envelope of the family of solutions
cannot be obtained from general solution by giving particular values of arbitrary constants.
The singular solution is related to the general solution by its being what is called the envelope
of that family of curves representing the general solution. An envelope is defined as the curve that is
tangent to a given family of curves.
Envelope : A curve which touches each member of a one-parameter family of curves and at each
point is touched by some member of the family, is called the envelope of the that parameter family of
curves.
Example 1 : y ' 4 has the general solution y ( x c )2 , which is a family of parabolas. The line
y ( x) 0 is also a solution of the differential equation, but it is not a member of the family constituting
the general solution. Hence line y ( x) 0 is singular solution of the differential equation.
Example 2 : Consider family of circles with centre lying on x-axis and radius 3 is. This family is
given by ( x a )2 y 2 9 . And corresponding differential equation can be obtained by eliminating a
parameter a.
( x a ) yy ' ( yy ') 2 y 2 9
y 2 y '2 1 9
also we can easily verify that y ( x) 3 and y ( x) 3 are the solution of y 2 y '2 1 9 which cannot
be obtained from general solution. These both are envelope to the general solution. Hence y ( x) 3
and y ( x) 3 are singular solution of the differential equation.
Def. Singular Solution : A singular solution is a solution of the differential equation which does not
contain any arbitrary constant but can not be derived from general solution.
In other words, it is not obtained just by providing some particular values to arbitrary constants in
general solution.
J.R. INSTITUTE OF MATHEMATICS
189/35 BEHIND RAILWAY STATION, VAISH COLLEGE ROAD, ROHTAK PIN-124001 (HARYANA)
E-mail us on - [email protected], [email protected] Mob. 8607383607, 9802177766
Page 11
If (3) satisfies the given equation, it is the required singular solution. If (3) does not satisfy (1), then
resolve F x, y into simpler factors.
Those factors which satisfy the equation (1), constitute the singular solution.
Method II : Let f (x, y, p) = 0 ......(1)
be the given differential equation and general solution of this equation be
(x, y, c) = 0 ......(2)
Differentiating (2) partially w. r. t. ‘c’, we get
0 ......(3)
c
Eliminating c between (2) and (3), we get c-discriminant relation, let it be
F x, y 0 ......(4)
12
If (4) satisfies the given equation, it is the required singular solution. If (4) does not satisfy (1), then
resolve it into simple factors.
Those factors which satisfy the equation (1) constitute the singular solution.
Remark : (1) If the given equation is quadratic in p, say
Ap 2 Bp C 0
then p-discriminant can be directly taken as
B 2 4 AC 0
(2) Similarly, if the general solution is quadratic in c, say
Ac 2 Bc C 0
then c-discriminant can be directly taken as
B 2 4 AC 0
a
Example 3 : Solve the differential equation y px and obtain its singular solution.
p
a
Solution : The given differential equation is y px .......(1)
p
which is Clairaut’s equation.
Thus its general solution is
a
y cx [Replacing p by c]
c
c 2 x cy a 0
which is quadratic in c.
Hence c-discriminant is given by
y 2 4ax 0
y2 = 4 ax
This is the required singular solution.
(It can be checked that it satisfies the given differential equation)
Exercise 8.3
Find the general and singular solution of the following differential equations :
1. 4 xp 2 (3x a) 2 0 2. y 2 2 pxy p 2 ( x 2 1) m 2
3. p 2 ( x 2 a 2 ) 2 pxy y 2 b 2 0 4. xp 2 2 yp x 2 y 0
J.R. INSTITUTE OF MATHEMATICS
189/35 BEHIND RAILWAY STATION, VAISH COLLEGE ROAD, ROHTAK PIN-124001 (HARYANA)
E-mail us on - [email protected], [email protected] Mob. 8607383607, 9802177766
Page 13
Answers
1. ( y c)2 x ( x a )2 ; x 0 is singular solution
x2 y 2
3. ( y cx)2 a 2c 2 b 2 , 1 is singular solution.
a 2 b2
1
4. y x x 2c ; y x 2 x is singular solution.
2c
---------------------------------------------------------------------------------------------------------------------------
Analytic function : If a function f x is single valued in a domain D, then its derivative f ' a
f a h f a
exists if the limit given by f ' a lim exists and is independent of the path along
h 0 h
which h 0 . This function f x is said to be analytic at x a if there exists a neighbourhood
x a , 0, at all points of which f ' x exists. If f ' x exists at every point in the domain D
except some exceptional points then f x is said to be analytic in D. These exceptional points are
Ordinary and singular points : Let us consider a differential equation of the form
y " P x y ' Q x y 0 then a point x x0 is called an ordinary point of the equation if the functions
If the point x x0 is not ordinary point of the equation then this is the singular point.
Example 1. x 2
1 y " xy ' y 0
x y
Convert into general form y " y ' 2 0
x 1 x 1
2
14
x 1
P x , Q x
x 1 x 1 x 1 x 1
Now both P x and Q x are analytic at x 0 then x 0 is the ordinary point. The function P x
x 2 x 1
Then x 1 is the singular point. Also x 1 P x and x 1 Q x
x 1 x 1
showing both x 1 P x and x 2 1 Q x are analytic at x 1 is a regular singular point.
7 x 1 3
Convert into general form y " 2 y0
2x 2x
7 x 1 3
P x , Q x 2
2x 2x
We can see that both P x and Q x are not analytic at x 0 so x 0 is a singular point. Also
7 x 1 2 3
x 0 P x and x 0 Q x which shows that both x 0 P x and
2 2
2
x 0 Q x are analytic at x 0 . Then x 0 is the regular singular point of the given differential
equation.
Exercise 8.4
1. (i) x 3 x 1 y " 2 x 1 y ' 3 xy 0 (ii) 3 x 1 xy " x 1 y ' 2 y 0
Answers
1
1. (i) x 0; irregular, x 1, regular (ii) x 0, regular
3
(iii) x 1, 1, regular (iv) x 0, irregular
2. (i) ordinary point
(ii) regular point
(iii) regular singular point
J.R. INSTITUTE OF MATHEMATICS
189/35 BEHIND RAILWAY STATION, VAISH COLLEGE ROAD, ROHTAK PIN-124001 (HARYANA)
E-mail us on - [email protected], [email protected] Mob. 8607383607, 9802177766
Page 1
3. x 2 6 y 0 4. x 2 6 0
1. x 2 0 2. y 2 0
3. x 4 0 4. x 2 0
6. Consider the differential equation
1
x 1 y " xy ' y 0 . Then
x
1. x 1 is the only singular point