154x Filetype PDF File size 0.16 MB Source: amsslive.files.wordpress.com
MATH168-DIFFERENTIALEQUATIONSI solutions to first order odes - exact differential equations 1 1 This is the first part of lecture four Lecturer: Dr. Peter Amoako-Yirenkyi - A series of lecture notes primarily taken from the text book (Elementary 2 differential Equation and boundary Recommended Textbook:Elementary differential Equation value problems by Boyce and Diprima. The laws of the universe are written largely in the language of math- After going through this lecture notes, ematics. Algebra is sufficient to solve many static problems, but the you should be able to: most interesting naturally phenomena involve change and are best de- • Youshouldbeable to understand scribed by equations that relate changing quantities. Many important the principle behind Exactness of and significant problems in engineering, the physical sciences, and the differential Equations social sciences such as economics and business when formulated in • Find the solution to exact differen- mathematical terms require the determination of a function satisfying tial equations an equation containing the derivatives of unknown function. Such • Solve non-exact differential equa- equations are called differential equation. tions using method of integrating factors 2 W.W. Boyce & R.C. DiPrima. Elemen- Exact Differential Equation tary Differential Equations and Boundary Value Problems. John Wiley & Sons, Inc., tenth edition, March 2012. ISBN Through out the previous weeks we have been discussing first order 978-0-470-45832-7 differential equations and a number of methods that are applicable to various classes of problems. The most important of these had been separable differential equations and linear coefficients of the differential form of some differential equations(both homogeneous and non-homogeneous). Although these methods can solve a number of differential equations, first order equations that can be solved by these methods are rather special since most first order differential equations cannot be solved in this way. The next class of equations we need to consider are known as exact differential equations which employs the idea of total differential. Given a function f(x,y), its total differential, df , is defined as: df = ∂f dx+ ∂f dy (1) ∂x ∂y This shows that the family of curves (or general solution) f(x,y) = c satisfies the differential equation df =0 (2) =⇒ ∂fdx+∂fdy =0 (3) ∂x ∂y So if there exists a function f(x,y) such that: M(x,y) = ∂f and (4) ∂x N(x,y) = ∂f (5) ∂y 3 3 then M(x,y)dx+N(x,y)dy is called an exact differential , and the In order to establish that a differential equation: equation is exact, recall that in partial M(x,y)dx+N(x,y)dy = 0 (6) derivatives if f(x,y) has continuous second order partial derivatives, then: is said to be an exact equation, whose solution is the family of func- ∂2f ∂2f tions f(x,y) = c . ∂x∂y = ∂y∂x (that is, the order of differentiation does not matter). So if ∂f = M and ∂f = N ∂x ∂y then: ∂M = ∂N (which is ∂2f = ∂2f ) ∂y ∂x ∂x∂y ∂y∂x Not only is this condition necessary for exactness, it also determines exactness. math168 - solutions to odes - exact differential equations and integration factors 2 Example 1. Solve the following differential equation: (1−2xy)dx+(4y3−x2)dy =0 Solution With (M(x,y) = 1−2xy) and N(x,y) = 4y3 −x2, we notice that: ∂M = −2x = ∂N So the given functions M(x,y) and N(x,y) passed ∂y ∂x 4 the test for exactness . All that is left to do is to find the function 4 f (x, y) such that: ∂f = M and ∂f = N. The way to accomplish this ∂x ∂y Example 2. Solve the differential equation is to integrate M with respect to x, integrate N with respect to y, and dy then merge the results: (3x2 +2xy)+(x+y2) =0 dx Solution Z Z Here M = 2x, N = 1,thegiven 2 x x M(x,y)dx = (1−2xy)dx = x−x y+h(y) (10) equation is not exact. To see that it Z Z cannot be solved by the procedure M(x,y)dx = (4y3 −x2)dy = y4 −x2y+g(x) (11) above, let us suppose that there is a function f(x,y) such that: These calculations imply that a function f(x,y) which satisfies both: fx = 3x2 +2xy, fy = x+y2 (7) Integrating the first of (7) gives: ∂f = M =1−2xy and ∂f = N =4y3−x2 is f(x,y) = x2y+y4 f (x, y) = x3 + x2y + ψ(y) (8) ∂x ∂y where ψ(y) is an arbitrary function of y Therefore, the exact equation we were given is satisfied by the family only. To try to satisfy the second of (7) of curves: x − x2y+y4 = c wecompute fy from (8), obtaining: f = x2+dψ =x+y2 or Example 3. Solve the differential equation (ycosx + 2xey) + (sinx + y dy x2ey +2)dy = 0 dψ 2 2 dx dy = x+y −x (9) Solution Since the right-hand side of (9) depends M((x,y) = ycosx+2xey), N(x,y) = sinx+x2ey+2 on x as well as y, it is impossible to ∂M = cosx+2xey, ∂N = cosx+2xey Therefore, the differential solve (9) for ψ(y). Thus there is no ∂y ∂x f (x, y) satisfying both of (7). So in the equation is exact. We know that ∂f = M(x,y); ∂f = N(x,y). Thus next section we will try to find a factor ∂x ∂y to make the equation exact and have a there is an f(x,y) such that function satisfying both (7) - the factor ∂f ∂f is called integrating factor. =ycosx+2xey and =sinx+x2ey+2 (12) ∂x ∂x 5 Integrating the first of these equations: R dy = R (ycosx +2xey)dx, 5 Note that either of equations (10) keeping y constant gives: having h(y) as a function of y alone or (11) with g(x) as a function of x alone f (x, y) = ysinx + x2ey +ψ(y) (13) maybechosenconveniently and hence differentiating the new function w.r.t y keeping x constant: ∂f = sinx+x2eydψ = sinx+x2ey +2 i.e. the differential is equal to ∂y dy R R the second equation in (12). Thus dψ = 2 =⇒ dψ = 2dy and dy ψ = 2y the constant of integration can be omitted; we do not require the most general one. Substituting for ψ(y) in (13) gives f(x,y) = ysinx+x2ey+2y. Before f(x.y) = c, constant. Hence c = ysinx+ x2ey +2y math168 - solutions to odes - exact differential equations and integration factors 3 Using Integrating factors for Non-Exact Equations (I f )M(x,y)dx + N(x,y)dy = 0 (14) is not exact, sometimes we can turn it into 6an exact differential equation by multiplying the whole equation by an appropriate factor, called an integrating factor, let say µ = µ(x,y), which is a function of x, y . For which we obtain: 6 Remark: Finding the integrating factor for a given equation can be very µM(x,y)dx+µN(x,y)dy =0 (15) difficult. Some of the important rules and procedures follow. and expect it to be exact. Then for (15) to be exact then: Theorem1. 1. Assume µ = µ(x) =⇒ ∂µ = µy = 0, then: µx = ∂y ∂(µM) ∂(µN) dµ M −N = = y x,whichisafunctionof ∂y ∂x dx N x alone, call this function ξ(x). Then Applying product rule we get: µ(x) = eR ξ(x)dx is an integrating factor of (14) 2. Assume µ = µ(y) =⇒ ∂µ = ∂M ∂µ ∂N ∂µ ∂x dµ M −N µ +M =µ +N 0, then: µy = = y x,whichis ∂y ∂y ∂x ∂x dy −M a function of y alone, call this function ∂µ ∂µ ψ(y). Then µ(y) = eR ψ(y)dy is an =⇒ M −N =µ(∂N∂x−∂M∂y) (16) integrating factor of (14) ∂y ∂x 3. If ∂M + ∂N 6= 0, then: µ(x,y) = ∂x ∂y 1 , is the integrating factor of (14) Simply put Mµ −Nµ =µ(N −M ), µ(M −N ) ∂M ∂N y x x y y x ∂x + ∂y Sometimes an integrating factor may be found by inspection, after 4. If M(x,y)dx + N(x,y)dy = 0 can grouping the terms in the equation, and recognizing a certain group be written in the form yg(x,y)dx + 7 xh(x,y)dy = 0 where h(x,y) 6= as being a part of an exact differential. The following tabulates g(x,y), then µ(x,y) = 1 possible terms and their respective Integrating factors and the exact xM−yN equations they produce: 7 The following observations are often helpful in finding integrating factors: 1. If a first-order differential equation contains the combination xdy + GroupofTerms Integrating factor Exact differential df(x,y) ydy = 1d(x2 +y2) try some function 2 22 ydx−xdy −1 xdy−ydx = d y x +y asamultiplier. x2 x2 x 2. If a first-order differential equation ydx−xdy 1 ydx−xdy = d x contains the combination xdy + y2 y2 y ydx = d(xy), try some function xy as ydx−xdy −1 xdy−ydx = d ln y a multiplier. xy xy x ydx−xdy − 1 xdy−ydx = d arctan y 3. If a first-order differential equation x2+y2 x2+y2 x 1 ydx+xdy contains the combination xdy−ydx, ydx+xdy xy xy =d(ln(xy)) try 1 or 1 as a multiplier. If h i x2 y2 ydx+xdy 1 , n > 1 ydx+xdy = d −1 neither of these works, try 1 or (xy)n (xy)n h(n−1)(xy)n−1 i xy ydy+xdx 1 or some function of these ydx+xdy 1 =d 1ln(x2+y2) x2+y2 x2+y2 x2+y2 2 expressions, as an integrating factor, ydx+xdy 1 , n > 1 ydy+xdx = dh −1 i remember that dtan−1 y = (x2+y2)n (x2+y2)n 2(n−1)(x2+y2)n−1 xdy−ydx y xxdy−ydx aydx+bxdy xa−1yb−1 xa−1yb−1(aydx+bxdy) = d(xayb) xy and tan−1 x = x2+y2 NOTE:Ifanon-exact equation has a solution, then an integrating factor is guaranteed to exist, but that does not mean it is easy to find. math168 - solutions to odes - exact differential equations and integration factors 4 Example 4. Solve the differential equation 8 8 dy Example 5. Solve (y2 −y)dx+xdy = 0 (3xy+y2)+(x2+xy) =0 dx Solution Solution M(x,y) = 3xy+y2, N(x,y) = x2+xy M(x,y) = y2 −y, N(x,y) = x M =3x+2yandN = 2x+yClearly, ∂M 6= ∂N,thereforethe M =2y−1, N =1 y x ∂y ∂x y x differential equation is not exact. Involve (16): Mµ − Nµ = µ(N − Hence the differential equal is not y x x exact and no integrating factor is My). Let µ be a function of x alone, that is, µ = µ(x). Applying 1 immediately apparent. Note, however, M −N 3x+2y−(2x+y) 1 that if terms are strategically regrouped, ξ(x) = y x = = the differential equation can be written N x2 +xy x as: R ξ(x)dx R 1dx lnx −(ydx−xdy)+ydx (17) µ(x) = e =e x =e =x The group of terms in parenthesis has many integrating factors in table . Multiplying µ(x) = x through the differential equation we obtain: (above). Trying each integrating factor separately, we find that the only one (3x2y+xy2)+(x3+x2y)dy =0 that makes the entire equating exact is: dx µ(x,y) = 1 , y2 which is a new equation in same form i.e.: M(x,y)+ N(x,y)y′ = 0 Using this integrating factor we can M(x,y) = (3x2y+xy2); N(x,y) = (x3+x2y) rewrite (17) as: −ydx−xdy +1dx=0 (18) M =3x2+2xy; N =3x2+2xy y2 y x which is exact. Let f = 3x2y+xy2. Integrating: Since (18) is exact, it can now be solved. x Alternatively, we note from table(above) Z that (18) can be rewritten as −d(x/y)+ f (x, y) = (3x2y+xy2)dx = x3y+ 1x2y2+ψ(y) 1dx = 0 or as dx = 1dx. Integrating, 2 y weobtain the solution: Find the derivative w.r.t y and equating with N(x,y) we obtain x = x+c or y = x y x+c =⇒ fy = x3+x2yψ′(y) = x3+x2y ψ′(y) = 0 =⇒ ψ(y)=c. Finally x3y+ 1x2y2 = c 2 Linear Equations and those Reducible to the form Afirst-order linear equation is defined as a differential of the form: A(x)y′ +B(x)y = C(x) (19) where A(x) 6= 0 such that the derivative of the dependent variable exist otherwise it is not a differential equation. Then (16) in the stan- dard form is written as y′ + p(x)y = q(x) (20) p(x) = B(x) and q(x) = C(x) A(x) A(x)
no reviews yet
Please Login to review.