1.1. Basic Concepts. Modeling
1.1. Basic Concepts. Modeling
First-Order ODEs
2 2
∂ u ∂ u
2
+ 2 =0
∂x ∂ y
PDEs have important engineering applications, but they are more complicated
than ODEs; they will be considered in Chap. 12.
( 4 ) F ( x , y , y ' ) =0
'
y =f ( x , y ) .
This is called the explicit form, in contrast to the implicit form (4). For
instance, the implicit ODE x−3 y ' −4 y 2 =0 (where x ≠ 0 ) can be written
explicitly as y ' =4 x 3 y 2.
Concept of Solution
A function
y=h ( x )
Is called a solution of a given ODE (4) on some open interval a< x<b if h ( x )
is defined and differentiable throughout the interval and is such that the
equation becomes an identity if y and y ' are replaced with h and h' ,
respectively. The curve (the graph) of h is called a solution curve.
Here, open interval a< x<b means that the endpoints a and b are not
regarded as points belonging to the interval. Also, a< x<b includes infinite
intervals −∞< x <b , a< x< ∞,−∞< x< ∞ (the real line) as special cases.
' dy
The ODE y = =cos x can be solved directly by integration on both sides.
dx
Indeed, using calculus, we obtain y=∫ cos x dx=sin x +c , where c is an
arbitrary constant. This is a family of solutions. Each value of c , for instance,
2.75 or 0 or −8 , gives one of these curves. Figure 3 shows some of them, for
c=−3 ,−3 ,−1,0,1,2,3,4.
Fig. 3. Solutions y=sin x +c of the ODE y ' =cos x
dy
y'= =0.2 e 0.2t=0.2 y
dt
Hence y is a solution of y ' =0.2 y (Fig. 4A). This ODE is of the form y ' =ky .
With positive-constant k it can model exponential growth, for instance, of
colonies of bacteria or populations of animals. It also applies to humans for
small populations in a large country (e.g., the United States in early times) and
is then known as Malthus’s law. We shall say more about this topic in Sec.
1.5.
(B) Similiarly, y ' =−0.2 (with a minus on the right) has the solution
−0.2 t
y=c e , (Fig. 4B) modeling exponential decay, as, for instance, of a
radioactive substance (see Example 5).
We see that each ODE in these examples has a solution that contains an
arbitrary constant c . Such a solution containing an arbitrary constant c is called a
general solution of the ODE.
(We shall see that c is sometimes not completely arbitrary but must be
restricted to some interval to avoid complex expressions in the solution.)
We shall develop methods that will give general solutions uniquely (perhaps
except for notation). Hence we shall say the general solution of a given ODE
(instead of a general solution).
In most cases, general solutions exist, and every solution not containing an
arbitrary constant is obtained as a particular solution by assigning a suitable value
to c . Exceptions to these rules occur but are of minor interest in applications; see
Prob. 16 in Problem Set 1.1.
( 5 ) , y ' =f ( x , y ) , y ( x 0 ) = y 0
dy
y'= =3 y , y ( 0 ) =5.7
dx
More on Modeling
The general importance of modeling to the engineer and physicist was emphasized
at the beginning of this section. We shall now consider a basic physical problem
that will show the details of the typical steps of modeling. Step 1: the transition
from the physical situation (the physical system) to its mathematical formulation
(its mathematical model); Step 2: the solution by a mathematical method; and
Step 3: the physical interpretation of the result. This may be the easiest way to
obtain a first idea of the nature and purpose of differential equations and their
applications. Realize at the computer (your CAS) may perhaps give you a hand in
Step 2, but Steps 1 and 3 are basically your work. And Step 2 equires a solid
knowledge and good understanding of solution methods available to you-you have
to choose the method for your work by hand or by the computer. Keep this in
mind, and always check computer results for errors (which may arise, for
instance, from false inputs).
Given an amount of a radioactive substance, say, 0.5g (gram), find the amount
present at any later time.
dy
(6 ), =−ky
dt
Where the constant k is positive, so that, because of the minus, we do get decay
(as in [B] of Example 3).
Now the given initial amount is 0.5 g, and we can call the corresponding
instant t=0 . Then we have the initial condition y ( 0 )=0.5. This is the instant at
which our observation of the process begins. It motivates the term initial
condition (which, however, is also used when the independent variable is not time
or when we choose a t other than t=0 ). Hence the mathematical model of the
physical process is the initial value problem
dy
(7) , =−ky , y ( 0 )=0.5
dt
( 8 ) , y (t )=c e−kt
We now determine c by using the initial condition. Since y ( 0 )=c from (8), this
gives y ( 0 )=c=0.5. Hence the particular solution governing our process is (cf.
Fig. 5)
( 9 ) , y (t )=0.5−kt , ( k >0 )
Always check your result-it may involve human or computer errors! Verify by
differentiation (chain rule!) that your solution (9) satisfies (7) as well as y ( 0 )=0.5:
dy −kt −kt 0
=−0.5 k e =−k ∙ 0.5 e =−ky , y ( 0 )=0.5 e =0.5 .
dt
Formula (9) gives the amount of radioactive substance at time t . It starts from the
correct initial amount and decreases with time because k is positive. The limit of
y as t → ∞ is zero.