Chapter 2 – First-Order Differential Equations

2.1 Autonomous Differential Equation

2.1.1 Autonomous First-Order Differential Equation

In the autonomous DE, the independent variable does not appear in the equation explicitly. The general and normal forms of such equations are:

(1)   \begin{equation*}F(y, y^\prime) = 0, \quad \quad  \frac{d y}{d x} = f(y).\end{equation*}


For the definitions of the general and normal forms of a differential equation, see the definitions in Chapter 1.

2.1.2 Critical Points of an Autonomous DE

The zeros of function f(y) in the above definition, are known as the critical points of an autonomous differential equation. These points are also known as equilibrium points or stationary points.

2.1.3 Attractors and Repellers of an Autonomous Differential Equation

The order of a differential equation is the order of the highest derivative in the equation. A second-order PDE and a third-order ODE are shown below.
This is explained in details in the following video.

2.2 Separable Differential Equations

2.2.1 Separable First-Order Differential Equation

A separable differential equation has a general form of

(2)   \begin{equation*}\frac{d y}{d x} = g(x) h(y).\end{equation*}


This type of differential equations can be solved by separating the variables as follows

(3)   \begin{equation*}\frac{d y}{h(y)} = g(x) dx \quad \Rightarrow \quad \int \frac{d y}{h(y)} = \int g(x) dx\end{equation*}

2.3 Linear Differential Equations

2.3.1 Linear First-Order Differential Equation

A linear differential equation in the dependent variable y has a general form of

(4)   \begin{equation*}a_1 (x) \frac{dy}{dx} + a_0 (x) y = g(x)\end{equation*}

2.3.2 Homogeneous and Non-Homogeneous Differential Equations

In the above differential equation, if g(x) = 0 the linear differential equation is said to be homogeneous, and if g(x) \neq 0 is said to be non-homogeneous.

2.3.3 The Standard Form of a Linear Differential Equation

A linear differential equation in the standard form will be as follows

(5)   \begin{equation*}\frac{dy}{dx} + P(x) y = f(x)\end{equation*}

2.3.4 The Solution of a Linear Differential Equation

The first step in solving such differential equations is to write it in the standard form. Then the integrating factor can be obtained as

(6)   \begin{equation*}I(x) = e^{\int P(x) dx}.\end{equation*}


The integration factor must be multiplied to both the left and right hand-side of the differential equation, which yields to

(7)   \begin{equation*}\frac{d}{dx} \left(I(x) y \right) = I(x) f(x).\end{equation*}


By integrating both sides of this equation, the solution to the linear differential equation will be obtained.

(8)   \begin{equation*}I(x) y = \int I(x) f(x) dx.\end{equation*}

2.4 Exact Equations

2.4.1 The Definition of an Exact Equation

A first-order differential equation of the following form is said to be exact, if the left hand-side of the expression is an exact differential.

(9)   \begin{equation*}M(x,y) dx + N(x,y) dy = 0.\end{equation*}

2.4.2 What is an exact differential?

The necessary and sufficient condition for the above differential equation to be an exact differential is

(10)   \begin{equation*}\frac{\partial M}{\partial y} = \frac{\partial N}{\partial x},\end{equation*}


provided that M(x,y) and N(x,y) are continuous and have continuous first partial derivatives in a rectangular region.

2.4.3 How to solve an exact equation?

For these type of equations, a function f exists for which

(11)   \begin{equation*}\frac{\partial f}{\partial x} = M(x,y).\end{equation*}


f can be found by the following integration

(12)   \begin{equation*}f(x,y) = \int M(x,y) dx + g(y),\end{equation*}


where g(y) is an arbitrary function of y and a constant, with respect to x, of integration. Then the differentiation of f with respect to y will be

(13)   \begin{equation*}\frac{\partial f}{\partial y} = \frac{\partial}{\partial y} \left( \int M(x,y) dx + g(y) \right) = \frac{\partial}{\partial y} \left( \int M(x,y) dx \right) + g^\prime (y) .\end{equation*}


We assume that \partial f / \partial y = N(x,y), then

(14)   \begin{equation*}g^\prime (y) = N(x,y) - \frac{\partial}{\partial y} \int M(x,y) dx .\end{equation*}


By integrating this equation with respect to y, we can find g(y) and the problem will be solved.

2.5 Solution by Substitution

2.5.1 When to Use this Method?

A differential equation is usually solved by recognizing it as a certain kind of equation. Often the first step in solving a given differential equation consists of transforming it into another differential equation by means of substitution.

2.5.2 Homogeneous Function

If a function f has a property of

(15)   \begin{equation*}f(tx,ty) = t^\alpha f(x,y)\end{equation*}


for some real number \alpha, then f is said to be a homogeneous function of degree \alpha.

2.5.3 Homogeneous Differential Equation

The following differential equation is said to be homogeneous if both coefficients M and N are homogeneous functions of the same degree.

(16)   \begin{equation*}M(x,y) dx + N(x,y) dy = 0,\end{equation*}


Either of the substitutions y=ux or x=vy will reduce a homogeneous differential equation to a separable first-order differential equation.

2.5.4 Bernoulli’s Equation

The following differential equation is known as Bernoulli’s equation

(17)   \begin{equation*}\frac{dy}{dx} + P(x) y = f(x) y^n,\end{equation*}


where n is any real number. This differential equation is reduced to a linear equation, using the substitution u = y^{1-n}, and will be solved as a linear differential equation.

2.5.5 Reduction to Separation of Variables

A differential equation of the form

(18)   \begin{equation*}\frac{dy}{dx} = f(Ax+By+C)\end{equation*}


can always be reduced to an equation with separable variables by means of the substitution u=Ax+By+C, where B \neq 0.

2.6 Linear Models

2.6.1 Growth and Decay

If we know the population at time t_0, then we can find the population at any time t, where t > t_0, by solving the following initial-value differential equation.

(19)   \begin{equation*}\frac{dx}{dt} = kx, \quad x(t_0) = x_0,\end{equation*}


where x is the population and k is called the constant of proportionality. This DE is a first-order linear separable differential equation, and can be solved by the method of seperable equations and/or the method introduced for the first-order linear differential equations.

2.6.2 Solution to the Growth and Decay Method as a Separable Equation

To solve this equation, as a separable DE, we need to first separate the variables,

(20)   \begin{equation*}\frac{dx}{x} = k dt,\end{equation*}


then we need to integrate it, at find the constant of integration by using the initial condition.

(21)   \begin{equation*}\int \frac{dx}{x} = \int k dt \quad \Rightarrow \quad \ln \left| x \right| = k t + c \quad \Rightarrow \quad x = e^{k t + c} = e^{kt} e^c = C e^{kt}\end{equation*}


(22)   \begin{equation*}x(t=0) = x_0 \quad \rightarrow \quad x_0 = C e^{k (0)} = C e^0 = C \quad \Rightarrow \quad C = x_0\end{equation*}


So, the population at any given time, t, will be

(23)   \begin{equation*}x(t) = x_0 e^{kt} .\end{equation*}

2.6.3 Newton’s Law of Cooling/Warming

The mathematical formulation for cooling or warming a substance is given as a linear first-order differential equation. This equation is

(24)   \begin{equation*}\frac{dT}{dt} = k \left( T - T_m \right),\end{equation*}


where T is the temperature of the substance at time t, and T_m is the ambient temperature.

2.7 Problems

  1. Find the critical points and phase portrait of the given autonomous first-order differential equation. Classify each critical point as asymptotically stable, unstable, or semi-stable.

    (a) \quad \quad \frac{dy}{dx} = y^2 - 3y
    Answer: y = 0 is an attractor and y = 3 is a repeller. Watch the solution here.

    (b) \quad \quad \frac{dy}{dx} = y \ln(y+2)
    Answer: y = -1 is an attractor and y = 0 is a repeller. Watch the solution here.

    (c) \quad \quad \frac{dy}{dx} = \frac{y e^y - 9y}{e^y}
    Answer: y = 0 is an attractor and y = \ln 9 is a repeller. Watch the solution here.

  2. Solve the given differential equations by separation of variables.

    (a2.2.7) \quad \quad \frac{dy}{dx} = e^{3x+2y}

    (b2.2.9) \quad \quad y \ln x \frac{dx}{dy} = \left( \frac{y+1}{x} \right)^2

  3. Find an implicit and an explicit solution to the given initial value problem.

    (a2.2.25) \quad \quad x^2 \frac{dy}{dx} = y - xy, \quad y(-1) = -1

  4. Find the general solution of the given differential equation.

    (a) \quad \quad \frac{dy}{dx} + 2y = 0
    Answer: y = c e^{-2x}. Watch the solution here.

    (b) \quad \quad  \left( x^2 - 1 \right) \frac{dy}{dx} + 2y = \left( x+1 \right)^2
    Answer: y = \frac{(x+c)(x+1)}{(x-1)}. Watch the solution here.

    (c2.3.12) \quad \quad (1+x) \frac{dy}{dx} - xy = x + x^2

    (d2.3.18) \quad \quad \cos^2 x \sin x \frac{dy}{dx} + \left( \cos^3 x \right) y = 1

    (e2.5.2) \quad \quad (x+y) dx + x dy = 0

    (f2.5.7) \quad \quad \frac{dy}{dx} = \frac{y-x}{y+x}

    (g2.5.27) \quad \quad \frac{dy}{dx} = 2 + \sqrt{y - 2x + 3}

    (h2.5.23) \quad \quad \frac{dy}{dx} = \left( x + y + 1 \right)^2

  5. Determine whether the given differential equation is exact. If it is exact, solve it.

    (a) \quad \quad (2x-1) dx + (3y+7)dy = 0.
    Answer: It is exact, and the solution is x^2 - x + \frac{3}{2} y^2 + 7y = c. Watch the solution here.

    (b2.4.16) \quad \quad (5y-2x) y^\prime - 2y = 0

  6. Solve the given initial value problem.

    (a2.4.21) \quad \quad (x+y)^2 dx + (2xy + x^2 - 1) dy = 0, \quad y(1) = 1

    (b2.4.38) \quad \quad (x^2 + y^2 - 5) dx = (y+xy) dy , \quad y(0) = 1

    (c2.5.14) \quad \quad y dx + x ( \ln x - \ln y - 1) dy = 0, \quad y(1) = e
Reference

Dennis G. Zill. Advanced Engineering Mathematics, 6^{th} edition. Jones & Bartlett Learning. 2016.