ODE: Linear First-Order Equations

Linear first-order ordinary differential equations are a cornerstone of mathematical modeling in engineering. They describe a direct relationship between a rate of change and the current state of a system, making them indispensable for analyzing dynamic processes in circuits, fluid flow, chemical reactions, and thermal systems. Mastering the integrating factor method to solve these equations is a fundamental skill that bridges pure mathematics and practical engineering design.

The Standard Form and Its Logic

A first-order linear ordinary differential equation is defined by its linearity in the unknown function $y$ and its derivative $dy/dx$. It must be writable in the standard form:

$$\frac{dy}{dx}+P(x)y=Q(x)$$

Here, $P(x)$ and $Q(x)$ are known functions of the independent variable $x$. The term "linear" refers to the fact that $y$ and $dy/dx$ appear to the first power and are not multiplied together. The function $Q(x)$ is called the forcing function or input. If $Q(x)=0$, the equation is homogeneous; otherwise, it is nonhomogeneous.

The challenge is that we cannot simply integrate both sides with respect to $x$ because of the presence of the term $P(x)y$. The left side is not the derivative of a simple product. The core insight is to multiply the entire equation by a specially chosen function, called an integrating factor, that transforms the left-hand side into the exact derivative of a product. This is the genius of the method: it ingeniously reverses the product rule of differentiation.

Deriving and Using the Integrating Factor

Our goal is to find a function $\mu (x)$ such that multiplying the standard form by it yields:

$$\mu (x)\frac{dy}{dx}+\mu (x)P(x)y=\frac{d}{dx}\left[\mu (x)y\right]$$

By the product rule, the right side expands to $\mu (x)\frac{dy}{dx}+\frac{d\mu }{dx}y$. Comparing terms, we see that for the equality to hold, we must have:

$$\frac{d\mu }{dx}=\mu (x)P(x)$$

This is a separable equation for $\mu (x)$. Solving it gives:

$$\int \frac{d\mu }{\mu }=\int P(x)dx\Rightarrow \ln |\mu |=\int P(x)dx$$

Therefore, the integrating factor $\mu (x)$ is:

$$\mu (x)=e^{\int P(x)dx}$$

When we multiply our original equation in standard form by this $\mu (x)$, the left side condenses into the derivative of the product $\mu (x)y$:

$$\frac{d}{dx}\left[\mu (x)y\right]=\mu (x)Q(x)$$

We can now simply integrate both sides with respect to $x$:

$$\mu (x)y=\int \mu (x)Q(x)dx+C$$

Finally, solving for $y$ gives the general solution formula:

$$y(x)=\frac{1}{\mu (x)}\left(\int \mu (x)Q(x)dx+C\right)$$

Critical Step: When computing the integrating factor, the constant of integration is always set to zero. Using any non-zero constant would simply multiply the entire equation by a scalar, canceling itself out in the final solution. The essential antiderivative is all that's required.

Applying Initial Conditions and Interpreting Solutions

The general solution contains an arbitrary constant $C$. An initial condition, typically given as $y(x_0)=y_0$, allows us to find a unique, particular solution. You substitute the initial condition into the general solution after the integration step to solve for $C$.

For nonhomogeneous equations ($Q(x)\ne 0$), the solution can often be decomposed into two conceptual parts: transient and steady-state components. The transient component is the part of the solution that decays to zero as $x$ (or more commonly time, $t$) increases. It is associated with the complementary solution to the homogeneous equation and contains the constant $C$. The steady-state component is the part that remains after the transient dies out; it is a particular solution that persists, often mirroring the form of the forcing function $Q(x)$. This decomposition is crucial in engineering for understanding the short-term versus long-term behavior of a system.

Application 1: Mixing (Concentration) Problems

A classic application involves a tank containing a solution (e.g., salt water). A brine of a certain concentration enters the tank at a fixed rate, the mixture is kept uniform by stirring, and it exits the tank at another rate. The goal is to find the amount of solute (e.g., salt) $A(t)$ in the tank at any time $t$.

The governing principle is the balance equation: the rate of change of the amount of solute equals the rate in minus the rate out.

$$\frac{dA}{dt}=(\text{concentration in})\times (\text{flow rate in})-\left(\frac{A(t)}{\text{volume in tank at time }t}\right)\times (\text{flow rate out})$$

The term $A(t)/\text{volume}$ is the time-varying concentration of the exiting mixture. The volume in the tank may be constant or changing depending on the input and output flow rates. This equation almost always fits the linear first-order standard form $\frac{dA}{dt}+P(t)A=Q(t)$, where $P(t)$ and $Q(t)$ are constructed from the flow rates and concentrations. Solving it with an integrating factor predicts the concentration history in the tank.

Application 2: Analysis of an RC Circuit

Consider a simple RC circuit with a resistor $R$, capacitor $C$, and a voltage source $V_s(t)$. Using Kirchhoff's voltage law, the sum of the voltages around the loop is zero: $V_R+V_C-V_s=0$.

The voltage across the resistor is $V_R=iR$. The voltage across the capacitor is $V_C=q/C$, where $q$ is the charge on the capacitor. Since current $i$ is the rate of flow of charge, $i=dq/dt$. Substituting these into the voltage law gives:

$$R\frac{dq}{dt}+\frac{1}{C}q=V_s(t)$$

This is a linear first-order ODE for the charge $q(t)$, with $P(t)=1/(RC)$ and $Q(t)=V_s(t)/R$. The integrating factor is $\mu (t)=e^{\int (1/RC)dt}=e^{t/(RC)}$. The constant $RC$ has units of time and is famously known as the time constant $\tau$ of the circuit. It dictates the speed of the circuit's response. The solution $q(t)$ will show a transient exponential approach from the initial charge to a steady-state charge determined by $V_s(t)$. The voltage across the capacitor, $V_C(t)=q(t)/C$, follows the same dynamic profile.

Common Pitfalls

1. Incorrect Standard Form: The most common error is misidentifying $P(x)$. The equation must be in the form $dy/dx+P(x)y=Q(x)$. If the coefficient of $dy/dx$ is not 1, you must divide the entire equation by that coefficient to put it in standard form before identifying $P(x)$. For example, in $2y^{\prime }+3xy=x^2$, first divide by 2 to get $y^{\prime }+(3x/2)y=x/2$, so $P(x)=3x/2$, not 3x.

1. Integrating Factor Constant: Adding a constant of integration when computing $\int P(x)dx$ for the integrating factor. Remember, $\mu (x)=e^{\int P(x)dx}$, not $e^{\int P(x)dx+C}$. Including $C$ creates $\mu (x)=e^C{e}^{\int P(x)dx}$, an extra multiplicative constant that complicates the algebra but cancels out in the end. It's an unnecessary step.

1. Forgetting the Constant of Integration $C$: After integrating $\frac{d}{dx}[\mu y]=\mu Q$, you must add the constant $C$ on the right side: $\mu y=\int \mu Qdx+C$. Forgetting $C$ at this stage will make applying an initial condition impossible and yields an incomplete general solution.

1. Misapplying the Initial Condition: The initial condition $y(x_0)=y_0$ must be applied to the general solution to solve for $C$. A frequent mistake is substituting it into an intermediate step (like before integrating $\mu Q$) or after simplifying an expression where $C$ has been conflated with other constants.

Summary

• A first-order linear ODE must be manipulated into the standard form $dy/dx+P(x)y=Q(x)$ before applying the integrating factor method.
• The integrating factor $\mu (x)=e^{\int P(x)dx}$ is designed to make the left side of the ODE a perfect derivative, leading to the general solution: $y(x)=\frac{1}{\mu (x)}\left(\int \mu (x)Q(x)dx+C\right)$.
• Initial conditions are used to find the constant $C$ and obtain a unique particular solution, which can often be interpreted in terms of decaying transient and persistent steady-state components.
• Engineering applications like mixing problems and RC circuit analysis are directly modeled by these equations, where the parameters (like flow rates or the time constant $RC$) define the system's dynamic behavior.
• Avoiding common algebraic mistakes, especially in putting the equation in standard form and handling constants of integration, is critical for reliable problem-solving.