# An Introduction to the Modelling of Real-World Problems by the Simplest Ordinary Differential Equations

## Stanisław Migórski

## 1 Introduction and preliminary terminology

This note is a preliminary version of a self-contained introduction to the modeling of selected real-life phenomena by simple ordinary di.erential equations.We restrict ourselves to a discussion on first order ordinary di.erential equations, leaving the detailson models described by higher order equations and systems of first-order equations to the sequel of the present work.

The text was designed for computer science students who have a calculus back ground1 and have taken prior physics courses (cf. [4]). It is our belief that computer science students should know how to model a selected problem, particularly in light of rapidly changing technologies. We begin with preliminary notations and defnitions.

By a dierential equation we mean a mathematical equation involving an unknown function of one or several variables and its derivatives of various orders. The order of the di.erential equation is the order of the highest derivative of the unknown func- tion involved in the equation. For instance, the equation

is called an ordinary differential equation (often abbreviated to ODE) of first order^{2}
in the normal form. It relates an independent variable t to an unknown function x and its first order derivative. The function f is given f: I × E → E, where I is an interval^{3}
and E is a prescribed space and we are looking for a function x: I → E. Of course the
equation (1) obliges the unknown function x to have some restrictions: it should be
differentiable in a suitable sense and the compose function f(t, x(t)) should have a
meaning. Classically, we admit the following

^{1}(I,E)

^{4}and x’(t) = f(t, x(t)) for all t ∈ I.

If the right hand side of the equation (1) can be expressed as a product g(t)h(x) where g
depends only on t and h depends only on x, then the differential equation (1) is called **separable**. If the right hand side f is independent of t, then the resulting differential
equation is called **autonomous**. Every autonomous differential equation is separable.

When the function f in (1) is (affine) linear with respect to the second variable, the first
order equation of the form x’ = α(t)x + β(t) with prescribed functions and is called
a nonhomogeneous linear differential equation. If β= 0, then the equation x=α(t)x
with a prescribed function, is called a **homogeneous linear** differential equation.

_{0}I and x

_{0}∈ E. A problem in which we are looking for the unknown function x: I → E of a differential equation (1) where the value of the unknown function at some point is known x(t

_{0}) = x

_{0}

^{5}is called an initial value problem (in short IVP) or a Cauchy

^{6}problem for (1).

Analogously as above, we have

_{0}∈ I and x

_{0}∈ E. We say that a function x: I → E is a solution of the IVP

if x is a solution to equation (1) on the interval I and it satisfies x(t_{0}) = x_{0}. Sometimes
if no initial condition is given, we call the family of all solutions to the differential
equation (1) the **general** solution.

## 2 Modelling

“How do we translate a physical phenomenon into a set of equations which describes it?” – this is certainly one of the most di.cult problems that scientists deal with intheir everyday research. This problem is a di.cult one since it is usually impossible to describe a phenomenon totally, so one often tries to reformulate a real-world problem as a mathematical one making certain simplyfying assumptions. As a result one usually describes the system approximately and adequately.

The problem of generating “good” equations is not an easy task. The set of
equations one deals with is called a **model** for the system. In general, once we have
built a set of equations, we compare the data generated by the equations with real data
collected from the system (by measurement). When the two sets of data “agree” (or are
“sufficiently” close), we gain con.dence that the set of equations will lead to a good description of the real-world system. If a prediction from the equations leads to some conclusions which are by no means close to real-world future behavior, we should modify
and “correct” the underlying equations. When creating a model, it is necessary to formulate
the problem under consideration into questions that can be answered mathematically.

The following are basic steps in building a model.

(1) Clearly state the assumptions on which the model will be based. These assumptions should describe the relationships between the quantities to be studied.

(2) Completely describe the parameters and variables to be used in the model.

(3) Use the assumptions (from step (1)) to derive mathematical equations relating the parameters and variables (from step (2)).

A large number of laws of physics, chemistry, economics, medicine, etc. can be for mulated as di.erential equations. They serve as models that describe the behavior of complex systems. The mathematical theory of differential equations first developed together with the sciences where the equations had originated and where the results found application. Whenever a mathematical model involves the rate of change of one variable with respect to another, a di.erential equation is apt to appear. As we will see in the following sections, diverse problems, sometimes originating in quite distinct scienti.c .elds, may give rise to identical di.erential equations. Whenever this happens, the mathematical theory behind the equations can be viewed as a unifying principle behind diverse phenomena.

In the following sections we provide examples of mathematical models of several rela- tively simple phenomena which are described by ordinary differential equations of the form (1). We mainly concentrate on modeling issues. The modern theory of ordinary differential equations together with known techniques, methods and applications can be found in references [1] and [2].

In the following sections we provide examples of mathematical models of several relatively simple phenomena which are described by ordinary differential equations of the form (1). We mainly concentrate on modeling issues. The modern theory of ordinary differential equations together with known techniques, methods and applications can be found in references [1] and [2].

## Newton’s law of cooling

*It is known from experimental observations that the surface temperature of an object
changes at a rate proportional to its relative temperature. The latter is the dierence
between its temperature and the temperature of the surrounding environment. This
is what is known as Newton’s law of cooling. Find the dependence of the temperature of the object on time.
*

We denote by T(t) the temperature of the object at time instant t and by M the
constant temperature^{7} of the environment. The Newton law of cooling gives a
differential equation of the first order

where k > 0. The proportionality constant in (13) is negative since the temperatures decrease in time and the derivative must be negative for T > M. The equation (13) is a separable ODE. We assume that the initial temperature of the object is prescribed

where T_{0} stands for the initial temperature of the object. Thus, the dependence of the
temperature of the object on time is described by the initial value problem (13) and
(14). Its solution has a form^{8}

_{0}= 34.5, M = 16, t = 60 and T(60) = 33.7. First, we find the value of the coefficient k as follows

In other words, we have determined k from the information on two temperatures
measured by Kowalski, at the point T_{0} = 34.5 and the point T(60) = 33.7.

Knowing the value of k, we are able to find a time instant td at which the crime happened.
To this end we take into account other two points T(0) = 37^{9} and T(t_{d}) = 34.5.

Therefore the murder occured about 172 minutes before noon, that is, around 9:08am.
We conclude that in order to find the time of death it is necessary to measure the
body temperature twice^{10}.

### References

[1] Z. Denkowski, S. Migórski and N.S. Papageorgiou, An Introduction to Nonlinear Analysis:
Theory. Kluwer Academic/Plenum Publishers, Boston, Dordrecht, London, New York (2003).

[2] Z. Denkowski, S. Migórski and N.S. Papageorgiou, An Introduction to Nonlinear Analysis:
Applications. Kluwer Academic/Plenum Publishers, Boston, Dordrecht, London, New York (2003).

^{1} Some universities, the Jagiellonian among them, make linear algebra a prerequisite for
differential equations. Many schools, especially engineering, only require calculus.

^{2} The most general (implicit) form of an ordinary differential equation of the first-order is as
follows: F(t, x, x’) = 0 with a prescribed function F.

^{3} The interval I may have a form [a, b], [a, b), (a, b], (a, b),-., b] and [a,+.).

^{4} A function x: I → E is C1, if x is dierentiable on I and its derivative x’ is a continuous function.

^{5} This condition is called an initial condition.

^{6} Augustin Louis Cauchy (1789–1857), a French mathematician.

^{7} The constant M is called the medium temperature since it is the temperature of the medium
an object is immersed in.

^{8} This can be easily checked by the reader.

^{9} Assuming the dead person was not sick and had a temperature of 37^{0}C. The interested
reader may estimate the time of death, if instead of 37^{0}C we suppose 36.6^{0}C.

^{10} Remember this method unless you are sergeant Colombo.

See also: