Survey
* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project
* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project
Plateau principle wikipedia , lookup
Mathematical optimization wikipedia , lookup
Generalized linear model wikipedia , lookup
Inverse problem wikipedia , lookup
Relativistic quantum mechanics wikipedia , lookup
Mathematical descriptions of the electromagnetic field wikipedia , lookup
Perturbation theory wikipedia , lookup
Simplex algorithm wikipedia , lookup
Navier–Stokes equations wikipedia , lookup
Routhian mechanics wikipedia , lookup
M.Sc. (Previous) Mathematics Paper –V Differential Equations BLOCK- I UNIT –I : Homogeneous linear differential equations with variable coefficients, Simultaneous differential equations and Total differential equations. UNIT -II: Picard’s method of integration, Successive approximation, Existence and uniqueness theorem. 1 BLOCK-INTRODUCTION The study of differential equations is a wide field in pure and applied mathematics, physics, meteorology, and engineering. All of these disciplines are concerned with the properties of differential equations of various types. Pure mathematics focuses on the existence and uniqueness of solutions, while applied mathematics emphasizes the rigorous justification of the methods for approximating solutions. Differential equations play an important role in modeling virtually every physical, technical, or biological process, from celestial motion, to bridge design, to interactions between neurons. This block contains two units: Unit I: This unit deals the concept, different methods to solve the Homogeneous linear differential equations with variable coefficients, Simultaneous differential equations and Total differential equations with several solved examples followed by exercise to check the progress of reader. Unit II: This unit deals the concept of numerical problems and their solutions especially in reference to Picard’s method of integration, Successive approximation, Existence and uniqueness theorem with several solved examples followed by exercise to evaluate reader by himself. At the end a list of reference books are given for the convenience to the reader. 2 UNIT –I : Homogeneous linear differential equations with variable coefficients, Simultaneous differential equations and Total differential equations. Unit Structure : 1.0 Object 1.1 Introduction 1.2 Linear differential equations and its kinds 1.3 Homogeneous Linear differential equations with variable coefficients 1.4 Algorithm to solve Homogeneous Linear differential equations with variable coefficients with solved examples 1.5 Equations reducible to Homogeneous Linear differential equations 1.6 Simultaneous differential equations 1.7 Algorithm to solve Simultaneous differential equations by different methods with solved examples 1.8 Total differential equations 1.9 Algorithm to solve Total differential equations by different methods with solved examples 1.10 Unit brief review (Summary) 3 1.1 Object: At the end of this unit students are in a position to find the solutions of Homogeneous Linear differential equations with variable coefficients, Simultaneous differential equations and Total differential equations in easy manner. 1.1 Introduction: Present chapter is deal with the study of Homogeneous Linear differential equations with variable coefficients, Simultaneous differential equations and Total differential equations. For the sack of convenience to students a brief summary related to the topic viz. linear differential equations, its kind and method to find the solutions in a simple manner is given. 1.2 Linear differential equations and its kind: 1.2.1 Linear differential equations: A differential equation of the form =X + where ………………. (1) and X are function of x or constants. In terms of D- notations (1) can be written as f(D)y = X or ………………(2) + where D = , , ……………. and so on. The general solution of equation (1) or (2) is given by y= complementary function (C.F.)+ particular integral (P.I.) Things to remember: (i) If X=0 in equation (1), the general solution is C.F. itself. 4 (ii) If all are constant in (1) then equation (1) is called linear differential equations with constant coefficients or otherwise it is known as linear differential equations with variable coefficients. 1.2.2 Algorithm to find complementary function (C.F.): Step I: find auxiliary equation (A.E.) f(m)=0 by writing D=m in f(D) of equation (2). Step II: find the roots of the A.E. i.e. values of m. Let the roots are m1, m2 ,…… , mn . Step III: required C.F. is obtained as per the roots stated below: Roots of A.E. Complementary function (C.F.) All roots m1 , m2 ,…… , mn are real and +……+ + different. m1 = m2 , but other roots are real and +……+ + different. If roots are imaginary If ( ), ( cos (say) ) repeated twice. sin ) or ) or ) corresponding part of C.F. is [ Solved examples: Example (1): Solve . Solution: Here A.E. is - 3m - 4 = 0 (m- 4)(m+1) = 0 Hence roots are 4 and -1, real and different. Therefore C.F. is y = general solution. + ,itself 5 Example (2): Solve . Solution: Here A.E. is =0 (m+ 2)(m- 2)2 = 0 Hence roots are 2,2 and -2, two are equal and one is different. Therefore C.F. is y= + , itself general solution. Example (3): Solve . Solution: Here A.E. is =0 (m- 2)(m2 + 2m + 4) = 0 Hence roots are 2, and -1 C.F. is y= + m = -1 , one is real and other is a pair of imaginary. Therefore cos sin 1.2.3 Particular integral (P.I.): If X ), itself general solution. 0 in equation (1) then P.I. = . In previous classes we studied the method to find P.I. either by resolving the f(D) into linear factors and partials fractions and then each factor of the form ∫ solve by ………… (3) or by shortcut methods given below: Type of function X X= then put D= a in f(D) i.e. Corresponding P.I. .If then (D-a) is one of the factor of f(D).This 6 factor is solve by formula (3) and rest is solve by shortcut method given here. X = sin ax (or cos ax) then put D2= -a2 in sin ax (or cos ax), provided f(D) i.e. or otherwise use following formula: sin ax = - or cos ax = X = xm then P.I. is [ Expand [ using binomial expansions and solve corresponding terms treating D as differentiation. X= , , solve V as per format given above. Solved examples: Example (1): Solve . Solution: Here A.E. is +4m +3 = 0 (m+3)(m+1) = 0 Hence roots are -3 and -1, real and different. Therefore C.F. is y = Now A.T.Q. Since f(a) P.I. = [here X= + ] on putting D=-2, we get P.I. y = =- =- Hence the general solution is y = C.F.+ P.I. i.e. y= + - . 7 Example (1): Solve . Solution: Here A.E. is -4 = 0 (m+2)(m-2) = 0 Hence roots are 2 and -2, real and different. Therefore C.F. is y = Now A.T.Q. P.I. = Therefore P.I. = . + [we know that cos 2x= 2 y = [ -1.] ] y = y = y = y = = Hence the general solution is y = C.F.+ P.I. i.e. y= + . Example (3): Find particular integral (P.I.) Solution: Here f(D) = = (D- 2)(D- 3) Therefore P.I. = = = . =( ) = = 8 = [since (1-D)-1 = 1+D+D2+…….] = [ D means differentiation wrt. x] =( ) = + . Check your progress: Solve the following differential equation: Q.1 . Q.2 [Ans. y = + + . [Ans. y = Q.3 . [Ans. y = ] + + + ] + ] Q.4 [Ans. y = Q.5 + + + ] [Ans. y = ]. 1.3 Homogeneous Linear differential equations with variable coefficients 1.3.1 Homogeneous Linear differential equations with variable coefficients: An equation of the form + =X …………. (1) + where are constant and X is function of x or constant, is called Homogeneous Linear differential equations. 9 Note: (i) Equation (1) above is called homogeneous because in each term the power of x is same as the order of the derivative in that term. (ii) Since the coefficient in each term is variable, therefore it is known as Homogeneous Linear differential equations with variable coefficients. (iii) This equation is also known as Cauchy’s linear equation. (iv) By the substitution of x = or z = log x, the above equation (1) is transferred into the linear differential equation with constant coefficient changing the independent variable x to z as below: If x = or z = log x then we have = ………… (2) = or ..……… (3) (say), ..………. (4) Now Hence ( )= ( Similarly we can obtain ( )= ( ) )= = ( ) , = and so on. On putting the equivalent values of etc. in equation (1),it reduces to linear differential equation with constant coefficients having independent variable z. This equation can be solved easily by the method given in (1.2). Important Note: Notation D stands for and . 1.4 Algorithm to solve Homogeneous Linear differential equations with variable coefficients with solved examples Algorithm: 10 Step I: Write the given equation in D-notation form. etc in the equation by Step II: Replace = , = and so on. Step III: Obtain equations is linear differential equation with constant coefficients, find C.F. and P.I. as per the method given in article 1.2 treating z as independent variable. Step IV: Lastly put back z = log x to get the required result. Solved examples : Example 1: Solve - . Solution: D-notation form of the given equation is On putting x = = and [ Here A.E. is ] , we have = 2z (as log ) = 0 gives m = 1,1. So C.F. is y = [ Now P.I. = = [1+2 +……….] = Hence the general solution is y = C.F.+ P.I. i.e. or y= y= . ( putting z = log x) Ans. 11 Example 2: Solve - . Solution: D-notation form of the given equation is On putting x = = and [ Here A.E. is , we have ]y=2 (as ) = 0 gives m = 2, 2. So C.F. is y = Now P.I. = = P.I. = (here X is of the form = = ∬ (here , V=1,art.1.2) ) = = Hence the general solution is y = C.F.+ P.I. i.e. y= or y= or y= Example 3: Solve . . + (putting z = log x) Ans. . Solution: D-notation form of the given equation is . 12 On putting x = = and , = we have [ ]y= Here A.E. is . =0 or Hence, C.F. = . Now P.I. = = First part of P.I. = =5 [On putting ] ………… (i) Now second part of P.I. = = On putting because it becomes zero, a failure case of f(a)=0 ∫ which can be solved by the formula , = = =2 P.I. = 5 ∫ ……….. (ii) =2 from (i) and (ii) =2 +2 Hence the general solution is y = C.F. + P.I. 13 i.e. y = or y= . (putting z = log x) ( y= ) Ans Check your progress: Solve the following differential equation: Q.1 + . [Ans. y = Q.2 [Ans. y = Q.3 [Ans. y = Q.4 + [Ans. y = + + - + ] ] + ] + + ] (Hint: Divide whole equation by x) Q.5 [Ans. y = + + ]. 1.5 Equations reducible to Homogeneous Linear differential equations An equation of the form + =X + …………. (1) where are constant and X is function of x or constant, can be reduced into Homogeneous Linear differential equations by putting 14 then we get b= =b and so we have ( ) , ( ) …………etc. On putting the above values in (1), equation reduces to Homogeneous Linear differential equations and easily solve by the same technique as per article (1.4). Solved Examples: Example 1: Solution: on putting , then ]y= , the equation becomes . , where , A.E. is , Hence, C.F. = Now P.I. = = = , (on writing ) = Hence the general solution is y = C.F. + P.I. i.e. y = or y= (back substitution ) Ans. 15 Example 2: Solution: on putting , then [ ]y= A.E. is Hence, , the equation becomes . , where , C.F. = Now P.I. = = = , (on writing 4 ) =2z , failure case) (since cos ax = Hence the general solution is y = C.F. + P.I. i.e. y = +2z y= +2 (back substitution ) Ans. Check your progress: Solve the following differential equation: Q.1 [Ans. y = [ ] Q.2 [Ans. y = }. 16 1.6 Simultaneous differential equations System of differential equations in which number of differential equations will be same as is the number of dependent variables, called the system of Simultaneous differential equations and the differential equations are called Simultaneous differential equations. Present chapter deal with two types of Simultaneous differential equations. 1.7 Algorithm to solve Simultaneous differential equations by different methods with solved examples 1.7.1 Type I: Simultaneous differential equations of first order and of the first degree with constant coefficients i.e. …… (1) …… (2) and Now to solve the above equations we apply elimination method. We first eliminate x (or y) with suitable operation and obtain linear differential equation with constant coefficients in y and t (or in x and t) which can be solved by the method, discussed in article (1.2). Then we find the value of other variable by substituting value of first variable in (1) or (2). Important note: The number of arbitrary constants in general solution is equal to the degree of D in determinant | |=0 Solved Examples: Example 1: Solve . Solution: On writing given equations in D-notation form, we have 17 (D- 7) x + y = 0 …. (i) -2x + (D- 5) y = 0 …. (ii) To eliminate y we operate (i) by (D- 5 ) and then subtracting, we get (D- 7) x + 2x=0 i.e. (D2-12D+37) x = 0 A.E. (m2-12m+37) = 0 , m=6 x = e6t(c1cos t +c2sin t) On substituting the above value of x in (i), we get y = -(D- 7) e6t(c1cos t +c2sin t) =7e6t (c1cos t +c2sin t)-D e6t(c1cos t +c2sin t) =7e6t (c1cos t +c2sin t)- 6e6t (c1cos t +c2sin t)- e6t(-c1 sin t +c2 cost) = e6t (c1cos t +c2sin t) - e6t (-c1 sin t +c2 cost) = e6t [ (c1 - c2 ) cost + (c1+ c2 ) sin t] Hence the general solution is x = e6t (c1cos t +c2sin t), y= e6t [ (c1 - c2 ) cost + (c1+ c2 ) sin t] . Example 2: Solve . Solution: On writing given equations in D-notation form, we have (D+5) x + y = …. (i) -x + (D+3) y = …. (ii) To eliminate x we operate (ii) by (D+ 5 ) and then adding, we get (D+5)(D+3)y + y =(D+5) i.e. (D2+8D+15) y = D 18 (D2+8D+15) y = 2 (D2+8D+15) y = A.E. (m2+8m+15) = 0 , m= -3, -5 C.F. y = (c1 e-3t +c2 e-5t) and P.I. y= ( =( ) = ( ) ( ) {as per method given in art.1.2} ) = y = (c1 e-3t +c2 e-5t)+ On substituting the above value of y in (ii), we get x = (D+3) { (c1 e-3t +c2 e-5t)+ =D{(c1 e-3t +c2 e-5t)+ = -3c1 e-3t -5c2 e-5t + }+3{(c1 e-3t +c2 e-5t)+ +3c1 e-3t +3c2 e-5t+ }- - = -2c2 e-5t + Hence the general solution is x = -2c2 e-5t + , y = (c1 e-3t +c2 e-5t) + . 1.7.2 Type II: Simultaneous differential equations of first order and of the first degree in the derivatives. Here we are giving the technique to solve the equations containing three variables. The general form of Simultaneous differential equations of first order and of the first degree containing three variables is 19 ……… (1) P1 dx + Q1 dy + R1 dz = 0, P2 dx + Q2 dy + R2 dz = 0 where the coefficients are function of x, y, z. Solving these equations simultaneously, we get, ……… (2) or Thus simultaneous equations (1) can always be put in the form (2). Methods to solve 1.7.3 First method: Suppose that any two fractions, directly integrated then after integration we find an integral. Same procedure we apply for other two fractions. The two integrals so obtained form the complete solution. Sometimes first integral may be used to simplify the other two fractions. Solved Examples: Example1: Solve . Solution: On taking first two fractions, we have =y , on integration, we get Similarly, taking second and third fractions, we get, Thus the complete solution is Example 2: Solve = c, + c or ,gives =c c’ c’ . Solution: On taking first two fractions, we have =y , on integration, we get + c or = c ...... (i) 20 Similarly, taking second and third fractions, ( cy ) on integration, we get c = , c’ or c’ ii (i) and (ii) form the complete solution. Note: It will be noted that this technique to solve the equations applied to Simultaneous differential equations containing any number of variables. 1.7.4 Second Method: We have . If l, m, n are such that then we get , . If it is exact differential equation (say) du = 0, then u=c is one part of the complete solution. Here l, m, n are known as multipliers. This method may be repeated to get another integral by choosing new multipliers l’, m’, n’. Note: Sometimes one integral we find by using method first and second integral using second method. Solved Examples: Example1: Solve . Solution: On choosing 1, 1, 1 as multipliers, we get or , gives x+ y+ z = c ………. (i) Again, choosing x, y, z as multipliers, we get , gives c‘ ii 21 (i) and (ii) form the complete solution. Example 2: Solve . Solution: Taking last two fractions, we have + , on integration, we get or = zc ...... (i) Now using x, y, z as multipliers, we get or On taking fractions , on integration, gives 0r ii (i) and (ii) form the complete solution. Example3: Solve . Solution: On choosing 1, -1, 0 as multipliers, we get or ... (i) Similarly, taking 0, 1, -1 and – 1, 0, 1 as multipliers, we get or …(ii) or … (iii) and From (i),(ii),and (iii) we have 22 Taking first two fractions of (iv), we get on integration gives i.e. ….(v) =c Similarly on taking last two fractions, we get ….. (vi ) = c’ (v) and (vi) form the complete solution . CHECK YOUR PROGRESS: Solve the following simultaneous differential equations: Q.1 [Ans. Q.2 [Ans. Q. 3 [Ans. Q.4 ]. ]. ]. [Ans. Q.5 √ [Ans. ]. ]. 1.8 Total differential equations: An equation of the form P ….. (1) where P, Q, R are the functions of x, y, z is called total differential equation. 23 Sometimes equation (1) can be directly integrable, if there exists a function S(x, y, z) whose total differential dS is equal to L.H.S. (1), otherwise we find the condition for which equation (1) be integrable. 1.8.1 Condition for integrability of Total differential equations P …. (1) The equation is P Let S(x, y, z) = c be a solution of equation (1). Then total differential dS must be equal to P …. (2) , but, we know that On comparing (1) and (2), we have = , say …. (3) so that now from first two equations, we have =P Since and , we have P Gives )=Q =Q =Q …. (4) P Similarly, )=R Q …. (5) and )=P R …. (6) Multiplying (4), (5) and (6) by R, P and Q respectively and adding, we get )+ )+ ) = 0. ….. (7) This is the required condition for integrability of equation (1). 24 Note: (I) If an equation satisfies condition (7), has an integral. In this case given equation is directly integrable. (II) Equation (1) is exact, if obtained by putting these conditions are in (4), (5) and (6). 1.9 Algorithm to solve Total differential equations by different methods with solved examples Let The equation is P …… (1) . Step I: First compare the given equation to P and find P, Q, and R. Step II: Check the condition of integrability )+ )+ ) = 0. .….. (2) Case I: If the given equation is exact then the coefficients of P, Q and R in (2) are zero. In this case given equation is directly integrable after regrouping the terms of the equation. Example 1: Solve (x - y)dx – x dy + z dz = 0 Solution: A.T.Q. , we have P = (x - y) , Q = – x, R = z. The condition of integrability is )+ )+ )–x )=0 )+ )=0 satisfied here. Therefore equation is exact, so on regrouping the terms of the given equation ,we have x dx –( y dx + x dy) + z dz = 0, on integration, we get , is the solution. Example 2: Solve (y + z) dx + dy + dz = 0 Solution: A.T.Q. , we have P = (y + z) , Q = 1, R = 1. 25 Here , = 0, , The condition of integrability is )+ )+ A.T.Q., we have )=0 ) +1 )+ ) = 0, satisfied here. Therefore equation is exact, so on rearranging the terms of the given equation, we have =0, integrating, we get Or = c , is the solution. Case II: If P, Q, R are homogeneous functions then substitute to separate one variable z (say) from the other variables. After separating variable z, equation is directly integrable. Satisfaction of condition of integrability is must in each problem. Students are advised to try the same in each question, and then go ahead as follows: Example 1: Solve (x - y) dx – x dy + z dz = 0 Solution: Here P = (x - y) , Q = – x, R = z. Clearly P, Q, R are homogeneous function, so on putting , we get and the equation becomes ) – uz Or Or Or +{ u- , integration gives, — — Example 2: Solve ( Solution: Here P = ( )+ =0 )} =0 — )+ 2 ) = c, solution is obtained by putting u = + yz ) dx +( xz + + yz ), Q = ( xz + dy + ,R= . - xy) dz = 0 - xy). The equation satisfies integrability condition. Also P, Q, R are homogeneous function, so on putting 26 , we get and the equation becomes ( ( +v +v )( ( )+( )+ + - + Or - ) dz = 0 ) dz = 0 , Or - , integration gives, Or , required solution is obtain by putting u = i.e. complete solution is or (x + z) y = c (y + z). Case III: Check first, which two terms can be readily solved then take remaining third term variable constant viz for our convenience let two terms P integrable, then we take the third variable z = constant , so we get readily . Put in the given equation and then integrate the remaining and take (z) as integration constant. Differentiate obtained relation with respect to x, y, z and then compare with the given equation to find (z), which on integration gives the value of (z) to complete the solution. NOTE: Satisfaction of condition of integrability is must in each problem. Students are advised to try the same in each question, and then go ahead as follows: Example 1: Solve (x - y) dx – x dy + z dz = 0 Solution: Rearranging the equation, we have x dx –( x dy+ y dx) + z dz = 0 …… (1) 27 now, if we put z=constant, then we have dz = 0 and the equation becomes x dx –( x dy+ y dx) = 0, integrating gives, on differentiating, we get (x - y) dx – x dy = comparing (1) and (3), we get ……. (2) …….. (3) dz gives +c putting this value in (2),the solution is , +c [ or Example 2: Solve ( Ans. ) dx + dy + 2 dz = 0 Solution: Taking x= constant, we have dx = 0 and the equation becomes dy + 2 dz = 0, integrating give, ….. (1) now on differentiating (1), we get dy + 2 dz or dy + 2 dz ….. (2) comparing it to given equation, we get = Or I.F. = , =-( = the general solution is ∫ [ by (1)] ), which is linear equation in ∫ . = = ∫( = ∫ ) , therefore ∫ = , putting this value of or . Ans. 28 Note: If none of the above three cases are applicable, then we apply the following case. Case IV: Compare (1) and (2) to obtain auxiliary equation (A.E.) in the form of simultaneous equations as below: . Solve the above equations by the methods discussed in (1.7.3) and (1.7.4) earlier. If two integrals are say, then by comparing with (1), we get the values of A and B and then the complete solution. Example 1: Solve ( + yz ) dx +( xz + Solution: here P = ( dy + + yz ), Q = ( xz + ,R= - xy) dz = 0. - xy) , Thus A.E. is or ….. (1) or Taking last two fractions, we get Also taking 1, 0, 1 as multipliers, we get Or or = v (say ). …. (2) Then A du + B dv =0 gives A )+B + ) =0 [ i.e. Comparing with the given equation, we get ( + yz ) i.e. ( +z) =u 29 And ( xz + Or i.e. = Hence (2) becomes du + u dv =0 i.e. on integration, we get Or Ans. CHECK YOUR PROGRESS: Solve the following total differential equation: Q.1 Q.2 ( +2x) -4 + + = 0. + Q.3 [Ans. + = 0. [Ans. . ] [ Q.4 Q.5 ] . . [ [Ans. ] ] ]. 1.10 Unit brief review (Summary): 30 UNIT –II: Picard’s method of integration, successive approximation, Existence theorem Uniqueness theorem, Existence and uniqueness theorem (All proof by Picard’s method) Unit Structure: 2.0 Object 2.1 Introduction 2.2 Picard’s method of integration 2.3 Algorithm to solve initial value problem by Picard’s method with solved examples 2.4 Existence theorem and its proof 2.5 Uniqueness theorem with proof 2.6 Existence and Uniqueness theorem with proof. 2.7 Unit brief review (Summary) 2.0 Object: This unit is presented in such an easy manner with solved examples that after going through the whole unit students is easily solved the problems related to I.V.P. 31 2.1 Introduction: A differential equation along with the subsidiary conditions on dependent variable and its derivatives at the some value of the independent variable is known as initial value problem (I.V.P.) while subsidiary conditions on dependent variable and its derivatives given at more than one value of the independent variable is called Boundary value problems (B.V.P.). In many of the engineering field and practical problems we are often confronted with the differential equation whose exact solution cannot be found by standard methods. Such problems cannot solve exactly and we find an approximate solution only. This process of finding approximate solution is called iteration method or successive approximation method. 2.2 Picard’s method of integration: Given differential equation ……... (1) ……... (2) = From (1), we have ∫ Or ∫ , ∫ [using (2)] ……... (3) , Now the integral on the right hand side can only evaluated if we known the value of y in terms of x. This is not known. As an initial approximation, we replace y by in R.H.S. of (3) and call the value of y on L.H.S. as first approximation of y, so that ∫ , In the same manner if we substitute y by in R.H.S. of (3) and call the value of y on L.H.S. as second approximation of y, we get ∫ , 32 Continuing in this way, we get better approximation in each step than the preceding one. At the nth step we get ∫ , This process is called Picard’s method of successive approximation. 2.3 Algorithm to solve initial value problem by Picard’s method with solved examples: Algorithm to solve initial value problem by Picard’s method: Step I: Find and by comparing given problem to equation (1) and (2). by putting n=1, 2, 3….. in Step II: Find ∫ and corresponding values given of and Solved Examples: Example 1: Apply Picard’s method up to third approximation to solve . Solution: Here and ∫ Picard’s method is given by So, first approximation ∫ A.T.Q. . We know that the nth approximation by is obtained by putting n=1 ∫ or ∫ [since ] or Similarly, second approximation is obtained by putting n=2 ∫ And, third approximation, ∫ ∫ or . Ans. 33 Example 2: Apply Picard’s method up to third approximation to solve . Solution: Here . The nth approximation by Picard’s and ∫ method is given by So, first approximation ∫ ………. (A) is obtained by putting n=1 in (A), therefore we have ∫ or ∫ [since ] Or Now the second approximation given by ∫ ∫ , - ∫ , + The third approximation given by ∫ Or - ∫ + + . Ans. Example 3: Apply Picard’s method up to third approximation to solve . Solution: Here method is given by and ∫ . The nth approximation by Picard’s ………. (A) 34 So, first approximation is obtained by putting n=1 in (A), therefore we have ∫ ∫ or ∫ [Since ] Or Now the second approximation given by ∫ ∫ ∫ + The third approximation given by ∫ ∫ ∫ Or + + + + . Ans. Example 4: Apply Picard’s method to approximate y up to second approximation and corresponding to for that particular solution of . Solution: Here and . The nth approximation by Picard’s method is given by ∫ and ∫ … (A) 35 So, first approximation is obtained by putting n=1 in (A), therefore we have ∫ or ∫ ∫ [Since ] Or …. and ∫ or (i) ∫ …. (ii) The second approximation given by ∫ ∫ ∫ Or ∫ ∫ Or …. (iii) and or ∫ ( ∫ ) or .... (iv) On putting x= 0.1in (i), (ii), (iii) and (iv), we get particular solution Ans. CHECK YOUR PROGRESS: Solve by Picard’s method up to three approximations (or iteration): 36 Q.1 . [Ans. Q.2 . ]. [Ans. Q.3 ]. . [Ans. ]. Q.4 . [Ans. ]. Q.5 . [Ans. ]. 2.4 Existence theorem and its proof: Existence theorem: The I.V.P. , ……... (1) = has at least one solution y(x), provided the function f(x, y) is continuous and bounded for all values of x in a domain D and there exist positive constants M and K such that | | ………. (2) and satisfies the Lipschitz condition | | | | ……….. (3) for all points in the domain D. Proof: As per Picard’s method, the iterative sequence is given by ∫ , ……….. (4) 37 The initial value problem has a solution if necessarily the sequence {yn(x)} converges to y(x) is either a solution of (1) or of an equivalent integral solution ∫ ……….. (5) , ensure the existence of the limiting function ………… (6) Let be the sum of successive difference ∑ Clearly sequence { …………. (7) } converges if the infinite series ∑ ∫ From (4) converges. , ∫ , = ∫ [ ………. (8) The relation (8) is true for k = 1, 2, 3, ………. , ∫ Also from (4), we have ∫ ……… (9) Now | | ∫ | Or | | | | | ∫ |, from (2) | ……….. (10) = ∫ [ Again from (8) | || | ∫ | ∫ || | || | = | | ∫ | | || from (10) |, from (3) ………(11) 38 Similarly | | | | | On continuing, we have | = | | | | ……… (12) By Mathematical induction we easily establish that (12) holds true for (k+1) i.e. | | | | ....…… (13) From (13) it is quite clear that the value of each term of the series (7) is getting smaller term by term and we get | Or | | [ | | | | | ], which converges for all values ( ) and so, . Taking limit , we get from (4) ∫ , ∫ = ∫ Thus the iterative sequence (4) converges to a solution, hence the theorem. 2.5 Uniqueness theorem with proof: Uniqueness theorem: The I.V.P. , = has a unique solution y(x), provided the function f(x, y) is continuous and bounded for all values of x in a domain D and there exist positive constants M and K such that | | ………. (1) and satisfies the Lipschitz condition 39 | | | | ……….. (2) for all points in the domain D. Proof: Let if possible the I.V.P. , ∫ Then = have two solutions y(x) and v(x). ∫ and ∫ [ | | ∫ | ( )| ∫ )| | | ( | |, from (1) = 2M| …….. (3) Again using (2), we have | | ∫ | || | || |= ……..(4) On combining (3) and (4), we get | | | ∫ | | …… (5) In view of (5) and (6), we have | |= | | Continuing in this way, we get | |= | | , n= 1, 2, 3,…… …..... (6) Right hand side of (6) tends to zero as n tends to infinity, thus | | = 0 gives, , solution is unique. 2.6 Existence and Uniqueness theorem with proof: 40 Existence and Uniqueness theorem: Let function f(x, y) is continuous and bounded for all values of x in a domain D and there exist positive constant M such that | | ………. (1) . Let the function f(x, y) satisfy the Lipschitz condition | | | | Where the constant K is independent of | lie in D, where Mh | | ……….. (2) . Let the rectangle R defined by | ……… (3) . Then the I.V.P. , = has a unique solution y(x). Proof: By Picard’s successive approximation, we have ∫ ∫ … … … … … … … …….. (4) ∫ ∫ First we claim that the function or | ∫ Or | | lies in R. We have | | | ∫ | Mh || | ∫ | |, from (1) by (3). This proves the desired result for n= 1.By Mathematical induction we assume that lies in R. Then from (4), we have | | ∫ | || | ∫ | |, from (1) 41 Or | | Shows that | | Since, we have | | by (3). | | lies in R. Secondly we claim that | | Mh | | Similarly | | | | On continuing, we have ……. (5) | ∫ | ∫ | | || | || | = | | ∫ || | | | = | | | |, from (2) | , conclude that (5) is true for n= 1, 2, 3, ….. . By (3) and (5), we have | | , for n= 1, 2, 3, ….. . ……. (6) Using the infinite series ……. (7) Mh+ [ + + …. ] which is convergent and so series (7) is convergent. Thus . Now we show that tends to +…………+ satisfies the differential , = . Since uniformly in R and by Lipchitz conditions, we have on taking limit , we get from (4) ∫ , 42 ∫ ∫ = …… (8) The integrand on right hand side of (8) being a continuous function, we come to conclusion that the integral has a derivative and thus y(x) satisfies the I.V.P. Hence the solution of given I.V.P. is exist. Uniqueness of solution already proved in (2.5), shall left to the reader. Solved Examples: Example 1. If rectangle R is defined by | | f (x, y)= | | , satisfy the Lipchitz condition. Find the Lipchitz constant. Solution: Here f (x, y)= function on R , show that the function = . Since f is real valued is exist and being continuous and bounded in R. | |=| | As we know that | | | Lipchitz condition. Thus | | | | | and | | | || , where where K= | | | | | +1} is positive constant. Hence f satisfies +1. Example 2: Examine existence and uniqueness theorem for the I.V.P. Solution: Here f (x, y)= and . Clearly f and (x , y). We consider rectangle R as | +1}. | | , = are both continuous for all | about the point (1, -1). Obviously in this rectangle R | | | = (2+2b) | | | || | | Clearly implies that Lipchitz condition is satisfied in R. Now let M= max. | h = min.{a, b/M}, then the problem posesses a unique solution in | | | and . 43 CHECK YOUR PROGRESS: 1. State and prove Picard’s theorem. 2. Show that if the solution of the I.V.P. is exist, is unique. 3. If rectangle R is defined by | | | | show that the function f(x, y)= , satisfies the Lipchitz condition. Find the Lipchitz constant. 2.7 Unit Summary: 44 45 M.Sc. (Previous) Mathematics Paper –V Differential Equations BLOCK- II UNIT –III : Dependence on initial conditions and parameters; Preliminaries, Continuity, Differentiability, Higher order Differentiability. Poincare-Bendixson theory –Autonomous systems, Index of a stationary point, Poincare-Bendixson theorem, Stability of periodic solutions, rotation point, foci, nodes and saddle points. UNIT -IV: Linear second order equations-Preliminaries, Basic facts. Theorems of Sturm, Sturm-Liouville Boundary value problems. Numbers of zeros, Nonoscillatory equations and principal solutions, Nonoscillation theorems. UNIT -V: Partial differential equations of first and second order, linear partial differential equation with constant coefficient. 46 BLOCK-INTRODUCTION Differential equations are mathematically studied from several different perspectives, mostly concerned with their solutions —the set of functions that satisfy the equation. An Ordinary Differential Equation is a differential equation in which the unknown function known as the dependent variable is a function of an independent variable. On the other hand a Partial Differential Equation is a differential equation in which the unknown function is a function of multiple independent variables and the equation involves its partial derivatives. The order of partial differential equations defined similarly to the case of ordinary differential equations. Both ordinary and partial differential equations are broadly classified as linear and nonlinear homogeneous and non homogeneous. This block contains three units: Unit III: This unit deals the concept of Dependence on initial conditions and parameters, Continuity, Differentiability, Higher order Differentiability. Poincare-Bendixson theorem, Index of a stationary point, Poincare-Bendixson theorem, Stability of periodic solutions, rotation point, foci, nodes and saddle points with several solved examples followed by exercise to check the progress of reader. Unit IV: This unit deals the concept of linear second order equations-Preliminaries, Basic facts, Theorems of Sturm, Sturm-Liouville Boundary value problems, Numbers of zeros, Nonoscillatory equations and principal solutions, Nonoscillation theorems with several solved examples followed by exercise to evaluate reader by himself. Unit V: This unit deals the solutions of Partial differential equations of first and second order, linear partial differential equation with constant coefficient by usual and shortcut methods followed by several solved examples and exercise to check the progress of reader. At the end a list of reference books are given for the convenience to the reader. 47 UNIT –III: Dependence on initial conditions and parameters; Preliminaries, Continuity, Differentiability, Higher order Differentiability. Poincare-Bendixson theory –Autonomous systems, Index of a stationary point, Poincare-Bendixson theorem, Stability of periodic solutions, rotation point, foci, nodes and saddle points. Unit Structure : 3.0 Object 3.1 Dependence on initial conditions and parameters; Basic facts 3.2 Continuity, differentiability theorems 3.3Poincare- Bendixson theory-Autonomous systems 3.4 Stability of periodic solutions, rotation point, foci, nodes and saddle points 3.5 Unit brief review (Summary). 3.0 Object: The main object of the present unit is to provide all the important contents in easy manner with several solved examples. 3.1 Dependence on initial conditions and parameters; Basic facts: 3.1.1 Dependence on initial conditions and parameters: If is defined on an open interval with the property that then initial value problem (I.V.P.) ….. (1) has a unique solution 48 defined on maximal interval of existence conditions where depends on initial . In general the I.V.P. depends upon a set of parameters. If , then I.V.P. …. (2) Where for each fixed z, then I.V.P. (2) has a unique solution Now condition (1) can be reduced to the system (2) by change of variables say and . Clearly I.V.P. (1) changes into the form …. (3) ( In which ) considered as parameters. (3.1.2) Lower and Upper Semi continuity: The lower semi-continuity of function ( ) at ( ( ) as ( Similarly upper semi-continuity of function ( ) is defined as ) at ( ( ) ( ) ) is defined as ) as ( ) ( ) 3.2 Continuity and differentiability theorems: 3.2.1 Continuity theorem: Let function property that for every ( be continuous on an open set E with the ) E then I.V.P. …..(1) Let be the maximal interval of existence of the solution , then 49 or function of and is a lower(upper) semi continuous ,( is continuous on the set ) E Proof: The I.V.P. (1) can be replaced by , , i.e. by the I.V.P. and ( Without loss of generality assume that function defined on an open interval ….. (2) ) by does not depend on z. Thus with the property that is Then I.V.P. (2) has a unique solution on the maximal interval of existence ( ), where Now we choose a sequence of points ( )( ( And ) ( ) ) ( ) where as . Since the solution of (2) is unique, we can apply the theorem if and be a sequence of continuous function defined on open set E such that as holds uniformly on each compact subset at E. Let ( And let ( of ) ) be the maximal interval of existence. Also ( ) 50 , then there exist a solution and a sequence of positive integers with the property that then and as uniformly for . In particular lim lub as n = n(k) . Applying the result of this theorem we get Implies is lower semi continuous. Similarly , thus Now we claim that is upper semi continuous. is continuous on ( In the above theorem we have ). as , therefore is considered for a fixed x. We know that is continuous for fixed x, so in the above reference is uniformly continuous i.e. for a given depending upon s.t. | If Hence | | | , then | is continuous on set ( | ). 51 3.2.2 Differentiability Theorem: Let f(x, y, z) be a continuous on an open (x, y, z) set E and possess continuous first order partial derivatives with respect to the components y and z. Then (i) The unique solution Is the class on its open domain ( (ii) Furthermore if J(x)= J y at …. (1) of ) where . is the Jacobian matrix …. (2) J(x)= J Then …. (3) is the solution of the I.V.P. Where and And { is the solution of Where = at of f(x, y, z) with respect to ... . (4) is the vector. and is given by ∑ …. (5) Proof of this theorem is beyond our scope. Reader may find the proof of the above theorem in various reference books given at the end of this block. 3.3 Poincare- Bendixson theory-Autonomous systems: 3.3.1Basic concepts (I) Autonomous system: A system of two first order equations of the form …. (1) 52 is said to be autonomous, when the independent variable t does not appear explicitly. The system (1) defines as …. (2) which is indeterminate only at points where vanish simultaneously. Such points of course represent singular solution of (1). (II) Phase Plane : Let and be continuously differentiable functions in some region R in the xy-plane. Then xy-plane is called Phase plane for the system …. (1) (III) Trajectory or The Orbit: A unique solution of the system …. (1) Is called trajectories or the orbit if it is exists in the phase plane for some open interval and any point ( ) of region R and satisfies the initial condition (IV) Critical Point: A critical point of the system …. (1) 53 is a point ( ) obtained by solving ( ) ( On the other hand if ( and and such that ) ) is a critical point of the system, then the constant valued function satisfying system (1) and are called equilibrium solution of the system (1). A point which is not critical is called regular point. (V) Index: An index is a significant property of a critical point P of the system …. (1) This is an integer which represents the net rotation of the direction field along a simple closed curve. (VI) Centers: A center (or vertex) is a critical point which is surrounded by a family of closed paths. It is not approached by any path as . (VII) Spirals : A spiral (or sometimes called a focus) is a critical point which is approached in a spiral like manner by a family of paths that wind around it an infinite number of times as (VIII) Almost Linear System: The non linear system is of the form ….(1) where are constants. The matrix form of above system is written as ( ) ( )( ) By ignoring the non linear terms ( ) …. (2) in (2), related linear system is 54 …. (3) If we assume that | | , for system (2) then clearly (0,0) is the critical point to related linear system (3). are continuous and have continuous partial derivatives for all (x,y). √ and √ Then (0, 0) is said to be simple critical point of the system (2) and the system is called almost linear. 3.4 Stability of periodic solutions, rotation point, foci, nodes and saddle points: 3.4.1 Stable and Unstable critical Point: A critical point ( that if the initial point ( remains close to ( ) of the almost linear system (2) is said to be stable provided ) is sufficiently close to ( ) t 0. In other words for given | Then critical point ( ) then the point | |( if ) | ) is called stable. Further a point which is not stable is called unstable. The critical point ( ) is called asymptotically stable if it is stable and every trajectory that begins sufficiently close to ( ) also approaches to ( ) as 55 | | . 3.4.2 Nature and Stability of the critical point (0,0) : The linear autonomous system is …. (1) If we assume that | | , then clearly (0,0) is the critical point to related linear system (1). The above system can be written in matrix form as ( ) ( )( ) Then the eigen values of the coefficient matrix equation | Or ( ) are the roots of the eigen | . Let the roots (eigen values) of the above eigen equation are the critical point (0,0) is depend on the values of . Then the nature of and describes as below: 56 Nature of roots Linear system Type Node Node Stability Unstable Almost linear system Type Node Asymptotically Node stable Saddle Unstable point Node Stability Unstable Asymptotically stable Saddle Unstable point Unstable Node or Unstable Spiral Node Spiral Spiral Asymptotically Node or Asymptotically stable Spiral Unstable Spiral Unstable Spiral Asymptotically Asymptotically stable Centre Stable stable stable Centre Indeterminate or Spiral Solved Exampes based on Nature and Stability of the critical point (0,0) : Example 1: For the system of the equations Verify that (0,0) is a critical point . Show that the system is almost linear and discuss the type and stability of the critical point (0,0). Solution: The given system can be written as 57 Where and . Now for critical point we must have and Solving the equations, we get . Thus (0,0) is a critical point. Also = √ and = √ =0 √ √ =0 Therefore the system is almost linear. Again the related linear system for the given system is , whose matrix form is, ( ) ( )( ) Now the eigen values are the roots of the eigen equation | | or | | or Therefore {Here } So from the table (3.4.2) it is clear that the critical point is a spiral and asymptotically stable. 58 Example 2: For the system of the equations Verify that (0,0) is a critical point . Discuss the type and stability of the critical point (0,0). Solution: For critical point we must have and Solving the equations, we get . Thus (0,0) is a critical point. Since ine given system , therefore we have = √ and = √ =0 √ =0 √ and the system is almost linear. Again the related linear system for the given system is , whose matrix form is, ( ) ( )( ) Now the eigen values are the roots of the eigen equation | | or | | or Therefore , here have same sign. So from the table (3.4.2) it is clear that the critical point is a node and unstable. 59 CHECK YOUR PROGRESS: Q.1 For the system of the equations (i) Show that (0,0) and(1,1) a critical point for the given system. (ii) Show that (0,0) is a saddle point and (1,1) is the centre of above system. Q.2 Determine the nature and stability of the critical point (0,0) for each of the following systems of the equations (i) , (ii) , (iii) , (iv) , (v) , . 3.4.3 Periodic Solutions: Let …. (1) is a non linear autonomous system. Then a solution periodic if neither function is constant, if both are defined for all is said to be and if such that 60 Then the number T with the above property is called periodic solution with period T. 3.4.4 Closed Paths: A non linear system of equation has a closed path if it has a periodic solution. On the other hand a linear system has closed path if the roots of its auxiliary equation are purely imaginary. We can easily understand the concept of closed path with the following example: …. (2) By using polar co-ordinates and taking we have and which gives x + and x .... (3) By (2) and (3) we get or .... (4) and .... (5) The system (2) has a single critical point at To find the path consider Separate the variables and integration we get √ …..(6) and The corresponding general solution of (2) is √ and √ … . (7) Geometrically if we analyze (6) we come to the conclusion that for c = 0 , we have 61 as If c 0 then . Also if c and 0 then r 1 and again . The above observations clearly shows that there is only a single closed circular path for . 3.4.5 Poincare- Bendixson Theorems: Theorem 1: A closed path of the system necessarily surrounds at least one critical point of this system. Proof: Let C be a simple closed curve and assume that C does not pass through any critical point of the given system. If P= (x, y ) is a point on C then V(x, y) = is a non zero vector and have a definite direction given by then change by 2 If P moves once around C , n is an integer and it is called the index of C. If C shrinks continuously to a smaller simple closed curve without passing over any critical point, then index change continuously but index is an integer which cannot change. Theorem 2: Let R be a bounded region of the phase plane together with its boundary, and assume that R does not contain any critical points of the system If C = [ x(t), y(t)] is a path of (1) that lies in R for some T and remains in R for all t then C is either itself a closed path or it spirals toward a closed path as , . Thus in either case system has a closed path in R. The proof of this theorem can be easily illustrated by an example we discuss in 3.5.4. 62 UNIT –IV: Linear second order equations-Preliminaries, Basic facts. Theorems of Sturm, Sturm-Liouville Boundary value problems. Numbers of zeros, Nonoscillatory equations and principal solutions, Nonoscillation theorems. Unit Structure : 4.0 Object 4.1 Linear second order equations-Preliminaries, Basic facts 4.2 Algorithm to solve linear second order differential equations with solved examples 4.3 Sturm-Liouville Boundary value problems 4.4 Algorithm to solve Sturm-Liouville Boundary value problems and solved examples 4.5 Number of zeros 4.6 Nonoscillatory equations and principal solutions. 4.7 Nonoscillation theorems 4.8 Unit brief review (Summary) 4.1 Object: 4.1 Linear second order equations-Preliminaries, Basic facts: 4.1.1 Linear second order equations: An equation of type , 63 Where P, Q, X are functions of x alone, is called the linear equation of second order. If the coefficients P and Q are constants, the equation can be solved by the methods which we have discussed in unit II, block I, otherwise there is no general method to solve such equations. We are giving a procedure in which integral belongs to the complementary function can be found. 4.1.2 Complete solution of linear equation of second order, when one integral belonging to the C.F. is known: , The given equation is Let y = ….. (1) be a known part of the C.F. i.e. a solution of equation , Let y =u Now y =u ….. (2) be the complete solution of equation (1), where u is the function of x. gives = and Substituting these values in (1), we get ( ) ( ) Or ( ) ( ) , , ….. (3) = Putting , in (3), we get - = Which is linear in t, therefore I.F. = ….. (4) ∫{ } ∫ Hence solution of (4) is 64 t (I.F.) = ∫ or Solving and putting back u=( ∫ )=∫ , and integrating we get ∫ ∫ ) [∫ Complete solution of (1) is y =u 4.1.3 Special integral part ∫ t( ] , where u is given by (5). : The integral part given below of C.F. to is obtain by putting and Condition .. in (2). Integral part is 1+P+Q=0 1-P+Q=0 + mP + Q = 0 P + Qx = 0 2+2Px + Q X +=0 m(m-1) + mPx + Q +=0 4.2 Algorithm to solve linear second order differential equations with solved examples: 4.2.1 Algorithm to solve linear second order differential equations: Step I: Put the equation in standard form i.e. Step II: Test for integral part Step III: Put y = u coefficient of is unity. of C.F. as per table given in 4.1.3 above. in step I equation and simplify as simplified in 4.1.2. The following method fully illustrates the method. Solved Examples: Example 1: Solve 65 Solution: Given differential equation can be written as Here P = ,Q= and R = 0; Clearly 1 + P + Q = 1 =0 + is a part of C.F. Now we assume that the Complete solution is y =u = 0, Putting , we get x = , then the given equation reduces to or = Integrating both sides, we get log t + log x = log c or t . x = c Putting back , we get . x = c or Integrating both sides, we get = c log + Hence complete solution is y =(c log + = . ) . Ans. Example 2: Solve Solution: Here P = ,Q= and R = x; =0 Clearly P + Q x = is a part of C.F. Assume that the Complete solution is y =u = x, Putting or ( ) = , we get x = ….. (I) , a linear equation in t whose I.F. = Solution of (I) is then the given equation reduces to ∫( ) ∫ 66 . ( ) / =. ( ) / ) [Putting / . This gives ( [ ) = / dx + . Hence complete solution is y = ] / . Integrating both sides, we get ] / or . Or ( . ∫ . / dx + . Ans. CHECK YOUR PROGRESS: Q.1 Solve [Ans. Q.2 Solve [Ans. dx + ]. + ]. Q.3 Solve [Ans. 4.2.2 Reduction of equation + ]. to its normal form: When we fail to obtain a part of C.F. we cannot apply method 4.2.1. In such cases the equation may be solved by reducing equation ….. (1) to its normal form i.e. by removing first derivative as below: 67 Let y =u gives = and Substituting these values in (1), we get ( ) ( Now we choose ( ) ) , ( ) in such a way that the coefficient of , i.e. ( ∫ Or ( Or ) = And , differentiating both sides gives = Putting .. (2) ….. (3) ∫ Again from (2) we have = (3), we get ( Or ) – ( – Or )= ) , where . ∫ ∫ ( – ) = . ∫ or ….. (4) Equation (4) is the required normal form of equation (1), can be easily integrated. 4.2.3 Algorithm to solve linear second order differential equations by reducing to its Normal form: Step I: Put the equation in standard form i.e. coefficient of is unity. 68 ∫ Step II: Find P, Q, X, ( – , Step III: Put the values of ) in normal form = . ∫ . . Step IV: Obtained equation is linear differential equation with constant coefficient and solve by finding C.F. and P.I. as given in Unit I, Block I. Step V: Required Solution is obtain by putting the values of Solved Examples: Example 1: Solve Solution: Here P = Q= ∫ = , ( – And to reduce in normal form we choose ) = ( – )= 2 Equation reduces to . C.F. = Here A.E. is And P.I. = = ( , therefore ) , And the complete solution is Ans. Example 2: Solve Solution: Here P = Q= To reduce in normal form we choose the complete solution of the given equation as , where ∫ = = , = 69 ( – And ) ( – )= 6 Equation reduces to . C.F. = Here A.E. is And P.I. = = ( ) ( , ) Therefore , And the complete solution is Ans. CHECK YOUR PROGRESS: Q.1 Solve [Ans. Q.2 Solve ]. [Ans. Q.3 Solve ]. [Ans. Q.4 Solve ]. [Ans. 4.2.4 Transformation of the equation ]. by changing the independent variable: Sometimes the equation is transformed into an integrable form by changing the independent variable. Let the equationbe ….. (1) Let the independent variable x be changed to z by taking z as the function of x. and ( ) ( ) , 70 Substituting these values in (1), we get [ ( ) , ( ) Or [ ] + , …. (2) Or Where ( , ) ( and ) ( ) . After obtaining equation (2) we like to choose z in such a way that (2) can be easily integrated. Case I: . We choose z to make the coefficient of ( or ) integrating again, we get or Integrating, we get log or which can be easily solved provided by in (2), equal to zero i.e. , this value of z reduces (2) to comes out to be a constant or a constant multiplied . Case II: We choose z such that i.e. ( ) or ( ) (constant), , integrating gives ∫ . 71 The above value of z reduces (2) to Which can be solved easily, if comes out to be a constant. Things To Remember: To apply these methods student are advised to remember equation (2) and the values of .Since we have assumed that z is f(x),so to get the required the required result replace z with its corresponding value in terms of x. Solved Examples: Example 1: Solve . Solution: Writing given equation in standard form, we have …. (i) Here P = ,Q= , X= Changing independent variable from x to z, equation becomes ….. (ii) Where, ( ) , ( Let us choose z such that i.e. Then ( ) ( or ( ) and ) ) ( ) . = - 2 (constant) , integrating gives and ( . ) ….. (iii) = 72 Hence equation (ii) transferred to Or Or By (iii) ….. (iv) Or Now C.F. = + and P.I. = = , Hence the solution of the given equation is y= + + Or y = + + Ans. [ As z = sin x ] Example 2: Solve Solution: Writing given equation in standard form, we have …. (i) Here P = ,Q= , X= Changing independent variable from x to z, equation becomes ….. (ii) Where, ( , ) ( ) and ( ) . Let us choose z such that A.T.Q. ( ) , putting then ( ) Separating the variables and integrating, we get 73 + Since or gives ( ) ( = ) ( Separating the variables and integrating, z = = 4 and ( ( ) ) , ) ( ( ) ) Hence the transformed equation is its solution is y = or + y= + [As z = ] Ans. CHECK YOUR PROGRESS: Q.1 Solve [Ans. y = Q.2 Solve ]. [Ans. y = Q.3 Solve Q.4 Solve + [Ans. y = [Ans. y = + ( ) + - ]. ]. ]. 4.3 Sturm-Liouville Boundary value problems: 4.3.1 Boundary value problems: A differential equation along with some conditions on the unknown function and its derivatives given for one specific value of the independent variable is called initial value 74 problem (I.V.P.). If conditions on the unknown function and its derivatives given for more than one value of the independent variable is called boundary value problem (B.V.P.). Ex. (1) y is an I.V.P. because for one value independent variable we have conditions on the unknown function derivatives of and its . Ex. (2) y is a B.V.P. because for two values of independent variable we have conditions on the unknown function its derivatives and . 4.3.2 Sturm-Liouville Boundary value problem: The boundary value problem given by the second order differential equation of the form [ [ ….. (1) on some interval [a , b] satisfying the conditions of the form (i) and (ii) , where are a real parameter, the real valued continuous functions of x . Also and are real constants at least one in each is non zero of conditions (2) is called the Sturm Liouville Problem. Here equation (1) is known as Sturm Liouville equation and equation (2) is called boundary conditions. Note (1): Two real functions f(x) and (x) are called orthogonal function on the interval [a , b] , if ∫ (2) If are two eigen functions of Sturm Liouville Problem (1) corresponding to two distinct eigen values respectively and their derivatives are continuous function in the same interval [a, b], then orthogonal . 4.3.3 Theorem: Eigen values of the Sturm Liouville Problem are all real. Proof: Sturm Liouville Problem is [ [ …….. (i) On the interval [a, b] satisfying the conditions 75 and We assume that p(x) 0 when eigen value . Let y(x) is an eigen function corresponding to an , then this function satisfies equations (i), (ii), and (iii), and it may be a complex function. Taking conjugate throughout in equations (i), (ii), and (iii), we get [ ] [ ] ….(iv) .… (v) …. (vi) and The above equations shows that is eigen function corresponding to the eigen value . Multiplying (i) by = [ ) = and (iv) by and subtracting, we get [ ] ] [ =*[ ] [ Integrating both sides from a to b, we get )∫ *[ )∫ Thus we have [ ] +[ + = 0, in view of (ii), (iii), (v) and (vi). )= 0 or as ∫ .Hence the theorem. 4.4 Algorithm to solve Sturm-Liouville Boundary value problems and solved examples: Algorithm to solve Sturm-Liouville Boundary value problems: Standard form of the equation is Step I: Consider three different cases for in each problem. Step II: solve each case separately by finding C.F. 76 Step III: Use boundary condition to find the values of constants taken in C.F. Conclusion: Non zero value of constant gives eigen value and corresponding eigen vectors. Solved Examples: Example 1: Findthe eigen values and eigen functions of the Sturm Liouville problem Solution: Case I. If , the given equation reduces to Integrating twice gives Using given boundary conditions, we have B = 0 and Therefore Case II. If or , which is not an eigen function. ,let the given equation reduces to which is a linear differential equation with constant coefficient, whose solution is Using given boundary conditions, we have A + B = 0 and Therefore Case III. If or , which is not an eigen function. ,let he given equation reduces to which is a linear differential equation with constant coefficient, whose solution is Using given boundary conditions, we have A = 0 and or If B = 0then we have, , which is not an eigen function. Therefore taking , we have 77 Therefore the eigen function are taking B = 1 and the eigen values are ……. Ans. Example 2: Find the eigen values and eigen functions of the Sturm Liouville problem Solution: Case I. If , the given equation reduces to Integrating twice gives and Using given boundary conditions, we have B = 0 and Therefore Case II. If . , which is not an eigen function. ,let the given equation reduces to which is a linear differential equation with constant coefficient, whose solution is and Using given boundary conditions, we have A + B = 0 and Therefore Case III. If or , which is not an eigen function. ,let he given equation reduces to which is a linear differential equation with constant coefficient, whose solution is and Using given boundary conditions, we have A = 0 and or If B = 0then we have, , which is not an eigen function. Therefore taking , gives 78 Therefore the eigen function are taking B = 1 and the eigen values ……. are Ans. CHECK YOUR PROGRESS: Q.1 For the eigen values and eigen functions of the Sturm Liouville problem [Ans. eigen function eigen values are …….] Q.2 Find the eigen values and eigen functions of the Sturm Liouville problem [Ans. eigen function …….] eigen values are Q.3 Find the eigen values and eigen functions of the Sturm Liouville problem [Ans. eigen function eigen values are …….] 4.5 Number of zeros: We shall now discuss the problem of determining the number of zeros of non trivial solutions of the general second order differential equation ….. (1) + q(t)y = 0 Where the functions and q(t) are continuous on some interval Before moving to the main result, we first understand the Prufer’s transformation. 4.5.1. Cor. I: let ….. (i) be a non trivial solution of existing on the interval , then the transformation , reduces (i) to ….. (ii) 79 ….. (iii) ( and ) ….. (iv) Proof: The proof of the above cor. is quite usual and it can be obtained by differentiating (iii) with respect to t and using relation y = 4.5.2 Let the coefficient function P(t) continuous on the interval and . and q(t) in be and let u(t) be a nontrivial solution of , Suppose u(t) has exactly n ( [a, b]. If . is a function defined by (ii) then and . t= for are the zeros of y(t), it follows from the second relation of (ii), at t = From the continuity of ) , k=1 , 2 , 3…… . for Proof: Since t = that ) zeros at t = .Thus from(iii), we have this implies that is increasing in the nbd of the points , gives the result. 4.5.3 The Sturm Theory: Let u(t) be a solution of ….. (A) having first derivative u (t) on zeroes on then Also suppose that u(t) has an infinite number of for all t on Proof: Since u(t) has an infinite number of zeros on theorem the set of zeroes has a limit point l [ then by Bolzano-weirstrass . Therefore there exists a sequence 〈 〉 of zeroes which converges to l. Now since u is continuous , 80 through any sequence of points on [ where Let through the sequence of zeroes 〈 〉. Then . Now since u (t) exists, then by the definition of differentiability, we have u (l) = , through any sequence of zeroes 〈 〉. where For such points Therefore u (l) = 0. is a solution of equation (A) such that therefore [ u (l) = 0. Since l , . This concludes the theorem. 4.5.4. Sturm Separation Theorem: Let u and v be real linearly independent solutions of ….. (A) on t [ . Then between any two consecutive zeroes of u, there is precisely one zero of v. are two consecutive zeroes of u on [ Proof: Let have v( ) v( ) . Then by theorem 4.5.3 we . Now we assume that v has no zero in the open interval . Then since solutions u and v have continuous derivatives on [ has a continuous derivative on the interval [ . Also u/v is zero at the end points of this interval, then by Rolle’s theorem c such that * + , at t = c . But , the quotient u/v * + point and therefore u and v are linearly independent on [a, b]. Thus * + , which is a contradiction. Hence v has a zero in v has precisely one zero in Now to claim that we assume that v has more than one zero in be two consecutive zeroes of v. Then by interchanging the role of 81 u and v and proceeding as above, we can easily claim that u must have at least one zero in This concludes the theorem. 4.6 Nonoscillatory equations and principal solutions: 4.6.1 Oscillatory and Nonoscillatory equations and principal solutions: Equation is of the form where is a real valued and continuous function on is called oscillatory equation if it’s all the non trivial solutions have an infinite number of zeros on These nontrivial solutions are known as oscillatory solutions. If some of them are oscillatory and remaining are non oscillatory, then such equations are known as nonoscillatory equation and their solutions are known as non-oscillatory solutions. Clearly the equation , is non oscillatory equation, if is a real valued and continuous function in the interval . In other word an equation is called non oscillatory in interval , if no solution of the equation can change its sign more than once in the interval. 4.7 Nonoscillation theorems: Theorem: If all the non trivial solutions of continuous, and are oscillatory, is , then some non trivial solutions of (1) are oscillatory. On the other hand, if some non trivial solution of equation (1) are non oscillatory and then some non trivial solutions of must be non oscillatory. 82 Proof: Let and are the non trivial solution of (1). Multiplying (1) by y and the equation by z, and subtracting we get [ [ Or Let and ……..(2) be two consecutive zeros of on the interval [ and let and that y(t) ]. By integrating equation (2) ranging from , we get ∫ [ Since and be two consecutive zeros of and ….(3) , we have = 0 with Therefore from (3), we get ∫ [ Now we claim that has no zero in [ and contradiction. Thus [ ….. (4) has a zero in the interval [ ].Then and ]. On the contrary we assume that does not change its sign. Since and are non negative in the interval [ has a zero in [ ] and ] leads a changes its sign in the interval ]. This concludes that between any two consecutive zeros of there is a Second part follows with the similar argument as we used in first part. 4.8 Unit brief review: 83 UNIT -V: Partial differential equations of first and second order, linear partial differential equation with constant coefficient. Unit Structure : 5.0 Introduction 5.1 Partial differential equations of order one and two with examples 5.2 Linear and non linear Partial differential equations 5.3 Linear partial differential equation with constant coefficient 5.4 Algorithm to solve linear partial differential equations with constant coefficient and solved examples. 5.5 Unit brief review (Summary) 5.0 Introduction: In the theory of partial differential equations, a variable (say) is a function of more than one independent variable. Here we are restrict ourselves to assume that is a function of two variables x and y; we write than derivatives of with respect to derivatives of are denoted by the partial of order one. Similarly the second order partial . A partial differential equation is a relation between dependent variable, independent variables and partial derivatives of dependent variables. 5.1 Partial differential equations of order one and two with examples: 5.1.1 Partial differential equations of order one: The order of partial differential equation is determined by the highest order partial derivative in it. Thus a Partial differential equation is called of order one if it contains highest derivatives of order one i.e. the equation contains Symbolically denoted as (1) (2) only. . For examples, , , 84 are partial differential equations. Clearly order of partial differential equation (1) is one and order of (2) is two. We already studied the formation of partial differential equations by eliminating arbitrary functions and arbitrary constants in earlier classes. Thus in the present chapter we restricted ourselves to discuss the solution of linear partial differential equations with constant coefficients. 5.2 Linear and non linear Partial differential equations: 5.2.1 A partial differential equation involving partial derivatives p and q and no higher derivatives is called of order one . In addition, the degree (or power) of p and q is unity (one), then it is a linear partial differential equation of order one. For example, 2xp+ 5yq= z and p are both linear partial differential equation of order one. On the other hand + 3q = z and x +y are both partial differential equation of order one but non linear. The solutions of linear partial differential equation of order one we already been studied in the previous classes using Lagrange’s and Charpit’s method. So we left this to the reader as practice exercise. 5.3 Linear partial differential equation with constant coefficients: 5.3.1 A partial differential equation is called linear partial differential equation with constant coefficient if it consist higher partial derivatives of z with respect to x and y but the power of each derivatives that occurs and variable z is one and the coefficient of various terms are constant quantities. The general form of such equation is denoted as {[ ] +[ ] + ………. +[ In (1) the notations and coefficients In brief (1) can be written as ] + N}z = taken for the operators ....(1) respectively and the are all constants. , where 85 {[ ]+[ + …+[ ] ] + N}. 5.3.2 Homogeneous linear partial differential equation: Equation (1) in 5.3.1 is called homogeneous linear partial differential equation if i.e. {[ ]} {[ ]}z = in (1) is in the form ….(2) The general solution of homogeneous linear partial differential equation consist two parts namely complementary function (C.F.) and particular integral (P.I.).ion Thus the general solution of above equation is z = C.F. + P.I. Here C.F. is the general solution of of while P.I. is the any particular solution . 5.3.3Solution of homogeneous linear partial differential equation: Complementary function(C.F.): We know that C.F. is the general solution of Let one of the linear factor of ( - ) or The above equation is clearly of Lagrange form is say ( - ).Then we have or p . Therefore the auxiliary equations of Lagrange’s are First two relations give relation gives ( – ) is (on integration). And the last Hence the required solution corresponding to factor In the same fashion solutions corresponding to all factors can be obtained. 86 Particular solution (P.I.): We know that P.I. is the solution of free from arbitrary constants. P.I. of any equation is taken as function can be treated as algebraic function of expanded in ascending powers of or Note: In solving the problem treat and . Here the symbolic and and can be factorized or to obtain the required P.I. as partial operator w.r.t. x and y while as integration w.r.t. x and y respectively. In the next section we will discuss the shortcut method to find P.I. in detail. 5.4 Algorithm to solve linear partial differential equations with constant coefficient and solved examples: 5.4.1. Algorithm to find complementary function (C.F.): Let the given equation is . We follow the steps given below to obtain the required C.F: Step I: Take and put to get auxiliary equation (A.E.). Step II: Find the roots of the equation i.e. values of m. Let the values of m are . Then Case I: If values of m i.e. all are distinct, then C.F. is given by Case II: If values of m are repeated, say . then, C.F. is given by SOLVED EXAMPLES: Example 1: Solve ( Solution: The auxiliary equation is )z=0 {put and =1} 87 Or gives Therefore required C.F. is Or Example 2: Solve ( Ans. )z = 0. Solution: The auxiliary equation is Or {put gives and =1} , Therefore required C.F. is Ans. Example 3: Solve Solution: Given equation is Here A.E. is or gives [repeated roots] Therefore required C.F. is ( ) ( ) Ans. Example 4: Solve ( Solution: A.E. is ( Or ) Or Or )= 0 or [repeated roots] , therefore required C.F. is Ans. 88 5.4.2 Shortcut Methods To Obtain Particular Integral (P.I.):Let the given equation is , then Case I: When is a function of , then to get the particular solution of the , proceed as below [where F equation is a homogeneous of degree n]: function of Algorithm: Step I: Put and integrate Step II: Put a for D and b for , n times with respect to t. to get . [nth integral of Step III: Required P.I. = obtained in step I], where . SOLVED EXAMPLES: Example 1: Solve ( )z= Solution: The auxiliary equation is Or Now =1} P.I. = Here is a function of the form function of degree 2 in We first put and [as is a homogeneous . , and then integrating twice with respect to t, we get ] Also putting 2 for We get and gives Hence C.F. is i.e. {put and 3 for in P.I. = = Therefore the complete solution is . Ans. 89 Example 2: Solve ( )z= Solution: The auxiliary equation is Or {put and =1} gives Hence C.F. is Now, P.I. = = ( =( ) ( ) Here ) ( [ ) is a function of the form homogeneous function of degree 2 in We first put and . , and then integrating i.e. [as Now on putting twice with respect to t, we get ]. for and P.I1. = is a ( in , we get …(A) ) The same procedure we apply for second part i.e. for and putting m for We get P.I2. = ( and for , we get in . ) Therefore the complete solution is P.I1.+ P.I2. ( ) ( ) Ans. 5.4.3 Exceptional case when F(a,b) = 0: If On Putting a for D and b for in , the above method gets fails. Then to evaluate P.I., proceed as below: 90 Step I: Differentiate with respect to D partially and multiply the expression by x, so that Step II: If is also zero, Differentiate with respect to D partially and multiply the expression by x again, so that Repeat the above procedure till the derivative of vanishes. If , then , solved by shortcut method above. Example 1: Solve ( )z=4 Solution: The auxiliary equation is Or {put and =1} gives Hence C.F. is Or Now P.I. = 4 the denominator becomes zero when , 2, therefore differentiating the denominator with respect to D partially and multiplying by x, P.I. = 4 , The denominator again vanishes when 2, therefore again differentiating the denominator with respect to D partially and multiplying again by x, P.I. = 4 The denominator does not vanishes when , 2, therefore 91 P.I. = cos [ = , Therefore the complete solution is . CHECK YOUR PROGRESS: Q.1 Q.2 Solve ( )z= [Ans. [Ans. Q.3 Solve ( )z= [Ans. Q.4 Solve ( )z=( ) [Ans. ]. ]. ]. ]. 5.4.4 General method of finding the P.I. Consider the equation . , which is of Lagrange’s form. Therefore This can be written as A.E. is The first two relations give …(1) . Taking On integration we have z =∫ Thus z =∫ or z =∫ , from (1) , Where constant c is replaced by , after integration. The above method be repeated for all the factors of . 92 Solved Examples: Example 1: . Solution: Given equation can be written as ( )z= The auxiliary equation is {put and =1} Giving Hence C.F. is For finding P.I., we use general method. P.I. = Or = Or = ∫ , [Because corresponding to the factor Or = [ Or = [ Or = [ Or = we have solution . , , replacing , ∫[ [Because corresponding to the factor we have solution . Therefore P.I. = [ Or =[ 93 =[ Or , Thus the complete solution is CHECK YOUR PROGRESS: Q.1 [Ans. Q.2 ( )z= ]. . [Ans. ]. 5.4.5 Non homogeneous linear equations: A linear partial differential which is not homogeneous is called non homogeneous linear equation. We consider the differential , Where is not necessarily homogeneous. We Classify mainly in two types, which shall be treated separately as : (i) is reducible i.e. form can be expressed as product of linear factors of the , where a and b are constants. (ii) is irreducible i.e. Case I: is not reducible as above. is reducible: Complementary function (C.F.) : Let be the factor of Lagrange’s equation concept and write as , then C.F. is easily obtain by applying , where is an arbitrary function. We now give the different cases here: (I) have distinct linear factors: If Where all the factors are distinct, then corresponding C.F. is given by 94 ( ) ( ) Things to remember(1): If the linear factor is is , then the corresponding C.F. . (II) have repeated linear factors: Let a factor occur twice in . Then corresponding C.F. is given by ( ) [ In general if occur n times in , then the corresponding C.F. is given by ( ) [ Things to remember(2): If the linear factor is twice, then the corresponding C.F. is and occurs [ . Particular integral (P.I.): Particular integral of non homogeneous partial differential equation can be found in a way similar to those of ordinary differential equations. Here we are giving some standard form of Form I: If = ’ , then corresponding P.I. in Required P.I. = Form II: If and their corresponding P.I. format. = obtain by putting is obtain by putting .e. , [Provided ]. , then corresponding P.I. and is in Required P.I. = .e. , If the terms of D and D’ remains in the denominator, then 95 (i) Rationalize the denominator and replace and (ii) Treat D and D’ as partial differential operator to get required P.I. Form III: If = , then corresponding P.I. in terms of Form IV: If . = , then corresponding P.I. ’ putting is obtain by expanding in is obtain by then solve according to the form of V Required P.I. = . SOLVED EXAMPLE: Example 1: Solve ( )( ) . Solution: According to 5.4.5 case I (I) of complementary functions, here are two distinct linear factors ( ) and ( ). After comparing to and . Therefore C.F. = we have . [See things to remember(I), 5.4.5] Now P.I. = , writing 2 for D, -1 for = , Hence the complete solution is = + Example 2: Solve ( ) Solution: The given equation is ( )( .. Ans. . )z= , 96 According to 5.4.5 case I (I) of complementary functions, C.F. = Now P.I. = ( . [See things to remember (I), 5.4.5] ) ( , writing ) = for , (-1).2 for , Rationalizing the denominator, we get = , writing . = = [ = Hence the complete solution is = [ Example 3: Solve ( ) Solution: The given equation is ( )( …Ans. . )z= , According to 5.4.5 case I (I) of complementary functions, C.F. = Now P.I. = . [See things to remember (I), 5.4.5] [ = [ 97 = [ = [ = [ = [ . [Taking D and means partial differentiation w.r.t. x and respectively] Hence the complete solution is = [ . …Ans. CHECK YOUR PROGRESS: Q.1 Solve Ans. [ Q.2 Solve ( [ )z = Ans. [ Q.3 Solve ( Case II: ]. )z = xy is irreducible i.e. Ans. [ ]. ]. is not reducible into linear factor: In this case we follow the few steps to find C.F. as below: Step I: Take trial solution = Step II: Put the above value in the given equation. Step III: Find the value of k in terms of h. Put the value of h in trial solution to get the required solution. NOTE: To find Particular Integral we adopt the same technique as in case I. Solved Examples: Example 1: Solve ( )z = 0. 98 ….. (1) Solution: Let the solution of the above equation be = Then from the equation ( ) =0 ( ( Putting ) =0 ) in (1), we get the required solution is … Ans. = Example 2: Solve ( )z= . ….. (1) Solution: Let the solution of the above equation be = Then from the equation ( ) =0 ( ( Putting ) =0 ) in (1), C.F. is = Now P.I.= ( ( = ) , writing ) for , After rationalizing the denominator, we get = 99 = = ( ) ( ) = ( = ( , writing for ) ) Hence the complete solution is = ( ) … Ans. CHECK YOUR PROGRESS: Q.1 Solve ( )z=0 [Ans. = Q.2 Solve ( )z=0 [ = ] ]. 100 LIST OF REFERENCES: 1. Evans, Lawrence C.(1-CA) Partial differential equations. (English. English summary) Graduate Studies in Mathematics, 19. American Mathematical Society, Providence, RI, 1998. 2. Folland, Gerald B.(1-WA) Introduction to partial differential equations. (English. English summary) Second edition. Princeton University Press, Princeton, NJ, 1995 3. Andrews, Larry C. Elementary Partial Differential Equations with Boundary Value Problems New York, NY: Academic Press, 1986. 4. Balachandra Rao S. & H.R. Anuradha Differential Equations 5. Sankara Rao, Introduction to Partial Differential Equations 6. Raisinghania, M D Ordinary And Partial Differential Equations book 7. Raisinghania, M D , Advanced Differential Equations book 8. Zafar Ahsan, Differential Equations and Their Applications, Prentice Hall of India, New Delhi 1999. 9. B. D. Sharma, Differential Equations, K. N. R. N. Delhi 1983 10. Kapoor, N. M. Differential Equations, P.P. Company Ltd. Delhi 1999. 11. Renardy, Michael Rogers, Robert C.(1-VAPI) An introduction to partial differential equations, Springer-Verlag, New York, 1993 101