In this post, we show how to use math to study population dynamics. The use of mathematical models is useful to study many phenomena that exist in life such as population dynamics. Using an easy way, we explain how to represent the density of a population by means of age.

As a matter of fact, you should know the basic notion of differential functions and integrals, as we shall work on continuous time.

### Math model of population dynamics by means of age

In mathematics, we characterize the evolution of any phenomenon by a function of tow variables $p(t,x)$. Here $t$ is a positive real number describing the time, and $x$ is a vector, spacial variable, belonging to a vector space $\mathbb{R}^n$.

In this post, we shall denote by $p(t,a)$ the density of a population of age $a$ at time $t$. Here $t\ge 0$ and $a\ge 0$ are time of the evolution and the age $a\ge 0$ (which is a time its self). In fact, $p(t,a)$ represent the number of individual of age $a$ at time $t$. We remark that $t$ and $a$ follow the same evolution. This means that the derivative with respect to $t$ is the same as the derivative with respect to $a$; hence

\begin{align*}\frac{d}{dt}p(t,a)=\frac{d}{da}p(t,a).\end{align*}

On the other hand, we should not count died people of age $a$ at time $t$. The number of died people of age $a$ depends of a coefficient of mortality rate caused by natural death, denoted by $\mu(a)$; so we shall deduce the number of mortality $\mu(a)p(t,a)$ from the above equation. Thus the equation becomes

\begin{align*}\frac{d}{dt}p(t,a)=\frac{d}{da}p(t,a)-\mu(a)p(t,a).\end{align*}

Moreover, we notice that at initial times $t=0,$ the population has already a density represented by a function $f(a)$ that give the number of each age $a$. The initial density is $p(0,a)=f(a)$. We get an other equation

\begin{align*}\begin{cases}\frac{d}{dt}p(t,a)=\frac{d}{da}p(t,a)-\mu(a)p(t,a)\\ p(0,a)=f(a)\end{cases}\end{align*}

### Equation of the new born

Until now we only discussed the density of existent individual in a population. But there is still a class of individuals missing that we have not yet added to the model. Indeed, these are the new born, individuals of age $a=0$, just start a life. We mention that the number of new born, which is $p(t,0)$, heavily depends on a coefficient of fertility of the population. This latter denoted by $\beta(a)$ and measure the fertility of persons of age $a$. Using this, determining the number of new born is to calculate the sum of $\beta(a)p(t,a)$ from $a=0$ to $a=\infty$. More precisely,

\begin{align*}p(t,0)=\int^{+\infty}_0 \beta(a)p(t,a)da.\end{align*}

Now the global modal the density of a population dynamic with respect to the age satisfies

\begin{align*}\begin{cases}\frac{d}{dt}p(t,a)=\frac{d}{da}p(t,a)-\mu(a)p(t,a)\\ p(0,a)=f(a),\\ p(t,0)=\int^{+\infty}_0 \beta(a)p(t,a)da. \end{cases}\end{align*}

This is a math model that allow to completly study population dynamics.

### How to solve the model of population dynamics

Although, the solvability of the aforementioned model is not a standard problem. This depends on knowledge of some theories like operator theory; semigroup theory; and spectral theory. In fact, the first point is to determine the space of work; which is the space of solutions. As we deal with the density of continuous ages, the ideal space is $L^1([0,+\infty))$. That is the space of integrable functions. Moreover, as we work with differential operators we shall also define the so-called Sobolev spaces. They are spaces of the form

\begin{align*} W^{1,1}([0,+\infty)):=\left\{f\;\text{absolutly continuous}, \;f’\in L^1\right\}.\end{align*}

On the other hand, we define these operators

\begin{align*} Af=f’-\mu f,\quad D(A)=\left\{f\in W^{1,1}:f(0)=\int^{+\infty}_0 \beta(a)f(a)da\right\}.\end{align*}

One can reformulate the population model as a standard Cauchy problem in $L^1$ as follow:

\begin{align*} \dot{p}(t)=Ap(t),\quad p(0)=f.\end{align*}

For instance, If the operator $A$ satisfies the Hille-Yosida or Lumer-Phillips theorems, the Cauchy problem has a unique solution.