Linear Equations

Definition and Standard Form

First Order Linear ODE

Form: dy/dx + P(x)y = Q(x). Variables: x (independent), y (dependent). Coefficients: P and Q continuous on interval I.

Linearity Explanation

Equation linear in y and dy/dx. No powers or products of y and its derivatives. Superposition principle applies to homogeneous case.

Examples

dy/dx + 3y = 6x, dy/dt - (2/t)y = t², dy/dx + sin(x)y = cos(x).

Integrating Factor Method

Purpose

Transforms non-exact linear ODE into exact differential. Enables direct integration.

Integrating Factor Formula

μ(x) = exp(∫P(x) dx). Always positive. Multiplies entire equation.

Procedure

Multiply ODE by μ(x). Left side becomes d/dx[μ(x)y]. Integrate both sides: ∫d/dx[μ(x)y] dx = ∫μ(x)Q(x) dx.

μ(x) = e^{∫P(x) dx}d/dx[μ(x)y] = μ(x)Q(x)Solution: y = (1/μ(x)) ∫ μ(x) Q(x) dx + C/μ(x) 

Example

dy/dx + (2/x)y = sin(x)/x, μ(x) = e^{∫2/x dx} = x².

Homogeneous Linear Equations

Definition

Q(x) = 0 in dy/dx + P(x)y = Q(x). Equation reduces to dy/dx + P(x)y = 0.

Solution

Separable: dy/dx = -P(x)y. Integrate: ln|y| = -∫P(x) dx + C.

Result

y = Ce^{-∫P(x) dx} 

Properties

Solution space: one-dimensional vector space. Superposition applies. Stability depends on sign of P(x).

Nonhomogeneous Linear Equations

General Form

dy/dx + P(x)y = Q(x), Q(x) ≠ 0. Solution: sum of complementary solution and particular solution.

Complementary Solution

Solution to associated homogeneous equation dy/dx + P(x)y = 0.

Particular Solution Methods

Integrating factor, variation of parameters, undetermined coefficients (for constant coefficients).

Example

dy/dx + y = e^x, complementary y_c = Ce^{-x}, particular y_p = x e^x, general y = y_c + y_p.

Initial Value Problems (IVP)

Definition

Find solution y(x) satisfying ODE and initial condition y(x_0) = y_0.

Existence and Uniqueness

Conditions: P(x), Q(x) continuous near x_0. Picard–Lindelöf theorem guarantees unique solution.

Procedure

Solve general solution, then impose y(x_0) = y_0 to find constant C.

Example

dy/dx + 2y = 4, y(0) = 1; solution y = Ce^{-2x} + 2; impose y(0)=1 → C = -1; final y = 2 - e^{-2x}.

General Solution Structure

Linear Superposition

General solution = complementary solution + particular solution.

Complementary Solution

General form: y_c = C e^{-∫P(x) dx}, arbitrary constant C.

Particular Solution

Any specific solution y_p satisfying full ODE without arbitrary constants.

Formula Summary

y = y_c + y_py_c = Ce^{-∫P(x) dx}y_p = (1/μ(x)) ∫ μ(x) Q(x) dx 

Variation of Parameters Method

Concept

Find particular solution by replacing constant C with function u(x). Suitable for variable coefficients.

Derivation

Assume y_p = u(x) y_c. Differentiate and substitute into ODE to solve for u'(x).

Formula

y_c = e^{-∫P(x) dx}u'(x) = Q(x) / y_cy_p = y_c ∫ [Q(x) / y_c] dx 

Advantages

Applicable when undetermined coefficients fail. Works with non-constant coefficients.

Applications in Science and Engineering

Electrical Circuits

RC circuits: voltage/current modeled by first order linear ODEs.

Population Dynamics

Simple growth/decay models with linear terms in population.

Chemical Kinetics

First order reaction rates represented as linear ODEs.

Mechanical Systems

Damped harmonic oscillator reduction under certain assumptions.

Heat Transfer

Newton’s law of cooling modeled by linear ODEs.

Relation to Exact Equations

Exact Equation Definition

Form M(x,y) + N(x,y) dy/dx = 0 with ∂M/∂y = ∂N/∂x.

Linear ODE as Exact

Multiplying linear ODE by integrating factor yields exact differential equation.

Verification

After multiplying by μ(x), left side d/dx[μ(x)y]. Equates to exact derivative.

Summary Table

Systems of Linear First Order ODEs

Definition

Vector form: dY/dx = A(x)Y + B(x), Y vector of functions, A matrix of coefficients.

Diagonal and Triangular Systems

Reducible to single ODEs or solved sequentially when A is triangular or diagonal.

Matrix Exponential Solution

Fundamental matrix Φ(x) satisfies dΦ/dx = A(x)Φ(x). General solution Y = Φ(x) C + Φ(x) ∫ Φ^{-1}(x) B(x) dx.

Applications

Coupled mechanical/electrical systems, chemical reaction networks, population models with multiple species.

Numerical Methods for Linear Equations

Euler’s Method

Explicit: y_{n+1} = y_n + h f(x_n, y_n). Simple but low accuracy.

Runge-Kutta Methods

Higher order accuracy: RK4 common. Balances accuracy and complexity.

Stiff Equations

Implicit methods (Backward Euler, trapezoidal) preferred for stability.

Software Tools

MATLAB ode45, Python SciPy integrate.solve_ivp, Mathematica NDSolve.

Common Mistakes and Misconceptions

Ignoring Domain of P and Q

P(x), Q(x) must be continuous on interval; otherwise solution may not exist or be unique.

Incorrect Integrating Factor

Sign errors in ∫P(x) dx lead to wrong μ(x) and incorrect solutions.

Forgetting Constant of Integration

Particularly in indefinite integrals; constant critical for general solution.

Confusing Particular and Complementary Solutions

Particular solution has no arbitrary constants; complementary solution is homogeneous solution with constants.

Misapplication of Variation of Parameters

Must use correct complementary solution and integrate carefully.

References

• Boyce, W.E. & DiPrima, R.C., Elementary Differential Equations and Boundary Value Problems, Wiley, 10th ed., 2012, pp. 45-78.
• Zill, D.G., A First Course in Differential Equations with Modeling Applications, Brooks Cole, 11th ed., 2017, pp. 112-150.
• Kreyszig, E., Advanced Engineering Mathematics, Wiley, 10th ed., 2011, pp. 360-400.
• Tenenbaum, M. & Pollard, H., Ordinary Differential Equations, Dover Publications, 1985, pp. 50-90.
• Butcher, J.C., Numerical Methods for Ordinary Differential Equations, Wiley, 2nd ed., 2008, pp. 25-60.