Google Search in:

2280 12:55pm Lectures Week 3 S2018

Last Modified: January 22, 2018, 12:45 MST.    Today: December 14, 2018, 01:28 MST.
  Sections 2.3, 2.4, 2.5, 2.6
  The textbook topics, definitions, examples and theorems
Edwards-Penney 1.2, 1.3, 1.4, 1.5 (14.9 K, txt, 17 Jan 2018)
Edwards-Penney 2.1, 2.2, 2.3 (15.8 K, txt, 17 Dec 2014)
Edwards-Penney 2.4, 2.5, 2.6 (11.3 K, txt, 18 Dec 2013)
PDF: Week 3 Examples (98.3 K, pdf, 17 Jan 2018)

Monday: Linear Differential Equations

Linear Differential Equation y'+p(x)y=q(x)
   Review of Section 1.5
     Definition: Linear DE y'+p(x)y=q(x)
     Test: y'=f(x,y) is linear if and only if the partial
           derivative f_y is independent of y.
      Testing linear DE y'=f(x,y) by f_y independent of y.
      Classifying linear equations and non-linear equations.
   Picard's theorem implies a linear DE has a unique solution.
      We don't check for equilibrium solutions or exceptions.
   THEOREM. A linear DE has an explicit general solution.
   Def. Integrating factor W=exp(Q(x)), where Q(x) = int(p(x),x)
   THEOREM. The integrating factor fraction (Wy)'/W replaces the
            two-termed expression y'+py.
   Application Examples: y'+2y=1 and y'+y=e^x.
     How to solve a linear differential equation
       Test the DE for linear
       Identify p(x), q(x) in the standard form y'+py=q.
       Determine an integrating factor W(x)=exp(int(p(x),x))
       Replace y'+py in the standard form y'+py=q by the quotient
          (Wy)' / W
       and then clear fractions to get the quadrature equation
           (Wy)' = qW
       Solve by the method of quadrature.
       Divide by W to find an explicit solution y(x).
       y'+(1/x)y=1, y'+y=e^x, y'+2y=1, y'+(1+x^2)y=0
   Classification: separable, quadrature, linear.
  Superposition Theory
    Superposition for y'+p(x)y=0.
    Superposition for y'+p(x)y=q(x)
  Methods for solving first order linear equations:
     Failsafe: Linear integrating factor method,
     Shortcut: Superposition + equilibrium solution for
         constant-coefficient linear DE
  Fastest way to solve y'+2y=1 is the Shortcut.
   Drill: The example y' + 2y = 6 and more.
TEXT: Examples 1st order linear DE, integrating factor methods (4.3 K, txt, 14 Dec 2014) Solved in class y'+3y=6, y'+y=e^x, and homogeneous equations like y'+3y=0, y'+2y=0. Solve for equilibrium solutions in strange examples like 2y' + Pi y = e^2. Examples and Applications Growth-Decay model y'=ky and its algebraic model y=y(0)exp(kx). Pharmacokinetics of drug transport [PK models, ibuprofen] Pollution models. Three lake pollution model [Erie, Huron, Ontario]. Brine tanks. One-tank model. Two-tank and three-tank models. Chemical engineering example, 3 tanks. See the example:
PDF: Distillation Column (190.5 K, pdf, 22 Jan 2016) Recycled brine tanks and limits of chapter 1 methods. Linear cascades and how to solve them. Method 1: Linear integrating factor method. Method 2: Superposition and equilibrium solutions for constant-coefficient y'+py=q. Uses the shortcut for homogeneous DE y'+py=0. LR and RC circuits Ref: EPbvp 3.7 Basic information on LR, RC and LC circuits.
LR circuit Here
RC circuit Here
LC circuit Here LR Circuit LI' + RI = E RC Circuit RQ' + Q/C = E Method: Linear integrating factor method, usually the shortcut for constant equations. Problem 1.5-34 The expected model is x'=1/4-x/16, x(0)=20, using units of days and millions of cubic feet. The answer is x(t)=4+16 exp(-t/16). Model Derivation Law: x'=input rate - output rate. Definition: concentration == amt/volume. Example: 2kg salt in 15 liters brine, conc=(2/15)kg/liter Use of percentages 0.25% concentration means 0.25/100 concentration References for linear applications
Manuscript: Applications of linear DE (484.2 K, pdf, 16 Jan 2014)
Slides: Brink tanks (86.6 K, pdf, 14 Mar 2016)
Slides: Home heating (99.3 K, pdf, 14 Mar 2016)

Monday-Tuesday: Linear DE. Autonomous systems Section 2.1

Introduction to Ch 2 topics
  2.1, 2.2: Autonomous DE y'=f(y)
    Solution of the Verhulst DE y'=(a-by)y
  2.3: Newton models, Jules Verne problem
  2.4,2.5.26: Numerical solutions of DE.
     No exercises, but a numerical project based on EPH exercises.
  Solution of various exponential models
Manuscript: Exponential Application Library (340.2 K, pdf, 21 Jan 2015) Drill on 1.5 Problems, as time allows There are two special methods for solving y'+py=q If p,q are constant then use the SUPERPOSITION METHOD y = y_p + y_h y_p = an equilibrium solution (set y'=0, solve for y) y_h = constant divided by the integrating factor If one of p or q depends on x, then use the STANDARD METHOD Replace the LHS, which is y'+p(x)y, by the integrating factor quotient (Wy)'/W, where W=exp(int p(x)dx)) is the integrating factor. Cross-multiply by W to clear fractions. Then apply the method of quadrature.
Slides: Linear integrating factor method (126.0 K, pdf, 15 Dec 2014)
Manuscript: Linear DE part I. Integrating Factor Method (303.6 K, pdf, 16 Jan 2014) General Verhulst DE, Section 2.1 Solving y'=(a-by)y by a substitution Let u=y/(a-by). Then substitution into the DE gives u'=au Solve u'=au to get u=u0 exp(ax). Back-substitute u(x) into u=y/(a-by), then solve for y. Solving y'=(a-by)y by partial fractions Divide the DE by (a-by)y Apply the method of quadrature. Find the constants in the partial fractions on the left. Integrate to get the answer a y0 y(x) = -------------------------- b y0 + (a - b y0) exp(-ax) where y0=y(0)=initial population size. Review of Partial fractions. Delayed to week 3. DEFINITION: partial fraction=constant/polynomial with exactly one root THEOREM: P(x)/q(x) = a sum of partial fractions Finding the coefficients. Method of sampling clear fractions, substitute samples, solve for A,B, ... Method of atoms clear fractions, multiply out and match powers, solve for A,B,... Heaviside's cover-up method partially clear fraction, substitute root, find one constant References on partial fractions
Slides: Partial Fraction Theory (148.6 K, pdf, 14 Dec 2014)
Manuscript: Heaviside coverup partial fraction method (290.2 K, pdf, 07 Jan 2014)
Manuscript: Heaviside's method and Laplace theory (352.3 K, pdf, 07 Jan 2014) Separation of variables and partial fractions (Delayed to Monday) Exercise solution problem 2.1-8 The equation y'=7y(y-13), y(0)=17 F(x) = 7, G(y) = y(y-13) Separated form y'/G(y) = F(x) Answer check using the Verhulst solution P(t) = aP_0/(bP_0 + (a-b P_0)exp(-at)) Separation of variables details. Partial fraction details for 1/((u(u-13)) = A/u + B/(u-13)

Fri-Mon: Autonomous Differential Equations and Phase Diagrams. Section 2.2

Lecture on 2.2:
  Theory of autonomous DE y'=f(y)
     Picard's theorem and non-crossing of solutions.
     Direction fields and translation of solutions
  Constructing Euler's threaded solution diagrams
     No direction field is needed to draw solution curves
       We throw out the threaded solution rule used in chapter 1,
       replace it by two rullectureslides/es from calculus and a theorem:
          1. If y'(x)>0, then y(x) increases.
          2. If y'(x)<0, then y(x) decreases.
          THEOREM. For y'=f(y), a threaded solution starting with
              y'(0)>0 must satisfy y'(x)>0 for x>0. A similar result
              holds for y'(0)<0.
     Definition: phase line diagram, phase diagram,
       Calculus tools: f'(x) pos/neg ==> increasing/decreasing
       DE tool: solutions don't cross
       Maple tools for production work.
  Stability theory of autonomous DE y'=f(y)
    Stability of equilibrium solutions.
    Stable and unstable classification of equilibrium solutions.
    funnel, spout, node,
  How to construct Phase line diagrams
  How to make a phase diagram graphic
    Inventing a graph window
    Invention of the grid points
    Using the phase line diagram to make the graphic
        calculus tools
        DE tools
 References for 2.1, 2.2, 2.3. Includes the rabbit problem, partial
 fraction examples, phase diagram illustrations.
Slides: Autonomous DE (101.8 K, pdf, 14 Jan 2014)
Manuscript: Verhulst logistic equation (115.5 K, pdf, 16 Jan 2018)
Manuscript: fish farming (384.5 K, pdf, 16 Jan 2014)
Manuscript: Phase Line and Bifurcation Diagrams. Includes Stability, Funnel, Spout, and bifurcation (765.5 K, pdf, 14 Jan 2014)
Text: Problem notes Chapter 2 (10.8 K, txt, 22 Dec 2014)
Transparencies: ch2 sections 1,2,3: 2.1-6,16,38, 2.2-4,10, 2.3-9,27+Escape velocity (357.6 K, pdf, 29 Jan 2006)
Text: ch2 DEplot maple example 1 for exercises 2.2, 2.3 (0.7 K, txt, 06 Sep 2009)
Text: ch2 DEplot maple example 2 for exercises 2.2, 2.3 (0.7 K, txt, 06 Sep 2009)
html: Problem notes S2018 (2.6 K, html, 07 Nov 2017) Midterm 1 sample exam is at the course web site:
PDF: Sample Exam 1 S2018 with solutions (169.4 K, pdf, 03 Nov 2017)