Differential Equations Notes

Differential Equations

9.1 Modeling with Differential Equations

  • A mathematical model often takes the form of a differential equation, which contains an unknown function and some of its derivatives.

Models of Population Growth
  • A population growth model assumes the population grows at a rate proportional to its size, which is reasonable for bacteria or animals under ideal conditions.

  • Variables:

    • tt = time (independent variable)

    • PP = number of individuals in the population (dependent variable)

  • The rate of growth of the population is the derivative dPdt\frac{dP}{dt}.

  • The assumption that the growth rate is proportional to population size can be written as: dPdt=kP\frac{dP}{dt} = kP, where k is the proportionality constant.

    • This is a differential equation because it contains an unknown function PP and its derivative dPdt\frac{dP}{dt}.

  • If k > 0, then \frac{dP}{dt} > 0 for all tt, meaning the population is always increasing.

  • As P(t)P(t) increases, dPdt\frac{dP}{dt} becomes larger, so the growth rate increases as the population increases.

  • A solution to the equation dPdt=kP\frac{dP}{dt} = kP is an exponential function of the form P(t)=CektP(t) = Ce^{kt}.

    • Taking the derivative: P(t)=C(kekt)=k(Cekt)=kP(t)P'(t) = C(ke^{kt}) = k(Ce^{kt}) = kP(t)

  • The family of solutions is P(t)=CektP(t) = Ce^{kt}.

  • Since populations have only positive values, we are interested in solutions with C > 0 and t > 0.

  • When t=0t = 0, P(0)=Cek(0)=CP(0) = Ce^{k(0)} = C, so CC is the initial population, P(0)P(0).

  • A more realistic model must account for limited resources, where the population levels off when it approaches its carrying capacity MM.

  • Assumptions for a refined model:

    • If PP is small, the growth rate is proportional to PP.

    • If P > M, then PP decreases if it ever exceeds MM.

  • A simple expression incorporating both assumptions is: dPdt=kP(1PM)\frac{dP}{dt} = kP(1 - \frac{P}{M})

  • If PP is small compared with MM, then PM\frac{P}{M} is close to 0, so dPdtkP\frac{dP}{dt} \approx kP.

  • If P > M, then 1PM1 - \frac{P}{M} is negative, so \frac{dP}{dt} < 0.

  • The equation dPdt=kP(1PM)\frac{dP}{dt} = kP(1 - \frac{P}{M}) is called the logistic differential equation.

  • The constant functions P(t)=0P(t) = 0 and P(t)=MP(t) = M are solutions because one of the factors on the right side of the equation is zero; these are equilibrium solutions.

  • If the initial population P(0)P(0) lies between 0 and MM, then \frac{dP}{dt} > 0 and the population increases.

  • If the population exceeds the carrying capacity (P > M), then 1PM1 - \frac{P}{M} is negative, so \frac{dP}{dt} < 0 and the population decreases.

  • If the population approaches the carrying capacity (PMP \rightarrow M), then dPdt0\frac{dP}{dt} \rightarrow 0, which means the population levels off.

A Model for the Motion of a Spring
  • Consider the motion of an object with mass mm at the end of a vertical spring.

  • Hooke’s Law: If the spring is stretched (or compressed) xx units from its natural length, then it exerts a force that is proportional to xx: restoring force =kx= -kx, where kk is a positive constant (the spring constant).

  • Ignoring external resisting forces, by Newton’s Second Law (force equals mass times acceleration), we have md2xdt2=kxm\frac{d^2x}{dt^2} = -kx.

  • This is a second-order differential equation because it involves second derivatives.

  • The equation can be rewritten as: d2xdt2=kmx\frac{d^2x}{dt^2} = -\frac{k}{m}x, which means the second derivative of xx is proportional to xx but has the opposite sign.

General Differential Equations
  • A differential equation is an equation that contains an unknown function and one or more of its derivatives.

  • The order of a differential equation is the order of the highest derivative that occurs in the equation.

    • dPdt=kP\frac{dP}{dt} = kP and dPdt=kP(1PM)\frac{dP}{dt} = kP(1 - \frac{P}{M}) are first-order equations.

    • md2xdt2=kxm\frac{d^2x}{dt^2} = -kx is a second-order equation.

  • The independent variable (e.g., tt) doesn’t always have to represent time.

  • For example, in the differential equation y=xyy' = xy, yy is an unknown function of xx.

  • A function ff is a solution of a differential equation if the equation is satisfied when y=f(x)y = f(x) and its derivatives are substituted into the equation.

    • ff is a solution of y=xyy' = xy if f(x)=xf(x)f'(x) = xf(x) for all values of xx in some interval.

  • To solve a differential equation means to find all possible solutions of the equation.

  • We have already solved simple differential equations of the form y=f(x)y' = f(x).

    • The general solution of y=x3y' = x^3 is y=x44+Cy = \frac{x^4}{4} + C, where CC is an arbitrary constant.

Example 1
  • Show that every member of the family of functions y=1+cet1cety = \frac{1 + ce^t}{1 - ce^t} is a solution of the differential equation y=12(y21)y' = \frac{1}{2}(y^2 - 1).

Solution

  • Using the Quotient Rule to differentiate yy:
    y=(1cet)(cet)(1+cet)(cet)(1cet)2=cetc2e2t+cet+c2e2t(1cet)2=2cet(1cet)2y' = \frac{(1 - ce^t)(ce^t) - (1 + ce^t)(-ce^t)}{(1 - ce^t)^2} = \frac{ce^t - c^2e^{2t} + ce^t + c^2e^{2t}}{(1 - ce^t)^2} = \frac{2ce^t}{(1 - ce^t)^2}

  • The right side of the differential equation becomes:
    12(y21)=12[(1+cet1cet)21]=12[(1+cet)2(1cet)2(1cet)2]=12[1+2cet+c2e2t(12cet+c2e2t)(1cet)2]=12[4cet(1cet)2]=2cet(1cet)2\frac{1}{2}(y^2 - 1) = \frac{1}{2}[(\frac{1 + ce^t}{1 - ce^t})^2 - 1] = \frac{1}{2}[\frac{(1 + ce^t)^2 - (1 - ce^t)^2}{(1 - ce^t)^2}] = \frac{1}{2}[\frac{1 + 2ce^t + c^2e^{2t} - (1 - 2ce^t + c^2e^{2t})}{(1 - ce^t)^2}] = \frac{1}{2}[\frac{4ce^t}{(1 - ce^t)^2}] = \frac{2ce^t}{(1 - ce^t)^2}

  • Therefore, for every value of cc, the given function is a solution of the differential equation.

  • When applying differential equations, we are usually interested in finding a particular solution that satisfies an additional requirement, rather than a family of solutions (the general solution).

  • In many physical problems, we need to find the particular solution that satisfies a condition of the form y(t<em>0)=y</em>0y(t<em>0) = y</em>0.

    • This is called an initial condition, and the problem of finding a solution of the differential equation that satisfies the initial condition is called an initial-value problem.