Tag: Carrying Capacity

  • The Population Model That Fails—and Why Its Equation Is Everywhere: Population, Interest, and Radioactive Decay

    “`html

    A remarkably simple differential equation appears in many different places in the real world. We will begin with population growth. The first model we try will have a serious problem: it predicts unlimited growth. Fixing that problem will lead us to the logistic equation.

    Then, surprisingly, we will return to our original equation and discover that it was not a bad equation at all. The same equation describes continuous compound interest, radioactive decay, and many other processes.

    1. The Simplest Population Model

    Let P(t) denote a population at time t.

    One of the simplest assumptions we can make is this: the rate at which the population grows is proportional to the population itself.

    dP dt = kP, P(0) = P0.

    Here k>0 is a constant. The idea seems reasonable. If there are twice as many individuals, we might expect approximately twice as many births. A larger population therefore grows faster.

    The equation is separable:

    dPP = kdt.

    Integrating gives

    lnP = kt+C,

    and therefore

    P(t) = P0 ekt.

    This is exponential growth.

    There is an immediate problem. If k>0, then

    P(t) → ∞ as t→∞.

    According to this model, the population eventually becomes arbitrarily large. That cannot continue indefinitely in the real world. Food, water, space, and other resources are limited.

    So our first population model is useful for describing growth over some periods, but it cannot be the whole story.

    “`
    “`html

    2. Introducing a Carrying Capacity

    Suppose the environment can sustainably support a maximum population K. This number is called the carrying capacity.

    We modify our original equation to

    dP dt = kP ( 1 − P K )

    This is the logistic equation.

    The new factor

    ( 1 − P K )

    is what changes everything. When the population is small compared with K, this factor is close to 1, so the population behaves approximately like ordinary exponential growth. As the population becomes larger, the factor becomes smaller and the growth slows down.

    We Can Predict the Solutions Without Solving the Equation

    This is one of the most useful ideas in differential equations: we do not always need an explicit formula to understand what the solutions will do.

    First suppose

    0 < P < K.

    Then

    ( 1 − P K ) > 0,

    and consequently

    dP dt > 0.

    So the population increases.

    Now suppose P=K. Then

    ( 1 − P K ) = 0,

    so

    dP dt = 0.

    The population remains constant at the carrying capacity.

    Finally, if P>K, then

    ( 1 − P K ) < 0,

    and therefore

    dP dt < 0.

    The population decreases toward the carrying capacity.

    Where Is the Population Growing Fastest?

    The growth rate is

    kP ( 1 − P K ).

    As a function of P, this is a downward-opening quadratic:

    kP − kP2 K .

    Its maximum occurs at

    P = K 2 .

    This tells us something important about the shape of the population curve.

    If P0 < K2 , the population initially grows faster and faster. When it reaches P = K2 , its growth rate is greatest. After that, the population continues to increase, but more and more slowly as it approaches K. This produces the familiar S-shaped logistic curve.

    If K2 < P0 < K , the population begins above the point of fastest growth. It still increases toward K, but it slows down from the beginning.

    Finally, if P0 > K , the population decreases toward K.

    Thus, before solving the logistic equation, we can already predict the three different types of solution curves shown in the next figure.

    “`
    “`html id=”p87fcf”

    3. Now Let Us Solve the Logistic Equation

    We now return to the logistic equation

    dP dt = kP ( 1 − P K ), P(0) = P0.

    We have already learned a great deal about its solutions without solving it. Now let us find the actual formula.

    First separate the variables:

    dP P ( 1 − P K ) = kdt.

    Since

    1 P ( 1 − P K ) = 1P + 1 K−P ,

    we can integrate:

    ∫ ( 1P + 1 K−P ) dP = ∫ kdt.

    This gives

    lnP − ln ( K−P ) = kt + C.

    Combining the logarithms,

    ln ( P K−P ) = kt + C.

    Exponentiating both sides gives

    P K−P = C e kt .

    Using the initial condition P(0) = P0 , we obtain

    C = P0 K − P0 .

    After solving for P, we obtain the logistic growth formula:

    P(t) = K 1 + K − P0 P0 e −kt .

    Now the formula confirms what we predicted from the differential equation. For positive initial populations, the population approaches the carrying capacity:

    P(t) → K as t→∞.

    So the carrying capacity is not merely a number inserted into the model. It becomes the long-term population predicted by the model.

    4. Was Our Original Equation Really So Bad?

    We rejected the equation

    dy dt = ky

    as a model of population growth over an unlimited period of time. But the equation itself is one of the most important differential equations in mathematics.

    The initial-value problem

    dy dt = ky, y(0) = y0

    has the solution

    y(t) = y0 e kt .

    What changes from one application to another is the meaning of y and the sign and meaning of the proportionality constant.

    5. Continuous Compound Interest

    Suppose an amount of money A(t) earns interest continuously at an annual rate r. The rate at which the account balance changes is proportional to the amount currently in the account:

    dA dt = rA, A(0) = A0.

    Therefore,

    A(t) = A0 e rt .

    The same equation that produced exponential population growth now describes the growth of money.

    6. Radioactive Decay

    Now consider a radioactive substance. The more radioactive nuclei that are present, the more nuclei are available to decay. Thus, the magnitude of the decay rate is proportional to the amount currently present.

    This time the quantity is decreasing, so we write

    dN dt = − λN, N(0) = N0,

    where λ>0 is the decay constant.

    The solution is

    N(t) = N0 e − λt .

    7. One Equation, Many Processes

    We began with perhaps the simplest population model imaginable: the rate of change of a population is proportional to the population itself.

    As a long-term population model, it failed. Unlimited exponential population growth is impossible in an environment with limited resources.

    Introducing a carrying capacity led naturally to the logistic equation. Even more importantly, we were able to predict the behavior of its solutions before solving the equation.

    But our original equation was far from useless. The same basic mathematical law appears in continuous compound interest and radioactive decay.

    The common idea is simple:

    The rate of change of a quantity is proportional to the amount of that quantity currently present.

    Population, money, and radioactive atoms seem like completely different things. Mathematically, however, they can obey the same law.

    That is one of the remarkable features of differential equations: the same mathematical equation can describe very different processes in the real world.


    Related: See another example where a simple mathematical model produces a surprising—and ultimately unrealistic—prediction: A Sliding Ladder: Is It Better to Slide or Jump? .