Study guide

Differential Equations

FE Fundamentals of Engineering Exam

Track your progress as you read

Finished sections get a check mark.

1 Introduction to the FE Exam

About the FE Exam FE Exam Study Strategies Navigating the FE Reference Handbook

2 Mathematics (8–12 Questions)

Analytic Geometry and Trigonometry Single-Variable Calculus Differential Equations Linear Algebra and Matrix Operations Numerical Methods

3 Probability and Statistics (6–9 Questions)

Probability Fundamentals Probability Distributions Statistical Estimation and Hypothesis Testing

4 Chemistry (5–8 Questions)

Chemical Reactions and Stoichiometry Acids, Bases, and pH Oxidation-Reduction and Corrosion

5 Ethics and Professional Practice (5–8 Questions)

Engineering Ethics Fundamentals

6 Engineering Economics (6–9 Questions)

Time Value of Money Economic Analysis Methods

7 Statics (9–14 Questions)

Force Systems and Equilibrium Trusses, Frames, and Machines Centroids and Moments of Inertia Static Friction

8 Dynamics (9–14 Questions)

Particle Kinematics Kinetics: Force, Mass, and Acceleration Work, Energy, Impulse, and Momentum Rigid-Body Rotation and Mechanical Vibrations

9 Strength of Materials (9–14 Questions)

Stress and Strain Fundamentals Beam Analysis: Shear and Moment Diagrams Bending Stress and Torsion Stress Transformation and Material Failure

10 Materials Science (6–9 Questions)

Material Properties and Selection Phase Diagrams and Alloy Systems

11 Fluid Mechanics (12–18 Questions)

Fluid Properties and Statics Fluid Dynamics and Bernoulli's Equation Pumps, Turbines, and Flow Measurement

12 Thermodynamics and Heat Transfer (9–14 Questions)

Laws of Thermodynamics and Properties Thermodynamic Cycles and Applications Heat Transfer: Conduction, Convection, and Radiation

13 Basic Electrical Engineering (6–9 Questions)

Electrical Fundamentals and DC Circuits AC Circuits and Three-Phase Power

14 Safety, Health, and Environment (6–9 Questions)

Industrial Safety and Hazard Management

15 Instrumentation and Controls (4–6 Questions)

Sensors, Data Acquisition, and Control Systems

16 Computational Tools and Computer Science

Spreadsheets, Structured Programming, and Flowcharts

Differential Equations

Key Takeaways Explain This Topic AI

First-order ODEs split into separable, linear (integrating factor μ = e^∫P dx), and exact (∂M/∂y = ∂N/∂x) types.

Second-order linear ODEs with constant coefficients are solved from the characteristic equation ar² + br + c = 0.

The discriminant b² − 4ac selects the solution form: distinct real roots, a repeated root, or complex conjugates.

Complex roots α ± βi give an oscillatory solution y = e^(αx)(C₁ cos βx + C₂ sin βx).

Laplace transforms turn an ODE into algebra: ℒ{f′} = sF(s) − f(0); look up the transform pairs in the Handbook.

An nth-order ODE needs n initial/boundary conditions to pin down all the arbitrary constants.

Last updated: June 2026

Differential equations model how engineering quantities change — RC/RL circuit currents, spring-mass vibrations, reactor concentrations, transient heat conduction. The FE expects you to identify the ODE type, pick the matching solution method, and apply initial conditions , leaning on the Laplace-transform and solution-form tables in the Handbook.

  • First-Order ODEs
  • Separable

Writable as dy/dx = f(x)·g(y); separate and integrate: ∫dy/g(y) = ∫f(x) dx.

Worked example — separable IVP. Solve dy/dx = 3x²y with y(0) = 2. Separate: dy/y = 3x² dx. Integrate both sides: ln|y| = x³ + C, so y = A·e^(x³). Apply y(0) = 2: A·e⁰ = A = 2. Thus y = 2e^(x³) . Note the exponent is x³ (the antiderivative of 3x²), not 3x — confusing these is the usual mistake.

First-Order Linear

Standard form dy/dx + P(x)y = Q(x). Multiply by the integrating factor μ(x) = e^∫P(x) dx, then y = (1/μ)∫μ Q dx.

Worked example — integrating factor. Solve dy/dx + 2y = 6. Here P = 2 so μ = e^(2x). Then y = e^(−2x)∫6e^(2x) dx = e^(−2x)(3e^(2x) + C) = 3 + Ce^(−2x) . As x → ∞ the transient Ce^(−2x) decays and y approaches the steady-state value 3.

Exact

M(x,y)dx + N(x,y)dy = 0 is exact when ∂M/∂y = ∂N/∂x; then a potential F(x,y) exists with ∂F/∂x = M and ∂F/∂y = N, and the solution is F(x,y) = C. To build F, integrate M with respect to x, then differentiate the result with respect to y and match it to N to recover the y-only term.

Recognizing the Type Quickly

On the FE the hardest part is classification , so scan in this order: (1) Can you split into f(x)·g(y)? → separable. (2) Is it linear in y of the form y′ + P(x)y = Q(x)? → integrating factor. (3) Otherwise test ∂M/∂y = ∂N/∂x for exactness. A first-order linear ODE with constant coefficient and constant forcing (like the integrating-factor example above) always settles to a steady state plus a decaying transient — a pattern that mirrors RC-circuit charging and first-order sensor response, both of which the FE poses as word problems.

Second-Order Linear ODEs (Constant Coefficients)

General form ay″ + by′ + cy = f(x). For the homogeneous case (f = 0), solve the characteristic equation ar² + br + c = 0. The discriminant b² − 4ac decides the solution shape:

  • Discriminant Roots General solution
  • b² − 4ac > 0 real, distinct r₁ ≠ r₂ y = C₁e^(r₁x) + C₂e^(r₂x)
  • b² − 4ac = 0 repeated r y = (C₁ + C₂x)e^(rx)
  • b² − 4ac < 0 complex α ± βi y = e^(αx)(C₁ cos βx + C₂ sin βx)

Worked example — distinct real roots. Solve y″ − 5y′ + 6y = 0. Characteristic equation r² − 5r + 6 = (r−2)(r−3) = 0 gives r = 2, 3. Solution: y = C₁e^(2x) + C₂e^(3x) .

For the non-homogeneous case, the full solution is y = yₕ + yₚ. Guess yₚ by undetermined coefficients :

  • f(x) Guess for yₚ
  • constant k A
  • polynomial degree n Aₙxⁿ + … + A₀
  • e^(αx) Ae^(αx)
  • sin βx or cos βx A cos βx + B sin βx

If a guessed term already appears in yₕ, multiply the guess by x (or x²) to keep it independent.

For a mechanical analogy, b² − 4ac < 0 corresponds to an underdamped vibration (oscillating decay), = 0 is critically damped , and > 0 is overdamped — useful for recognizing which root case a physical problem implies.

Laplace Transforms

The Laplace transform converts a time function f(t) into F(s) = ∫₀^∞ f(t)e^(−st) dt, turning calculus into algebra. The Handbook prints the transform pairs and properties; your job is to apply them.

  • f(t) F(s) f(t) F(s)
  • 1 1/s e^(at) 1/(s−a)
  • t 1/s² sin ωt ω/(s²+ω²)
  • tⁿ n!/s^(n+1) cos ωt s/(s²+ω²)

Derivative properties: ℒ{f′} = sF(s) − f(0); ℒ{f″} = s²F(s) − s f(0) − f′(0). s-shift: ℒ{e^(at)f(t)} = F(s−a).

Procedure

Laplace-transform both sides of the ODE.

Substitute the initial conditions.

Solve algebraically for Y(s).

Inverse-transform (match a table pair) to get y(t).

Worked example — Laplace solution. Solve y″ + 4y = 0 with y(0) = 1, y′(0) = 0.

Transform: s²Y − s(1) − 0 + 4Y = 0.

Collect: Y(s²+4) = s, so Y(s) = s/(s²+4).

Match the pair s/(s²+ω²) with ω² = 4 → ω = 2.

Inverse: y(t) = cos 2t .

The transform method shines for problems with discontinuous forcing (step or impulse inputs) where undetermined coefficients become awkward.

Quick transform recall. For f(t) = e^(3t), match e^(at) → 1/(s−a) with a = 3, giving F(s) = 1/(s−3), valid for s > 3. The sign trap: it is s minus a, so e^(+3t) → 1/(s−3), not 1/(s+3).

Why Engineers Use Laplace

Beyond the algebra, the Laplace domain exposes a system's poles — the values of s that make the denominator of Y(s) zero. Poles in the left half-plane (negative real part) mean the time response decays and the system is stable ; poles on the imaginary axis give sustained oscillation; right-half-plane poles mean the response grows without bound. This connects directly to the second-order root cases above: complex characteristic roots α ± βi become poles at exactly s = α ± βi, so the discriminant test and the pole location carry the same information. The FE may ask you to read stability straight off the denominator of a transfer function rather than to fully invert the transform.