Ordinary Differential Equations: First-Order ODEs, Initial Value Problems, Separable and Linear Equations with Biological Applications

Ordinary Differential Equations: First-Order ODEs, Initial Value Problems, Separable & Linear Equations

Ordinary Differential Equations (ODEs) are one of the most useful mathematical tools for understanding how biological systems change with time. Population size, bacterial growth, drug concentration, enzyme activity, radioactive decay, disease spread, nutrient consumption and many other biological processes can be represented using differential equations.

The main idea is simple: instead of describing only the present value of a biological quantity, an ordinary differential equation describes its rate of change. Once we know how fast a quantity changes and provide suitable starting information, we can often predict its future behaviour.

In this chapter, we will develop the basic ideas of first-order ordinary differential equations in a step-by-step and exam-friendly manner. The discussion covers first-order differential equations, initial value problems, separable equations, linear differential equations and important applications to biological systems.

📚 Contents

1. Introduction to Ordinary Differential Equations
2. Basic Terms and Concepts
3. First-Order Differential Equations
4. Initial Value Problems
5. Separable Differential Equations
6. Linear Differential Equations
7. Applications to Biological Systems
8. Population Growth Models
9. Logistic Growth Model
10. Decay and Elimination Models
11. Drug Concentration Model
12. Simple Biological Reaction Models
13. Basic Disease and Infection Models
14. Important Comparison Table
15. Common Mistakes
16. Important Exam Facts
17. MCQs
18. Frequently Asked Questions
19. Last-Minute Revision
20. Summary

1. 📐 Introduction to Ordinary Differential Equations

A differential equation is an equation containing an unknown function and one or more of its derivatives. A derivative tells us how a quantity changes with respect to another variable.

For example, if N(t) represents the number of microorganisms in a culture at time t, then the derivative dN/dt represents the rate at which the population is changing.

dN/dt = rate of change of population

If the population increases rapidly when the population itself is large, a simple mathematical model is:

dN/dt = rN

Here, r is the intrinsic growth rate. This equation is a first-order ordinary differential equation because the highest derivative appearing in it is the first derivative.

Function A quantity whose behaviour we want to determine, such as N(t).

Independent Variable The variable with respect to which differentiation is performed, often time t.

Dependent Variable The quantity that changes in response to the independent variable.

What does "ordinary" mean?

An ordinary differential equation contains derivatives with respect to one independent variable. For example:

dy/dx = 3x + 2

Here, y depends only on x, so this is an ordinary differential equation.

A partial differential equation (PDE), in contrast, involves partial derivatives with respect to two or more independent variables. For example, temperature may depend on both position and time.

Feature	ODE	PDE
Independent variables	Usually one	Two or more
Derivative	Ordinary derivative	Partial derivative
Example	dy/dt = ky	∂u/∂t = D∂²u/∂x²
Common biological use	Population or concentration over time	Diffusion over space and time

2. 🔬 Basic Terms and Concepts

2.1 Order of a Differential Equation

The order of a differential equation is the order of the highest derivative present in the equation.

Equation	Order
dy/dx = x²	First order
d²y/dx² + y = 0	Second order
d³y/dx³ + dy/dx = 5	Third order

2.2 Degree

The degree is the power of the highest-order derivative after the equation has been expressed as a polynomial in derivatives.

For example:

(d²y/dx²)² + dy/dx + y = 0

The highest derivative is the second derivative and its power is 2. Therefore, the equation has order 2 and degree 2.

2.3 Solution of an ODE

A solution is a function that satisfies the differential equation over the relevant interval.

For example:

dy/dx = 2x

Integrating both sides gives:

y = x² + C

where C is an arbitrary constant.

💡 Remember: A first-order differential equation normally produces one arbitrary constant in its general solution.

3. 🧮 First-Order Differential Equations

A first-order ODE contains the first derivative of the unknown function but no higher derivative.

The general form can be written as:

dy/dx = f(x,y)

or sometimes:

F(x,y,dy/dx) = 0

First-order equations are especially important in biology because many biological processes can be approximated by describing the rate of change of one quantity with respect to time.

Examples

dN/dt = rN

dC/dt = -kC

dP/dt = rP(1 - P/K)

These equations can represent population growth, concentration decay and density-dependent population growth respectively.

4. 🎯 Initial Value Problems

An initial value problem (IVP) consists of a differential equation together with a condition specifying the value of the unknown function at a particular starting point.

For example:

dy/dt = ky,    y(0) = y₀

The differential equation describes the rate of change, while the condition y(0) = y₀ tells us where the solution starts.

Step 1: Solve the differential equation

Find the general solution containing an arbitrary constant.

Step 2: Apply the initial condition

Substitute the given starting value to determine the unknown constant.

Step 3: Obtain the particular solution

The resulting equation satisfies both the differential equation and the initial condition.

Example: Simple Growth IVP

Consider:

dN/dt = rN,    N(0) = N₀

The solution is:

N(t) = N₀e^rt

The initial population N₀ determines the particular growth curve.

📌 Biological interpretation: If N represents the number of bacteria, N₀ is the initial number of bacteria and r describes the intrinsic per-capita growth rate under the assumptions of the model.

5. 🔀 Separable Differential Equations

A differential equation is called separable when the variables can be rearranged so that all terms involving one variable are on one side and the other variable is on the other side.

A common form is:

dy/dx = g(x)h(y)

We rearrange it as:

1/h(y) dy = g(x) dx

Then integrate both sides.

General Method

Step 1 — Identify the variables

Determine whether the equation can be written as a product of a function of x and a function of y.

Step 2 — Separate them

Move all y-related terms with dy to one side and all x-related terms with dx to the other side.

Step 3 — Integrate

Integrate both sides independently.

Step 4 — Apply the initial condition

If an initial value is given, use it to determine the integration constant.

Example 1: Exponential Growth

Consider:

dN/dt = rN

Separate the variables:

dN/N = r dt

Integrate:

ln|N| = rt + C

Exponentiating gives:

N = Ae^rt

If N(0) = N₀, then A = N₀. Therefore:

N(t) = N₀e^rt

Example 2: Decay

Suppose a substance is removed at a rate proportional to its current concentration:

dC/dt = -kC

Separating and integrating:

dC/C = -k dt

ln|C| = -kt + C₁

Therefore:

C(t) = C₀e^-kt

💡 Key idea: The negative sign in a decay equation indicates that the quantity decreases with time when k is positive.

6. 📈 Linear Differential Equations

A first-order linear differential equation has the standard form:

dy/dx + P(x)y = Q(x)

The equation is linear in y because y and its derivative appear only to the first power and are not multiplied together or placed inside nonlinear functions.

Integrating Factor Method

The most important method for solving a first-order linear ODE is the integrating factor method.

For:

dy/dx + P(x)y = Q(x)

the integrating factor is:

IF = e^∫P(x)dx

Multiplying the entire differential equation by the integrating factor converts the left side into the derivative of a product.

IF × dy/dx + IF × P(x)y = IF × Q(x)

Then:

d/dx [IF × y] = IF × Q(x)

Integrate both sides and solve for y.

Step-by-Step Example

Consider:

dy/dx + 2y = 4

Step 1: Identify P(x)

P(x) = 2.

Step 2: Find the integrating factor

IF = e^∫2dx = e^2x

Step 3: Multiply the equation by e^2x

e^2xdy/dx + 2e^2xy = 4e^2x

Step 4: Recognize the product derivative

d/dx(y e^2x) = 4e^2x

Step 5: Integrate

ye^2x = 2e^2x + C

Step 6: Solve for y

y = 2 + Ce^-2x

7. 🧬 Applications of ODEs to Biological Systems

Differential equations are particularly useful in biology because biological systems are dynamic. The amount of a substance, number of cells or population of organisms can change continuously with time.

🦠 Microbial Growth ODEs can describe bacterial or microbial population increase.

🌱 Ecology Population interactions, competition and resource limitations can be modeled mathematically.

💊 Pharmacokinetics Drug absorption, distribution and elimination can be represented using concentration equations.

🧪 Biochemistry Reaction rates and concentration changes can be modeled using differential equations.

🦠 Epidemiology ODE models can describe changes in susceptible, infected and recovered populations.

🌿 Plant Biology Growth, nutrient uptake and biochemical processes can be represented by dynamic models.

8. 🦠 Population Growth Models

One of the simplest biological applications of differential equations is population growth.

Suppose N(t) represents the population size. If every individual has the same average contribution to population growth and resources are not limiting, we can write:

dN/dt = rN

This is called the exponential growth model.

Interpretation of r

The parameter r is the intrinsic per-capita growth rate. When r is positive, the population increases. When r is negative, the population decreases.

Value of r	Basic interpretation
r > 0	Population tends to increase
r = 0	No net change in the simple model
r < 0	Population tends to decrease

Solution

N(t) = N₀e^rt

This model predicts increasingly rapid growth when r is positive. In real biological systems, however, unlimited growth cannot continue indefinitely because food, space, nutrients and other resources become limiting.

⚠️ Important: Exponential growth is an idealized model. It is most appropriate over a period during which environmental limitations are weak enough that the assumptions are reasonable.

9. 🌱 Logistic Growth Model

The logistic growth model introduces environmental limitations into population growth.

The differential equation is:

dN/dt = rN(1 - N/K)

Here:

• N = population size
• r = intrinsic growth rate
• K = carrying capacity
• t = time

The carrying capacity represents the approximate population level that the environment can sustain under the model assumptions.

How the model behaves

When N is much smaller than K, the term:

1 - N/K

is close to 1. Therefore, growth is approximately exponential.

As N approaches K, the factor becomes smaller and population growth slows.

When N = K:

dN/dt = 0

Thus K is an equilibrium population in the basic logistic model.

10. ☢️ Decay and Elimination Models

The same mathematical framework can describe processes in which a biological quantity decreases with time.

If the rate of decrease is proportional to the amount currently present:

dC/dt = -kC

The solution is:

C(t) = C₀e^-kt

This type of equation can describe idealized models of drug elimination, radioactive decay and loss of a substance from a compartment.

Half-Life

For first-order exponential decay, the half-life is:

t_1/2 = ln(2)/k

This means that the time required for the quantity to fall to half its initial value is constant in the simple first-order decay model.

📌 Exam point: First-order decay produces an exponential curve. The concentration decreases quickly at first and then approaches zero gradually.

11. 💊 Drug Concentration Model

Differential equations are widely used in pharmacokinetics to describe how drug concentration changes over time.

Consider a simplified one-compartment elimination model:

dC/dt = -kC

where C is drug concentration and k is the elimination rate constant.

The solution is:

C(t) = C₀e^-kt

This is a simplified model. Real pharmacokinetic systems can involve multiple compartments, absorption, distribution, metabolism and excretion, but the basic first-order equation provides an important mathematical foundation.

Continuous Drug Input

Suppose a drug enters a compartment at a constant rate R and leaves proportionally to its concentration. A simple model is:

dC/dt = R - kC

This is a first-order linear differential equation.

At equilibrium, the concentration stops changing, so dC/dt = 0:

0 = R - kC

Therefore:

C = R/k

This illustrates an important biological idea: a constant input combined with proportional removal can lead to a stable steady-state concentration.

12. 🧪 Simple Biological Reaction Models

Differential equations can also represent biochemical reactions. For a simple conversion of a substrate S into a product P:

S → P

If the reaction rate is proportional to substrate concentration:

dS/dt = -kS

Then substrate concentration follows exponential decay:

S(t) = S₀e^-kt

The negative sign indicates that substrate is being consumed.

The corresponding product formation rate in a simplified closed system can be related to substrate loss.

💡 Biological connection: In real enzyme kinetics, reaction rates may not remain proportional to substrate concentration over all concentration ranges. More detailed models, such as the Michaelis-Menten framework, may be required.

13. 🦠 Basic Disease and Infection Models

ODEs are also used in epidemiology. A classic conceptual framework divides a population into compartments such as:

S — Susceptible Individuals who can potentially become infected.

I — Infected Individuals currently infected or infectious according to the model definition.

R — Recovered/Removed Individuals moved out of the infectious compartment according to the model assumptions.

In a basic SIR model, a common set of equations is:

dS/dt = -βSI/N

dI/dt = βSI/N - γI

dR/dt = γI

Here β represents an effective transmission parameter and γ represents the rate of removal from the infectious compartment. N represents the total population under the model assumptions.

⚠️ Important: Mathematical disease models are simplified representations. Real epidemics can involve age structure, spatial movement, immunity, vaccination, reporting differences, multiple transmission routes and many other factors.

14. 📊 Important Comparison Table

Type	General Form	Main Method	Example
First-order ODE	dy/dx = f(x,y)	Depends on form	dy/dt = rN
Separable	dy/dx = g(x)h(y)	Separate and integrate	dN/dt = rN
Linear	dy/dx + P(x)y = Q(x)	Integrating factor	dC/dt + kC = R
Decay	dC/dt = -kC	Separation	Drug elimination
Exponential growth	dN/dt = rN	Separation	Idealized microbial growth
Logistic growth	dN/dt = rN(1-N/K)	Separation or standard logistic solution	Density-dependent growth

15. ⚠️ Common Mistakes While Solving ODEs

1. Forgetting the integration constant: Indefinite integration requires an arbitrary constant.
2. Incorrect separation: Make sure dy remains with the y terms and dx with the x terms.
3. Ignoring the initial condition: An IVP requires the initial value to determine the particular solution.
4. Wrong sign in decay: A simple decay model normally has a negative rate term.
5. Using the wrong integrating factor: For y' + P(x)y = Q(x), use e^∫P(x)dx.
6. Confusing order and degree: Order refers to the highest derivative; degree refers to its power after appropriate algebraic arrangement.
7. Ignoring biological assumptions: A mathematical solution is meaningful only in the context of the assumptions of the model.
8. Assuming exponential growth continues forever: Real populations usually encounter environmental limitations.

16. 📝 Important Exam Facts

Fact 1A first-order ODE contains a first derivative as its highest derivative.

Fact 2An IVP combines a differential equation with an initial condition.

Fact 3Separable equations are solved by separating variables and integrating.

Fact 4The integrating factor for y' + P(x)y = Q(x) is e^∫P(x)dx.

Fact 5dN/dt = rN represents exponential population growth.

Fact 6Logistic growth includes the carrying capacity K.

Fact 7For exponential growth, N(t) = N₀e^rt.

Fact 8For first-order decay, C(t) = C₀e^-kt.

Fact 9For first-order decay, half-life is ln(2)/k.

Fact 10ODEs are widely used in population biology and pharmacokinetics.

Fact 11At N = K in the basic logistic model, dN/dt = 0.

Fact 12Negative growth parameters can represent decline in simple population models.

17. 🎯 Multiple Choice Questions

Q1. What is the order of the equation dy/dx = x + y?

A. Zero

B. First

C. Second

D. Third

Answer: B. First

Q2. Which condition makes an ODE an initial value problem?

A. A boundary condition only

B. An initial value for the unknown function

C. A second derivative

D. A numerical coefficient

Answer: B. An initial value for the unknown function

Q3. Which equation is separable?

A. dy/dx = xy

B. dy/dx + y = x

C. y'' + y = 0

D. x + y = 1

Answer: A. dy/dx = xy

Q4. What is the integrating factor for y' + P(x)y = Q(x)?

A. P(x)

B. Q(x)

C. e^∫P(x)dx

D. e^Q(x)

Answer: C. e^∫P(x)dx

Q5. Which equation represents simple exponential population growth?

A. dN/dt = r + N

B. dN/dt = rN

C. dN/dt = rN(1-N/K)

D. dN/dt = -rN

Answer: B. dN/dt = rN

Q6. What does K represent in the logistic growth model?

A. Mutation rate

B. Carrying capacity

C. Decay constant

D. Initial population

Answer: B. Carrying capacity

Q7. Which equation represents first-order exponential decay?

A. dC/dt = kC

B. dC/dt = -kC

C. dC/dt = C²

D. dC/dt = k + C

Answer: B. dC/dt = -kC

Q8. The solution of dN/dt = rN with N(0) = N₀ is:

A. N = N₀ + rt

B. N = N₀e^rt

C. N = rt/N₀

D. N = N₀e^-rt

Answer: B. N = N₀e^rt

Q9. In the simple decay model, what does k usually represent?

A. Carrying capacity

B. Decay or elimination rate constant

C. Initial concentration

D. Population size

Answer: B. Decay or elimination rate constant

Q10. Which biological process is commonly modeled using ODEs?

A. Population growth

B. Drug concentration changes

C. Disease transmission models

D. All of the above

Answer: D. All of the above

Q11. In the logistic equation, what happens to the growth rate when N approaches K?

A. It increases without limit

B. It approaches zero

C. It becomes equal to K

D. It always becomes negative

Answer: B. It approaches zero

Q12. Which method is particularly useful for first-order linear ODEs?

A. Integrating factor

B. Matrix inversion only

C. Fourier transform only

D. Factorial expansion

Answer: A. Integrating factor

Q13. What does dN/dt describe biologically?

A. Total time

B. Rate of change of population

C. Carrying capacity

D. Initial population only

Answer: B. Rate of change of population

Q14. For dC/dt = -kC, k > 0, the concentration generally:

A. Increases exponentially

B. Decreases exponentially

C. Remains constant

D. Becomes exactly zero immediately

Answer: B. Decreases exponentially

Q15. In the simple logistic model, dN/dt becomes zero when:

A. N = 0 or N = K

B. N = 1 only

C. r = K always

D. t = 1 only

Answer: A. N = 0 or N = K

18. ❓ Frequently Asked Questions

What is an ordinary differential equation?

An ordinary differential equation is an equation involving an unknown function and its ordinary derivatives with respect to one independent variable.

What is a first-order ODE?

A first-order ODE is an ODE in which the highest derivative is the first derivative.

What is an initial value problem?

An initial value problem contains a differential equation together with a specified value of the unknown function at a particular starting point.

What is a separable equation?

A separable equation is one in which variables can be rearranged so that the x-dependent and y-dependent terms occur on opposite sides of the equation.

What is an integrating factor?

For a first-order linear equation y' + P(x)y = Q(x), the integrating factor is e^∫P(x)dx. It converts the left-hand side into the derivative of a product.

Why are ODEs important in biology?

Biological systems are dynamic. ODEs allow researchers to describe and predict changes in populations, concentrations, drug levels, infection numbers and other quantities over time.

What is the difference between exponential and logistic growth?

Exponential growth assumes that the per-capita growth rate remains constant, whereas logistic growth introduces a carrying capacity and therefore includes density-dependent limitation.

What does the negative sign mean in dC/dt = -kC?

It indicates that the quantity decreases with time when k is positive.

19. ⚡ Last-Minute Revision Table

Concept	Remember This
ODE	Differential equation involving one independent variable.
First order	Highest derivative is first derivative.
IVP	ODE + initial condition.
Separable ODE	Separate variables and integrate.
Linear ODE	y' + P(x)y = Q(x).
Integrating factor	IF = e∫P(x)dx.
Exponential growth	dN/dt = rN.
Exponential solution	N = N₀ert.
Decay	dC/dt = -kC.
Decay solution	C = C₀e-kt.
Half-life	t1/2 = ln(2)/k for first-order decay.
Logistic growth	dN/dt = rN(1-N/K).
K	Carrying capacity.
Biological use	Population, pharmacokinetics, biochemistry and epidemiology.

20. 🧠 Easy Memory Tricks

SEPARATE → INTEGRATE → APPLY
For a separable ODE: separate the variables, integrate both sides, then apply the initial condition if one is provided.

LINEAR → PIF
For a first-order linear ODE, remember: P(x) → calculate the integrating factor → integrate → solve.

GROWTH = +rN
Positive proportional rate gives exponential growth.

DECAY = -kC
The negative sign represents loss or elimination in the basic first-order decay model.

LOGISTIC = GROWTH + LIMIT
The term (1 - N/K) introduces environmental limitation into the simple exponential model.

21. 📌 Chapter Summary

Ordinary differential equations provide a mathematical language for describing changing biological systems. A first-order ODE focuses on the first derivative, which commonly represents the rate at which a biological quantity changes.

An initial value problem adds starting information to a differential equation, allowing a general solution to be converted into a particular solution. Separable differential equations can be solved by rearranging the variables and integrating both sides.

A first-order linear differential equation has the form:

dy/dx + P(x)y = Q(x)

and is commonly solved using the integrating factor:

IF = e^∫P(x)dx

Biology provides many practical examples of differential equations. Exponential growth can describe an idealized population under suitable assumptions, while logistic growth incorporates a carrying capacity. First-order decay can represent simplified elimination processes. ODEs also provide the foundation for mathematical models in pharmacokinetics, biochemical reactions, epidemiology, ecology and systems biology.

The most important lesson is that an ODE does not simply produce a mathematical formula. It connects a rate of change with a biological mechanism. The quality of a biological model therefore depends not only on solving the equation correctly but also on understanding the assumptions behind the model.

📐 Mathematics Understand derivatives, integration and differential equations.

🧬 Biology Translate biological processes into rate equations.

🔬 Modeling Use the mathematical solution to understand or predict system behaviour.

⭐ Final Exam Tip: Before solving any first-order ODE, first identify its form. Ask: Is it separable? Is it linear? Is an initial condition given? Choosing the correct method is often the most important step.

⬆ Back to Top