M244 2016 Solutions to Final Exam Review

M244: Solutions to Final Exam Review 1. Solve (use any method if not otherwise specified): (a) (2x − 3x2 ) dx dt = t cos(t) SOLUTION: As a separable DE: Z Z 2x − 3x2 dx = t cos(t) dt ⇒ x2 − x3 = cos(t) + t sin(t) + C You can leave your answer in implicit form. (b) y 00 + 2y 0 + y = sin(3t) SOLUTION: Get the homogeneous part then the particular solution (or use Laplace): r2 + 2r + 1 = 0 ⇒ (r + 1)2 = 0 ⇒ r = −1, −1 ⇒ yh (t) = e−t (C1 + C2 t) For the particular solution, (Undet Coefs), yp = A cos(3t) + B sin(3t). Substitute to get: (A − 6B − 9A) sin(3t) + (B + 6A − 9B) cos(3t) = sin(3t) so that A = −3/50, B = −2/25. Altogether, y(t) = e−t (C1 + C2 t) − 2 3 cos(3t) − sin(3t) 50 25 (c) y 00 − 3y 0 + 2y = e2t SOLUTION: Same technique here. The roots to the characteristic equation are r = 1, 2, so the homogeneous part of the solution is: yh (t) = C1 et + C2 e2t Initially, we guess that yp (t) = Ae2t , but that is part of yh , so multiply by t: yp = Ate2t . Now substitute into the D.E. to get: A = 1. The full solution is y(t) = C1 et + C2 e2t + te2t √ (d) y 0 = te−t − y. √ SOLUTION: This is Ra linear differential equation, with integrating factor; y 0 + y = te−t . The integrating factor is e 1 dt = et . Therefore,   √ 2 2 3/2 t + C e−t (yet )0 = t ⇒ yet = t3/2 + C ⇒ y = 3 3 (e) x0 = 2 + 2t2 + x + t2 x. SOLUTION: This is a linear differential equation: x0 − (1 + t2 )x = 2(1 + t2 ). The integrating factor is: R 2 3 e −(1+t ) dt = e−t−(1/3)t So we solve the following (to integrate, let u = t + (1/3)t2 )  3 xe−t−(1/3)t 0   3 3 = 2(1 + t2 )e−t−(1/3)t ⇒ xe−t−(1/3)t = 3 −2e−t−(1/3)t + C ⇒ x = −2 + Cet+(1/3)t (f) 3 x01 = 2x1 + 3x2 (Typo corrected) x02 = 4x1 + x2 SOLUTION: The characteristic equation is λ2 − 3λ − 10 = 0, which factors as (λ − 5)(λ + 2) = 0. Solving for eigenvectors: 1 • For λ = −2, we have (2 − −2)v1 + 3v2 = 0, or v1 = (−3/4)v2 . Putting this in vector form, we simplify it by multiplying by 2: v1 = [−3, 4]T . • For λ = 5, we have (2 − −5)v1 + 3v2 = 0, or v1 = v2 . Putting this answer in vector form, we have v2 = [1, 1]T . Overall, the general solution is x(t) = C1 e−2t   −3 2 + C2 e5t  1 1  (g) t2 y 00 − 2ty 0 − 10y = 0 This is an Euler equation. We can take the ansatz y = tr to get the characteristic equation: r(r − 1) − 2r − 10 = 0 ⇒ r = −5, 2 The solution is: y(t) = C1 t5 + C2 t−2 2. Show that with the proper substitution, the following equation becomes separable. NOTE: You do not need to solve the differential equation. 3x − 4y dy = dx y − 2x SOLUTION: Divide numerator and denominator by x, and let v = xy (the equation is a homogeneous equation), therefore y = xv and y 0 = xv 0 + v. Although we don’t need to solve the DE, we should go ahead and separate the variables. xv 0 + v = 3 − 4v v−2 ⇒ xv 0 = 3 − 4v − v(v − 2) 3 − 2v − v 2 = v−2 v−2 ⇒ v−2 1 dv = dx 2 3 − 2v − v x 3. Show that with the proper substitution, the following equation becomes linear. NOTE: You do not need to solve the differential equation. dy x + 3xy = 2 dx y SOLUTION: This is a Bernoulli equation. Multiplying through by y 2 gives: y 2 y 0 + 3xy 3 = x Let v = y 3 . Then v 0 = 3y 2 y 0 , or 13 v 0 = y 2 y 0 . Therefore, we get: 1 0 v + 3xv = x 3 ⇒ v 0 + 9xv = 3x which is linear in v. 4. Obtain the general solution in terms of α, then determine a value of α so that y(t) → 0 as t → ∞: SOLUTION: y 00 − y 0 − 6y = 0, y(0) = 1, y 0 (0) = α The general solution (before initial conditions): y(t) = C1 e3t + C2 e−2t With the initial conditions, 1 = C1 + C2 α = 3C1 − 2C2 ⇒ C1 = Therefore, 2+α , 5 2 + α 3t 3 − α −2t e + e 5 5 For y(t) → 0, we must have α = −2 (to zero out the first term). y(t) = 2 C2 = 3−α 5 5. If y 0 = y(1 − y)(2 − y)(3 − y)(4 − y) and y(0) = 2.5, determine what y does as t → ∞. SOLUTION: The idea here is to figure out where the equilibria are and if they are stable or unstable. This is a fifth degree polynomial which increases after y = 4 and decreases to −inf as y decreases. From a quick sketch of the polynomial, we see that y = 1, 3 are stable and y = 0, 2, 4 are unstable equilibria. If the solution starts at y = 2.5, it will increase to the next stable equilibrium y = 3. 6. If y1 , y2 are a fundamental set of solutions to t2 y 00 − 2y 0 + (3 + t)y = 0 and if W (y1 , y2 ) = 3, find W (y1 , y2 )(4). SOLUTION: Before using Abel’s Theorem, put the equation in standard form as: y 00 − 2 0 3+t y + 2 y=0 t2 t so that the Wronskian between any two solutions is: R −2 Ce2 t dt = Ce−2/t Given that the Wronskian is 3 at t = 2, we have: Ce−1 = 3 and now ⇒ C = 3e √ W (y1 , y2 )(4) = 3ee−2/4 = 3 e 7. What is the Wronskian? How is it used? SOLUTION: The Wronskian is an operation that is performed on two functions f, g: W (f, g) = f (t) g(t) f 0 (t) g 0 (t) The Wronskian is primarily used to determine if two solutions to a second order linear homogeneous differential equation form a fundamental set of solutions- That is, if an arbitrary solution to an IVP could be written as a linear combination of the functions in the set. 8. Explain Abel’s Theorem: SOLUTION: Abel’s Theorem states that, given any two solutions to a second order linear homogeneous DE (on the interval on which the existence and uniqueness theorem holds): y 00 + p(t)y 0 + q(t)y = 0 Then the Wronskian can be computed by: W (y1 , y2 )(t) = Ce− R p(t) dt This implies that either the Wronskian is always zero on the interval, or never zero (the two functions form a fundamental set). 9. Give the three Existence and Uniqueness Theorems we have had in class. SOLUTION: We had the general existence and uniqueness theorem, then one for linear first order, then one for linear second order equations: • Let y 0 = f (t, y) with y(t0 ) = y0 . If there is a open rectangle R that contains (t0 , y0 ), and on which f and ∂f /∂y are continuous, then there exists an  > 0 for which a unique solution exists to the DE, valid for t0 −  < t < t0 + . 3 • Let y 0 + p(t)y = g(t), with y(t0 ) = y0 . If there is an open interval I containing t0 on which both p, g are continuous, then there is a solution to the IVP, valid on all of I. • Let y 00 + p(t)y 0 + q(t)y = g(t), with y(t0 ) = y0 , y 0 (t0 ) = v0 . If there is an open interval I containing t0 on which both p, q, g are continuous, then there is a solution to the IVP, valid on all of I. 10. Let y 00 − 6y 0 + 9y = F (t). For each F (t) listed, give the form of the general solution using undet. coeffs (do not solve for the coefficients). SOLUTION: Before we start, we should go ahead and solve the homogeneous equation: r2 − 3r + 9 = 0 (r − 3)2 = 0 ⇒ ⇒ r = 3, 3 Therefore, yh = C1 e3t + C2 te3t (a) F (t) = 2t2 SOLUTION: Guess the full quadratic, yp = At2 + Bt + C (b) F (t) = te−3t sin(2t): SOLUTION: We need a polynomial of degree 1 with each sine and cosine: yp (t) = e−3t ((At + B) cos(2t) + (Ct + D) sin(2t)) (c) F (t) = t sin(2t) + cos(2t) SOLUTION: Similar to the last one, but we can take both sine and cosine together: yp = (At + B) cos(2t) + (Ct + D) sin(2t) (d) F (t) = 2t2 + 12e3t SOLUTION: Break this one up into two solutions- One for g1 = t2 and one for g2 = 12e3t . Then: yp1 (t) = At2 + Bt + C yp2 = At2 e3t where I multiplied the second guess by t2 so that it looks like no term of the homogeneous part of the solution. 11. Newton’s Law of Cooling states that the rate of change of the temperature of a body is proportional to the difference between the body’s temperature and the temperature of the environment. In notation, if u(t) is the temperature at time t, and T is the (constant) environmental temperature, then: du = −k(u − T ) dt (The negative sign out front implies that k can be assumed to be positive). 12. A spring is stretched 0.1 m by a force of 3 N. A mass of 2 kg is hung from the spring and is also attached to a damper that exerts a force of 3 N when the velocity of the mass is 5 m/s. If the mass is pulled down 0.05 m below its resting equilibrium and released with a downward velocity of 0.1 m/s, determine its position u at time t. SOLUTION: From what is given, we can compute the spring constant. mg − kL = 0 ⇒ 3− k =0 10 ⇒ k = 30 Now, we also assume that damping is proportional to the velocity, so: γv = 5γ = 3 ⇒ γ= 3 5 Now we can construct the model (downward is positive): 3 2u00 + u + 30u = 0 5 u(0) = 0.05 u0 (0) = 0.1 We’ll go ahead and use some numerical approximations for this- The roots of the characteristic equation are approximately: r = −0.15 ± 3.87 i 4 (Sorry! This came out of Section 3.7 so I thought the numbers would work out better...) Therefore, the solution is: e−0.15t (C1 cos(3.87 t) + C2 sin(3.87 t)) If we put in the initial conditions, we get: − 0.15C1 + 3.87C2 = 0.1 C1 = 0.05 ⇒ C2 ≈ 0.065 which gives the solution. 13. Let y(x) be a power series solution to y 00 − xy 0 − y = 0, x0 = 1. Find the recurrence relation and write the first 5 terms of the expansion of y. Same idea as the previous exercise, but be careful to use powers of (x − 1)! ∞ X n(n − 1)cn−2 (x − 1)n−2 − x n=2 ∞ X ncn (x − 1)n−1 − n=1 ∞ X cn (x − 1)n = 0 n=0 For the middle sum, recall our trick: x = 1 + (x − 1), which gives us a way of incorporating the x into the sum: ∞ ∞ ∞ ∞ X X X X x ncn (x − 1)n−1 = (1 + (x − 1)) ncn (x − 1)n−1 = ncn (x − 1)n−1 + ncn (x − 1)n n=1 n=1 n=1 n=1 Now shift the index of every sum to match, and you should get the recurrence relation: 1 (Ck+1 + Ck ) k+2 Ck+2 = We get: C0 + C1 C0 + 3C1 2C0 + 3C1 (x − 1)2 + (x − 1)3 + (x − 1)4 + . . . 2 6 12 We can check our answer using the next exercise... y(x) = C0 + C1 (x − 1) + 14. Let y(x) be a power series solution to y 00 − xy 0 − y = 0, x0 = 1 (the same as the previous DE), with y(1) = 1 and y 0 (1) = 2. Compute the first 5 terms of the power series solution by first computing y 00 (1), y 000 (1), y (4) (1). First, let’s compute the derivatives: y 00 = xy 0 + y at x = 1 ⇒ y 00 (1) = y 0 (1) + y(1) y 000 = xy 00 + 2y 0 at x = 1 ⇒ y 000 (1) = y 00 (1) + 2y 0 (1) y (4) = xy 000 + 3y 00 at x = 1 ⇒ y (4) (1) = y 000 (1) + 3y 00 (1) so that: and: From this, we see that: y 00 (1) = 3, y 000 (1) = 7, y (4) (1) = 16 Writing out the solution: 3 7 16 y(x) = 1 + 2(x − 1) + (x − 1)2 + (x − 1)3 + (x − 1)4 + . . . 2 6 24 15. Use the definition of the Laplace transform to determine L(f ):  3, 0≤t≤2 f (t) = 6 − t, 2 < t SOLUTION: Z ∞ f (t)e 0 −st Z dt = 3 2 e −st Z dt + ∞ (6 − t)e−st dt = 2 0 3 e−2s (1 − e−2s ) + 2 (4s − 1) s s 5 16. Determine the Laplace transform: (a) t2 e−9t ⇒ 2 (s + 9)3 e−5s s22 + 4 (c) e3t sin(4t) ⇒ (s − 3)2 + 16 (b) u5 (t)(t − 2)2 ⇒ 6 s2 + 9 s
(d) et δ(t − 3) ⇒ e−3s+3 17. Find the inverse Laplace transform: √   √ √ 2 3 2s − 1 s−2 3 2t (a) 2 . Rewrite: 2· +√ · The inverse is then e 2 cos( 2t + √ sin( 2t) s − 4s + 6 (s − 2)2 + 2 2 (s − 2)2 + 2 2 7 7 2 −3t (b) ⇒ t e (s + 3)3 2 (c) e−2s (4s + 2) . You might rewrite this as e−2s H(s), where (s − 1)(s + 2) H(s) = 4s + 2 2 2 = + (s − 1)(s + 2) s+2 s−1 Now, h(t) = 2e−2t + 2et , and the solution is u2 (t)h(t − 2). 3s − 2 (d) We might rewrite this as: (s − 4)2 − 3 √ s−4 10 10 3 3· +√ · = 3F (s − 4) + √ G(s − 4) 2−3 (s − 4)2 − 3 (s − 4) 3 3 √ s 3 where F (s) = 2 , G(s) = 2 . The inverse is (Item 14 from the Table): s −3 s −3     √ √ 10 10 e4t 3f (t) + √ g(t) = e4t 3cosh( 3t) + √ sinh( 3t) 3 3 18. Solve the given initial value problems using Laplace transforms: (a) y 00 + 2y 0 + 2y = 4t, y(0) = 0, y 0 (0) = −1. The Laplace transform: Y (s) = 4 − s2 2 2 2s + 1 =− + 2 + 2 = + 2s + 2) s s s + 2s + 2 s2 (s2 2 2 s+1 1 − + 2 +2 − s s (s + 1)2 + 1 (s + 1)2 + 1 so that y(t) = −2 + 2t + e−t (2 cos(t) − sin(t)) (b) y 00 − 2y 0 − 3y = u1 (t), y(0) = 0, y 0 (0) = −1 Use partial fractions:   1 1 1 1 1 1 1 1 1 1 Y (s) = − · + · + e−s − · + · + · 4 s+1 4 s−3 3 s 4 s + 1 12 s − 3 Think of this second term as e−s · H(S), where 1 1 1 h(t) = − + e−t + e3t 3 4 12 and the solution is: 1 1 y(t) = − e−t + e3t + u1 (t)h(t − 1) 4 4 6 (c) y 00 − 4y 0 + 4y = t2 et , y(0) = 0, y 0 (0) = 0. Y (s) = 2 (s − 2 t 1)3 (s − 2)2 2 1 · = F (s)G(S) 3 (s − 1) (s − 2)2 = 2t Where f (t) = t e and g(t) = te . Therefore, y(t) = t2 et ∗ te2t 19. Consider t2 y 00 − 4ty 0 + 6y = 0 (a) Using y1 = t2 as one solution, find y2 by computing the Wronskian two ways. SOLUTION: We are verifying the solution above by using the Wronskian. Using this technique, we have: t2 y2 = t2 y20 − 2ty2 2t y20 And, using Abel’s Theorem (with p(t) = − 4t ): W = Ce R 4 t dt = Ct4 Therefore, t2 y20 − 2ty2 = Ct4 This has an integrating factor of t−2 :  y 0 2 =C ⇒ t2 ⇒ y2 = Ct + C1 t2 2 y20 − y2 = Ct2 t ⇒ y2 = Ct3 + C1 t2 We already have y1 = t2 , so y2 = t3 (as before). (b) Using y1 = t2 as one solution, find y2 by using reduction of order. SOLUTION: In the reduction of order, we use the following substitutions: y2 = vy1 = t2 v y20 = t2 v + 2tv y200 = t2 v 00 + 4tv 0 + 2v Substituting into the DE gives: t2 y200 −4ty20 +6y2 0 = t4 v 00 + 4t3 v 0 + =− 4t3 v 0 − = = t4 v 00 2t2 v 8t2 v 6t2 v ⇒ v 00 = 0 ⇒ v = C1 t + C2 from which: y2 = t2 v = C1 t3 + C2 t2 and again, since y1 = t2 , we only require that y2 = t3 . 20. For the following differential equations, (i) Give the general solution, (ii) Solve for the specific solution, if its an IVP, (iii) State the interval for which the solution is valid. y(0) = 1. This is a linear (integrating factor) differential equation; the solution (a) y 0 − 21 y = e2t will be valid for all time t. y(t) = 7 2 2t 1 1 t e + e2 3 3 (b) y 0 = 12 y(3 − y). SOLUTION: Separable. We’ll need partial fractions to integrate. Z 1 1 3 1 1 1 · + · dy = t + C ⇒ ln |y| − ln |3 − y| = t + C2 3 y 3 3−y 2 2 Now solve for y: 3 y = Ae3/2t ⇒ y = 3−y 1 + Be−3/2t If the initial condition is positive, this is valid for all time (Draw the phase diagram and direction field for this autonomous DE to see why). If the initial condition is negative, we would need to find where (in positive time) the solution has a vertical asymptote. NOTE: We are assuming that t ≥ 0. (c) y 00 + 2y 0 + y = 0, y(0) = α, y 0 (0) = 1 y(t) = αe−t + (1 + α)te−t This solution is valid for all t. (d) 2xy 2 + 2y + (2x2 y + 2x)y 0 = 0 This is an exact equation: ∂ ∂ (2xy 2 + 2y) = 4xy + 2 = (2x2 y + 2x) ∂y ∂x Recall that the solution will be (implicit) F (x, y) = C, where Fx = 2xy 2 + 2y ⇒ F (x, y) = x2 y 2 + 2xy + h(x) and Fy = 2x2 y + 2x ⇒ F (x, y) = x2 y 2 + 2xy + g(y) Comparing, we see F (x, y) = x2 y 2 + 2xy, and the implicit solution is: x2 y 2 + 2xy = C Here we will not be able to give an interval on which the solution is valid unless we isolate y, although we would have a requirement that 2x2 y + 2x 6= 0, so that y 0 would be defined. (e) y 00 + 4y = t2 + 3et , y(0) = 0, y 0 (0) = 1. y(t) = 1 19 1 1 3 sin(2t) − cos(2t) + t2 − + et 5 40 4 8 5 The solution is valid for all t. 21. Suppose y 0 = −ky(y − 1), with k > 0. Sketch the phase diagram. Find and classify the equilibrium. Draw a sketch of y on the direction field, paying particular attention to where y is increasing/decreasing and concave up/down. Finally, get the analytic (general) solution. Your graph should be an upside down parabola (vertex up). There are equilibrium solutions at y = 0 (unstable) and y = 1 (stable). 22. True or False (and explain): Every separable equation is also exact. If true, is one way easier to solve over the other? True: dy = f (y)g(x) dx ⇒ −g(x) dx + 1 dy = 0 f (y) So, letting −g(x) = M (x, y) and 1/f (y) = N (x, y), then My = Nx = 0. Therefore, every separable equation is exact. To see if either is simpler, check to see what integrals must be performed. Treating the equation as separable, we have to compute the two integrals: Z Z dy g(x) dx f (y) 8 Treating the equations as exact, Z Z M (x, y) dx = − Z g(x) dx + H(y) Z N (x, y) dy = dy + G(x) f (y) so we have to compute the same integrals either way (so neither is easier than the other). 23. Let y 0 = 2y 2 + xy 2 , y(0) = 1. Solve, and find the minimum of y. Hint: Determine the interval for which the solution is valid. This is separable: y 0 = y 2 (2 + x) ⇒ y −2 dy = (2 + x) dx, so y(x) = x2 −2 + 4x − 2 √ √ This has √ vertical asymptotes at x = −2 ± 6, so that the solution is valid only when −2 − 6 < x < −2 + 6, or when x is approximately between −4.45 and 0.45. Between these vertical asymptotes, y has a minimum where its derivative is 0, y 0 = y 2 (2 + x) = 0 ⇒ y = 0 or x = −2 From our solution, we see that y 6= 0, so the minimum occurs at x = −2, and the minimum is: 1/3. 24. A sky diver weighs 180 lbs and falls vertically downward. Assume that the constant for air resistance is 3/4 before the parachute is released, and 12 after it is released at 10 sec. Assume velocity is measured in feet per second, and g = 32 ft/sec2 . (a) Find the velocity of the sky diver at time t (before the parachute opens). SOLUTION: Recall our model: mv 0 = mg − γv And, since 180 = mg = 32m, we get that mass is 45/8. Therefore, the model becomes: v 0 = 32 − 3 8 2 · v = 32 − v 4 45 15 v(0) = 0 Then, for the first 10 seconds, the velocity is given by: 240 − 240e−(2/15)t At 10 seconds, the model changes to: v 0 = 32 − 12 · 32 8 v = 32 − v 45 15 The general solution (for t > 10) is: 15 + C1 e−(32/15)t and we would find C1 so that the velocity function is continuous (a bit of a mess algebraically). (b) If the sky diver fell from an altitude of 5000 feet, find the sky diver’s position at the instant the parachute is released. SOLUTION: Solve for C using −5000 (down is positive in our model): s(t) = 1800e−(2/15)t + 240t − 6800 so at t = 10, the position s is approximately: −3925.52, which is about 3926 feet high. (c) After the parachute opens, is there a limiting velocity? If so, find it. (HINT: You do not need to re-solve the DE). SOLUTION: Yes- The equilibrium solution to the model equation (after the parachute is pulled) is 15 feet per second (which we can also see from the solution). 9 25. Rewrite the following differential equations as an equivalent system of first order equations. If it is an IVP, also determine initial conditions for the system. (a) y 00 − 3y 0 + 4y = 0, y(0) = 1, y 0 (0) = 2. SOLUTION: Let x1 = y and x2 = y 0 : x01 x02 = x2 = −4x1 + 3x2 (b) y 000 − 2y 00 − y 0 + 4y = 0 SOLUTION: Same idea as before, but x3 = y 00 . Then: x01 x02 x03  ⇒ x(0) = 1 2  now we need three variables. Let x1 = y, x2 = y 0 and = = = −4x1 x2 x2 x3 +2x3 (c) y 00 − yy 0 + t2 = 0 SOLUTION: This equation is nonlinear and non-autonomous, but the substitutions work just as before: Let x1 = y, x2 = y 0 . Then: x01 = x2 x02 = x1 x2 − t2 26. Convert one of the variables in the following systems to an equivalent higher order differential equation, and solve it (be sure to solve for both x and y): x0 y0 = 4x + y = −2x + y SOLUTION: Solve the first equation for y, substitute into the second equation to get: x00 − 5x0 + 6x = 0 ⇒ (r − 3)(r − 2) = 0 Therefore, x(t) = C1 e3t + C2 e2t . To find y, we have to compute x0 − 4x, which is: y(t) = 3C1 e3t + 2C1 e2t − 4(C1 e3t + C2 e2t ) Therefore, putting these together as a parametric solution:       x 1 1 = C1 e3t + C2 e2t y −1 −2 27. Solve the previous system by using eigenvalues and eigenvectors. SOLUTION: As a hint, we should be able to read the eigenvalues and eigenvectors off of the given solution, but we’ll go through it the long way just to be double check the method.  A= 4 −2 1 1  Tr(A) = 5 det(A) = 6 ⇒ λ2 − 5λ + 6 = 0 And this is indeed the same characteristic function as before. We’ll work with λ = 3 first, and:   (4 − 3)v1 + v2 = 0 v1 = v1 1 ⇒ ⇒ −2v1 + (1 − 3)v2 = 0 v2 = −v1 −1 A similar computation shows that [1, −2]T is an eigenvector for λ = 2. Therefore, we get the same general solution as before. 28. Verify by direct substitution that the given power series is a solution of the differential equation: y= ∞ X (−1)n+1 n x n n=1 10 (x + 1)y 00 + y 0 = 0 SOLUTION: If you get stuck on this one, you might think about what this function is without the sum: ∞ X 1 1 1 1 1 (−1)n+1 n y = x − x2 + x3 − x4 + x5 − x6 + · · · = x 2 3 4 5 6 n n=1 Therefore, ∞ X y 0 = 1 − x + x2 − x3 + x4 − x5 + x6 · · · = (−1)n+1 xn−1 n=1 y 00 = −1 + 2x − 3x2 + 4x3 − 5x4 + · · · = ∞ X (−1)n+1 (n − 1)xn−2 n=2 00 0 So you can check that (1 + x)y + y = 0 for these terms- However, we should be able to show it in general using the sum notation. We need to show that the following polynomial is zero for all x: (1 + x) ∞ X (−1)n+1 (n − 1)xn−2 + n=2 ∞ X (−1)n+1 xn−1 n=1 which is easiest if we have a single sum: ∞ X (−1)n+1 (n − 1)xn−1 + n=2 ∞ X ∞ X (−1)n+1 (n − 1)xn−2 + n=2 (−1)n+1 xn−1 n=1 Combining to have xk with k = 0, 1, 2, · · ·, we have: ∞ X
k(−1)k+2 + (k + 1)(−1)k+3 + (−1)k+2 xk k=0 The coefficients become: (k + 1)(−1)k+2 + (k + 1)(−1)k+3 so, if k is even, the expression is 0 and if k is odd the expression is zero- Every coefficient is zero. 29. Convert the given expression into a single power series: ∞ X n(n − 1)an xn + 2 n=2 ∞ X nan xn−2 + 3 n=2 ∞ X an xn n=1 SOLUTION: The first sum starts with x2 , the second sum starts with x0 and the third sum starts with x1 . It will probably be easiest to make the combined sum start with x1 , so we’ll do that by first making the individual sums do that: " # ∞ ∞ ∞ X X X n n−2 n(n − 1)an x + 4a2 + 2 nan x +3 a n xn n=1 n=3 n=1 If we make m our new index, then in the first and last sums, n = m and in the middle sum, m = n − 2: 4a2 + ∞ X [m(m − 1)am + 2(m + 2)am+2 + 3am ] xm m=1 11 30. Find the recurrence relation for the coefficients of the power series solution to y 00 − (1 + x)y = 0 at x0 = 0. SOLUTION: Substitute the power series into the DE to get: ∞ X n(n − 1)an xn−2 − (1 + x) n=2 ∞ X an xn = 0 n=0 Break up the second sum: ∞ X n(n − 1)an xn−2 − n=2 ∞ X a n xn − n=0 ∞ X an xn+1 = 0 n=0 We want to collect these into a single sum- Looks like we might need to peel the constant terms off of the first and second sums to get all series to start with the first power: ! ∞ ∞ ∞ X X X n−2 n 2 · 1 · a2 + n(n − 1)an x − a0 + an x − an xn+1 = 0 n=3 n=1 n=0 Now we can write everything as a single sum: (2a2 − a0 ) + ∞ X ((k + 2)(k + 1)ak+2 − ak − ak−1 ) xk = 0 k=1 from which we get the recurrence relation: a2 = 1 a0 2 and ak+2 = ak + ak−1 (k + 2)(k + 1) for k = 1, 2, 3, · · · 31. Find the first 5 non-zero terms of the series solution to y 00 − (1 + x)y = 0 if y(0) = 1 and y 0 (0) = −1 (use derivatives). SOLUTION: Taking derivatives, y 00 = (1 + x)y y (3) = y + (1 + x)y 0 y (4) = 2y 0 + (1 + x)y 00 and so on. From this we see that (remember to divide by the appropriate factorial): 1 1 1 y(x) = 1 − x + x2 − x4 + x5 + · · · 2 24 40 EXTRA: You might notice that the derivatives have the relation: y (k+2) (0) = ky (k−1) (0) + y (k) (0) Substituting from the Taylor formula y (k) k! = ak , or y (k) = k!ak we get: (k + 2)!ak+2 = k(k − 1)!ak−1 + k!ak Simplifying, we get ak+2 = ak + ak−1 (k + 1)(k + 2) which you might recognize from the previous problem! 32. We have two tanks, A and B with 20 and 30 gallons of fluid, respectively. Water is being pumped into Tank A at a rate of 2 gallons per minute, 2 ounces of salt per gallon. The well-mixed solution is pumped out of Tank A and into Tank B at a rate of 4 gallons per minute. Solution from Tank B is entering Tank A at a rate of 2 gallons per minute. Water is being pumped into Tank B at k gallons per minute with 3 ounces of salt per gallon. The solution is being pumped out of tank B at a total rate of 5 gallons per minute (2 of them are going into tank A). 12 • What should k be in order for the amount of solution in Tank B to remain at 30? Use this value for the remaining problems. k=1 • Write the system of differential equations for the amount of salt in Tanks A, B at time t. Do not solve. Let A(t), B(t) be the amount (in ounces) of salt in Tanks A, B respectively. Then: 2 4 dA =4+ B− A dt 30 20 dB 4 5 =3+ A− B dt 20 30 • Find the equilibrium solution and classify it. Solve: 1 4 − 51 A + 15 B =0 1 3 + 5 A − 16 B = 0 ⇒ A= 130 3 B = 70 Notice that these quantities are in total ounces in each tank. We can classify this equilibrium by looking at the matrix associated with the system:  1  1 − 5 15 1 − 16 5 The trace is negative, the determinant is positive, and the discriminant is positive. The equilibrium is a SINK (which is what we would expect from the physical problem). 33. Solve, and determine how the solution depends on the initial condition, y(0) = y0 : y 0 = 2ty 2 y(t) = −y0 y0 t2 − 1 If y0 > 0, then the solution will only be valid between ± √1y0 . If y0 < 0, the solution will be valid for all t. 34. Solve the linear system x0 = Ax using eigenvalues and eigenvectors, if A is as defined below:   2 8 (a) A = SOLN: λ = ±2i −1 −2 For λ = 2i, the equation (A − λI)v = 0 simplifies to:   v1 = (−2 − 2i)v2 −2 − 2i ⇒ v= v2 = v 2 1 Now we compute eλt v before writing the solution:     −2 − 2i (−2 cos(2t) + 2 sin(2t)) + i(−2 cos(2t) − 2 sin(2t)) (cos(2t) + i sin(2t)) = 1 cos(2t) + i sin(2t) The solution to the DE is now:     −2 cos(2t) + 2 sin(2t) −2 cos(2t) − 2 sin(2t) x(t) = C1 + C2 cos(2t) sin(2t)   2 3 (b) A = SOLN: λ2 − 3λ − 4 = 0, so λ = −1, 4. 2 1 If λ = −1, then v1 + v2 = 0, or   −1 v= 1 13 If λ = 4, then 2v1 + 3v2 = 0, or  3 2  −1 1  v= Therefore, the general solution is: −t x(t) = C1 e  +e 4t  3 2    3 −18 2 −9 You should see that λ = −3, −3. Putting λ = −3 into the equation for the eigenvector, we get:   6v1 − 18v2 = 0 3 ⇒ quadv = 2v1 − 6v2 = 0 1 (c) A = Computing the generalized eigenvector w, we get: 6w1 − 18w2 2w1 − 6w2 =3 =1  ⇒ quadw = 1/2 0  where we choose one convenient solution for w. The general solution is then:         3 3 1/2 x(t) = e−3t C1 + C2 t + 1 1 0 35. For each system x0 = Ax, the matrix A depends on the parameter α. Find how the classification of the origin changes depending on α. Tr(A) = 0 det(A) = −4 + 5α ∆ = −4(−4 + 5α) SOLUTION: On the Poincare Diagram, we’re located on the determinant axis. Therefore, we have a center if the determinant is positive (α > 4/5), we have “uniform motion” if α = 54 , and we have a saddle if α < 4/5.   Tr(A) = α + 1 α 2 det(A) = α − 6 (b) 3 1 ∆ = α2 − 2α + 25 SOLUTION: We would probably need to do a sign analysis like we did in class, but if we look at the discriminant, we see that it is always positive (it has complex roots). Therefore, we stay outside of the parabola in the Poincare Diagram. We might check the other signs- The trace changes sign at α = −1 and the determinant changes sign at α = 6. Putting these on a number line:  (a) 2 α −5 −2  Tr(A) det(A) −− ++ −− −− α < −1 −1 < α < 6 ++ ++ α>6 Now we can summarize: α < −1 α = −1 −1 < α < 6 α=6 α>6 Saddle (Quadrant III) Saddle Saddle (Quadrant IV) Line of unstable fixed points Source (Quadrant I) 14 15