Calculus: Differential Equations Introduction

Differential equations are the mathematical language of change. They describe how quantities evolve over time, space, or any other variable, making them indispensable for modeling real-world phenomena from physics and engineering to biology and economics. Mastering first-order differential equations provides you with the foundational toolkit to translate a verbal description of a dynamic system into a solvable mathematical model and predict its future behavior.

What is a Differential Equation?

A differential equation is an equation that relates a function to its derivatives. In essence, it describes a rule for how a quantity changes. When the function depends on only one independent variable (often time $t$), the equation is called an ordinary differential equation (ODE). The order of an ODE is the order of the highest derivative present. This article focuses on first-order ODEs, which involve only the first derivative $dy/dt$ or $y^{\prime }$.

The general form is often written as $\frac{dy}{dt}=f(t,y)$. The solution is not a single number, but a function $y(t)$ that satisfies the equation. An initial condition, like $y(t_0)=y_0$, specifies the state of the system at a starting point, allowing you to find a unique, particular solution from the infinite family of possible solutions.

Solving by Separation of Variables

The first and often simplest technique is solving separable equations. An ODE is separable if you can algebraically manipulate it to isolate all terms involving the dependent variable $y$ with $dy$ on one side and all terms involving the independent variable $t$ with $dt$ on the other.

The standard form is $\frac{dy}{dt}=g(t)h(y)$. The solution process is methodical:

1. Separate: Rewrite as $\frac{1}{h(y)}dy=g(t)dt$.
2. Integrate: Integrate both sides: $\int \frac{1}{h(y)}dy=\int g(t)dt$.
3. Solve: Solve the resulting equation for $y$ as a function of $t$, if possible. Include the constant of integration, $C$, from one of the integrals.

Example - Radioactive Decay: The rate of decay of a radioactive substance is proportional to the amount present. This translates to $\frac{dA}{dt}=-kA$, where $A(t)$ is the amount at time $t$ and $k>0$ is the decay constant.

1. Separate: $\frac{1}{A}dA=-kdt$.
2. Integrate: $\int \frac{1}{A}dA=\int -kdt\Rightarrow \ln |A|=-kt+C$.
3. Solve: $A=e^{-kt+C}=e^C{e}^{-kt}$. Letting $A_0=e^C$ be the initial amount, the solution is $A(t)=A_0{e}^{-kt}$. This exponential decay model is a core application.

Solving Linear Equations with an Integrating Factor

Not all first-order ODEs are separable. A crucial solvable type is the first-order linear differential equation, which can be written in the standard form $\frac{dy}{dt}+P(t)y=Q(t)$. The powerful technique for solving these involves the integrating factor.

The integrating factor, $\mu (t)$, is defined as $\mu (t)=e^{\int P(t)dt}$. The genius of this method is that multiplying the entire standard form equation by $\mu (t)$ makes the left side a perfect derivative: $\frac{d}{dt}[\mu (t)y]=\mu (t)Q(t)$. You then solve by integrating both sides with respect to $t$.

Step-by-Step Process:

1. Put the equation in standard form: $y^{\prime }+P(t)y=Q(t)$.
2. Compute the integrating factor: $\mu (t)=e^{\int P(t)dt}$.
3. Multiply through: $\mu (t)y^{\prime }+\mu (t)P(t)y=\mu (t)Q(t)$.
4. Recognize the left side as $\frac{d}{dt}[\mu (t)y]$.
5. Integrate: $\mu (t)y=\int \mu (t)Q(t)dt+C$.
6. Solve for $y$: $y=\frac{1}{\mu (t)}\left(\int \mu (t)Q(t)dt+C\right)$.

Example - Mixing Problem: A tank contains 100L of brine with 10kg of dissolved salt. Pure water enters at 5 L/min, and the well-stirred mixture drains at 5 L/min. Find the amount of salt $S(t)$ in the tank. The rate of change of salt is $\frac{dS}{dt}$ = (rate in) - (rate out). Salt enters at 0 kg/min. It drains at a concentration of $\frac{S(t)}{100}$ kg/L times 5 L/min. This gives the linear ODE: $\frac{dS}{dt}=-\frac{5}{100}S=-0.05S$.

1. Standard form: $S^{\prime }+0.05S=0$. (Here, $Q(t)=0$).
2. Integrating factor: $\mu (t)=e^{\int 0.05dt}=e^{0.05t}$.
3. Multiply: $e^{0.05t}{S}^{\prime }+0.05e^{0.05t}S=0$.
4. Left side is $\frac{d}{dt}[e^{0.05t}S]=0$.
5. Integrate: $e^{0.05t}S=C$.
6. Solve: $S(t)=C{e}^{-0.05t}$. Using $S(0)=10$, we find $C=10$, so $S(t)=10e^{-0.05t}$ kg.

Graphical and Numerical Approaches: Slope Fields and Euler's Method

When an ODE is difficult or impossible to solve analytically, we can still understand its solutions graphically or approximate them numerically.

A slope field (or direction field) is a visual tool. For $\frac{dy}{dt}=f(t,y)$, you plot a short line segment with slope $f(t,y)$ at numerous points $(t,y)$ in the plane. The resulting field shows the "flow" of solutions. Any curve that follows these slopes tangentially is a solution curve. This is excellent for visualizing families of solutions and equilibrium solutions.

Euler's Method is a fundamental numerical algorithm for approximating the solution to an initial value problem. It uses the tangent line to extrapolate forward in small steps. Given $\frac{dy}{dt}=f(t,y)$, $y(t_0)=y_0$, and a step size $h$:

1. Set $t_1=t_0+h$, then $y_1=y_0+f(t_0,y_0)\cdot h$.
2. Next, $t_2=t_1+h$, and $y_2=y_1+f(t_1,y_1)\cdot h$.
3. Continue: $y_{n+1}=y_n+f(t_n,y_n)\cdot h$.

While simple, Euler's method introduces error that accumulates with each step. However, it clearly illustrates the concept of building a solution from local rate-of-change information.

Key Modeling Applications

The power of differential equations is revealed in their applications. Newton's Law of Cooling states that the rate of change of an object's temperature is proportional to the difference between its temperature and the ambient temperature. If $T(t)$ is the object's temperature and $T_a$ is ambient temperature, the model is $\frac{dT}{dt}=-k(T-T_a)$, which is both separable and linear, leading to an exponential approach solution: $T(t)=T_a+(T_0-T_a)e^{-kt}$.

Population growth can be modeled in multiple tiers. The simple Malthusian model assumes unlimited resources: $\frac{dP}{dt}=kP$, leading to exponential growth $P(t)=P_0{e}^{kt}$. The more realistic logistic growth model accounts for carrying capacity $M$: $\frac{dP}{dt}=kP(1-\frac{P}{M})$. This separable ODE produces an S-shaped curve where growth slows as the population approaches $M$.

Common Pitfalls

1. Misapplying Separation of Variables: Attempting to separate variables in non-separable equations like $y^{\prime }=t+y$ is a classic error. Always check if you can write the equation as a product of a function of $t$ and a function of $y$. If not, another method, like checking for linearity, is required.
2. Incorrect Integrating Factor: For linear equations $\frac{dy}{dt}+P(t)y=Q(t)$, the integrating factor is $e^{\int P(t)dt}$. A common mistake is to use $P(t)$ itself or to forget to put the equation in standard form first. If the coefficient of $y^{\prime }$ is not 1, you must divide the entire equation by it to identify the correct $P(t)$.
3. Forgotten Constant of Integration: When integrating, the $+C$ is essential. However, a frequent subtle error is adding a constant on both sides during separation of variables. You only need one constant, as combining two constants ($C_2-C_1$) yields a single arbitrary constant $C$.
4. Confusing General and Particular Solutions: The general solution contains an arbitrary constant $C$ and represents a family of curves. A particular solution satisfies both the ODE and a given initial condition, resulting in a specific value for $C$. Failing to use the initial condition to find $C$ when the problem requests a specific solution is a critical oversight.

Summary

• Differential equations model dynamic systems by relating a function to its rate of change. First-order ODEs form the essential foundation for this study.
• Separable equations $\left(\frac{dy}{dt}=g(t)h(y)\right)$ are solved by isolating $dy$ and $dt$ on opposite sides and integrating. This method directly applies to classic models like radioactive decay and unrestricted population growth.
• First-order linear equations $\left(y^{\prime }+P(t)y=Q(t)\right)$ are solved using an integrating factor $\mu (t)=e^{\int P(t)dt}$, which transforms the equation into an easily integrable form. This is key for mixing problems.
• When analytic solutions are challenging, slope fields provide a qualitative graphical picture of solution families, while Euler's Method offers a straightforward numerical approximation.
• These tools unlock powerful modeling applications, including Newton's Law of Cooling (exponential approach) and logistic population growth (S-curve with a carrying capacity), demonstrating how differential equations translate real-world behavior into predictive mathematics.