Rational functions with polynomial denominators are common in control systems, circuit work, and calculus, but rarely convenient to manipulate directly. Partial fraction decomposition breaks these tricky expressions into simpler ones, using numerator coefficients, denominator roots, and factor type—making analysis possible. You’ll see this mainly for inverse Laplace transforms, transfer function breakdowns, and when you need to integrate something that can’t be handled another way. This page lays out the formulas, a worked control example, the reasoning behind each method, and a straight-shooting FAQ.
What is partial fraction decomposition?
Partial fraction decomposition is a way to take a tough rational function—a polynomial divided by another polynomial—and express it as a sum of easier-to-work-with fractions. It’s about breaking a single complicated problem into several manageable ones.
Simple Explanation
If you’ve got a fraction that’s messy to integrate or invert—think of trying to process a whole watermelon at once—partial fractions let you cut it into more workable pieces. Each piece lines up with a denominator factor, and all of them together add back up to the original.
📐 Browse all 1000+ Interactive Calculators
Table of Contents
Visual Diagram
How to Use This Calculator
- Select your decomposition type from the dropdown — choose from distinct linear factors, repeated roots, irreducible quadratic, mixed, improper fraction, or cover-up method.
- Enter the numerator coefficients for your rational expression in the fields that appear.
- Enter the denominator roots (or quadratic coefficients) depending on the mode selected.
- Click Calculate to see your result.
Partial Fractions Interactive Calculator
This calculator is intended for education, concept evaluation, and preliminary design. Results are based on the equations and assumptions described on this page, but cannot account for every real-world load case, tolerance, material property, environmental condition, installation detail, safety factor, code, or regulatory requirement. Verify all inputs, assumptions, units, and results independently before selecting components or using the result in a real application. Safety-critical, structural, medical, lifting, transportation, or regulated applications must be reviewed by a qualified engineer.
Partial Fractions Interactive Visualizer
Decompose complex rational functions into simpler fractions by adjusting numerator coefficients and denominator roots. Watch how different factor types create distinct partial fraction patterns in real-time.
COEFFICIENT A
-12
COEFFICIENT B
17
ORIGINAL DEGREE
1
FIRGELLI Automations — Interactive Engineering Calculators
Mathematical Formulas for Partial Fraction Decomposition
Use the formula below to calculate the partial fraction decomposition for each factor type.
Distinct Linear Factors
P(x) / [(x - r₁)(x - r₂)] = A/(x - r₁) + B/(x - r₂)
Where:
P(x) = numerator polynomial (degree < 2)
r₁, r₂ = distinct real roots of denominator
A, B = constants to be determined
Degree of P(x) must be less than degree of denominator (proper fraction)
Repeated Linear Factor
P(x) / (x - r)ⁿ = A₁/(x - r) + A₂/(x - r)² + ... + Aₙ/(x - r)ⁿ
Where:
P(x) = numerator polynomial (degree < n)
r = repeated root with multiplicity n
A₁, A₂, ..., Aₙ = constants determined by substitution or coefficient comparison
n = multiplicity of the repeated root (n ≥ 2)
Irreducible Quadratic Factor
P(x) / [(x - r)(x² + bx + c)] = A/(x - r) + (Bx + C)/(x² + bx + c)
Where:
x² + bx + c = irreducible quadratic (b² - 4c < 0, no real roots)
r = real root of linear factor
A = constant for linear factor
B, C = constants for quadratic factor (both required)
Discriminant b² - 4c must be negative for true irreducibility
Cover-Up Method (Quick Coefficient Finding)
A = limx→r₁ [(x - r₁) · P(x)/Q(x)]
Where:
Q(x) = denominator polynomial with factor (x - r₁)
The limit is evaluated by substituting x = r₁ after canceling (x - r₁)
Method applies only to distinct linear factors
Provides fastest coefficient computation for simple cases
Improper Fraction Long Division
P(x)/Q(x) = D(x) + R(x)/Q(x)
Where:
deg(P) ≥ deg(Q) = improper fraction condition
D(x) = quotient polynomial from long division
R(x) = remainder polynomial with deg(R) < deg(Q)
R(x)/Q(x) = proper fraction for partial fraction decomposition
Must perform division before applying partial fractions
General Coefficient Determination
P(x) = A·Q₁(x) + B·Q₂(x) + ... (multiply both sides by denominator)
Where:
Q₁(x), Q₂(x) = complementary factors of denominator
Coefficients found by: (1) Strategic substitution of x values, or (2) Equating coefficients of like powers
Strategic x values: choose values that zero out certain terms
Coefficient comparison: match x², x¹, x⁰ terms on both sides
Simple Example
Decompose (5x + 7) / [(x − 1)(x − 2)] using distinct linear factors mode.
- a₁ = 5, a₀ = 7
- r₁ = 1, r₂ = 2
- A = (5·1 + 7) / (1 − 2) = 12 / (−1) = −12
- B = (5·2 + 7) / (2 − 1) = 17 / 1 = 17
- Result: −12/(x − 1) + 17/(x − 2)
Theory & Engineering Applications of Partial Fraction Decomposition
Partial fraction decomposition is a staple in engineering analysis when you’re dealing with rational expressions you can’t work with directly. It’s especially common when solving differential equations, breaking down transfer functions, handling inverse Laplace transforms, or just integrating something awkward. The main point is straightforward: any proper rational function P(x)/Q(x) (with numerator degree less than denominator degree) can be re-written as a sum of terms, where each denominator is a factor (possibly to a power) that you can trace back to the original denominator.
Mathematical Foundation and Uniqueness
Partial fraction decomposition depends on how the denominator factors—linears, repeated roots, or irreducible quadratics (where real factorization isn’t possible). Each case sets the form of the numerators: constants for unique linear terms, multiple terms for repeated roots (each with increasing power in the denominator), and linear numerators for irreducible quadratics.
Here’s a limitation that gets overlooked: you can’t decompose unless you’re able to factor the denominator, and as degree goes up, that’s not a given. For fifth-order and higher polynomials, you may need numerical root-finding, which brings in approximation error—not always ideal when precision matters.
Control Systems and Transfer Function Analysis
If you work with control systems, you’ll end up decomposing transfer functions like H(s) = N(s)/D(s). The poles—roots of D(s)—set the behavior. Partial fractions let you see how each pole affects the time-domain response. For first-order factors, each term is a simple exponential decay. For quadratics, you get oscillations determined by damping ratio and natural frequency.
By doing the decomposition, you can read off resonance, settling times, and more—without needing to grind out the inverse Laplace explicitly. For example, with H(s) = (3s + 7)/(s² + 3s + 2), partial fractions get you 4/(s + 1) - 1/(s + 2). That’s two simple exponentials: one decays with time constant 1 s, the other 0.5 s.
Laplace Transform Inversion in Circuit Analysis
Any circuit engineer using Laplace knows the drill: you map the circuit to Laplace, solve algebraically, then decompose before inversion. If you have an RLC or similar system, breaking the function down gives you components you know how to invert by table. For a transfer function like Vout(s)/Vin(s) = ωₙ²/(s² + 2ζωₙs + ωₙ²), multiply by 1/s for a step input, decompose, and you immediately see both the steady-state and transient portion. The component with 1/s is the final value, the quadratic term gives the oscillations and damping. It’s a routine but essential step—manual or automated—before you can make sense of the time-domain signal.
Integration Techniques and Calculus Applications
For integrals involving rational functions, breaking things into partial fractions is still the most practical approach. Each decomposed term translates to a standard, known integral: logarithmic for simple linears, power law for repeated, and a mix of arctangent/log for quadratics. This also comes up on definite integrals from 0 to infinity—sometimes producing telescoping terms or letting you compute tough results with basic calculus. In some fields (like probability), this is the only practical way to assess characteristic functions and related integrals.
Signal Processing and Filter Design
Partial fractions play a key role in digital filtering—especially when turning a transfer function H(z) into a parallel structure. Each decomposed term can be mapped to a first-order or second-order digital filter. This isn’t just theoretical: using the right decomposition improves both numerical stability and ease of implementation. For example, in IIR filter design, you typically want to split a fourth-order transfer function into two cascaded or parallel second-order sections—easier to tune, and less chance of amplifying numerical errors if your processing uses fixed point.
Advanced Applications: Heat Transfer and Vibration Analysis
Outside of circuits, you see the same decompositions in vibration and heat problems. When you want to analyze or control vibration in a structure (like a satellite panel), you break the frequency response into modal contributions—each mode corresponding to a partial fraction term. In heat transfer, splitting the equation this way gives you access to the thermal time constants and shows you which dynamics dominate the system’s response. You won’t always get closed-form answers for everything, but the structure tells you how the different physical effects add up.
Worked Example: Control System Step Response
Take a transfer function H(s) = 12/(s² + 7s + 12). If the system gets a unit step input, here’s the process:
Step 1: Form the output expression
Y(s) = H(s)·U(s) = 12/[s(s² + 7s + 12)] = 12/[s(s + 3)(s + 4)]
Step 2: Set up partial fraction form
12/[s(s + 3)(s + 4)] = A/s + B/(s + 3) + C/(s + 4)
Step 3: Apply cover-up method for coefficient A
Multiply both sides by s and substitute s = 0:
A = 12/[(0 + 3)(0 + 4)] = 12/12 = 1.0000
Step 4: Apply cover-up method for coefficient B
Multiply both sides by (s + 3) and substitute s = -3:
B = 12/[(-3)(-3 + 4)] = 12/(-3)(1) = -4.0000
Step 5: Apply cover-up method for coefficient C
Multiply both sides by (s + 4) and substitute s = -4:
C = 12/[(-4)(-4 + 3)] = 12/[(-4)(-1)] = 12/4 = 3.0000
Step 6: Write complete decomposition
Y(s) = 1.0000/s - 4.0000/(s + 3) + 3.0000/(s + 4)
Step 7: Compute inverse Laplace transform
y(t) = 1.0000·u(t) - 4.0000·e^(-3t) + 3.0000·e^(-4t)
Step 8: Verify steady-state value
As t → ∞: y(∞) = 1.0000 - 0 + 0 = 1.0000 ✓ (matches final value theorem: lim[s→0] sY(s) = 1)
Step 9: Determine settling time (2% criterion)
Dominant pole at s = -3 gives τ = 1/3 = 0.3333 seconds
Settling time ts = 4τ = 4(0.3333) = 1.3333 seconds
Physical Interpretation: The response starts at y(0) = 1 - 4 + 3 = 0, creeping up to 1.0, with the e^(-4t) part dying out fastest while e^(-3t) controls the longer settling. If B is negative (like -4), you see an undershoot right after t=0, typical for overdamped systems. You can estimate what’s important (like 75% of transient energy in the s = -3 term here) just by looking at the coefficients.
Partial fraction decomposition isn’t only a calculation trick—it lets you pick apart what each part of your system does, giving quick insight into which dynamics matter. Mastering it gets you from theory to practical understanding fast. If you’re working repeatedly with transfer functions or tough integrals, this is a tool to keep sharp. More tools are available in the engineering calculators hub.
Practical Applications
Scenario: Electronics Engineer Designing Active Filter
Marcus, a senior electronics engineer at an audio equipment company, is designing a fourth-order Butterworth lowpass filter for a high-fidelity preamplifier. The filter's transfer function in the s-domain is H(s) = ω₀⁴/(s⁴ + 2.613ω₀s³ + 3.414ω₀²s² + 2.613ω₀³s + ω₀⁴) with cutoff frequency ω₀ = 2π(20000) rad/s. To implement this as cascaded second-order Sallen-Key stages (which are more stable and easier to tune than a single fourth-order stage), Marcus uses the partial fractions calculator to decompose the transfer function into two second-order sections. He enters the numerator and denominator coefficients, and the calculator reveals the optimal pole pairing: two sections with Q-factors of 0.541 and 1.307. This decomposition tells Marcus exactly which resistor and capacitor values to use in each stage, and importantly, which stage should come first (lower Q) to minimize signal clipping in the intermediate stages. The partial fraction coefficients directly translate to component values that achieve the desired 80 dB/decade rolloff beyond 20 kHz while maintaining flat response in the passband.
Scenario: Mechanical Engineer Analyzing Structural Vibration
Dr. Yuki Tanaka, a vibration specialist at an aerospace manufacturer, is investigating resonance issues in a satellite solar panel deployment mechanism. Accelerometer data shows a complex frequency response function with multiple peaks, represented as a rational function H(ω) = (2.3ω² + 15.7)/(ω⁴ + 18.6ω³ + 892ω² + 7440ω + 14400). Traditional frequency-domain analysis shows peaks, but Yuki needs to identify the individual modal contributions—each representing a physical vibration mode with its own frequency and damping. She uses the partial fractions calculator to decompose this into four first-order complex conjugate pairs. The decomposition reveals that the system has natural frequencies at 6.8 Hz and 12.3 Hz, with damping ratios of 0.12 and 0.08 respectively. The residue magnitudes (partial fraction coefficients) show that the 12.3 Hz mode carries 73% of the vibrational energy, identifying it as the critical mode to address. Based on these coefficients, Yuki recommends adding a 47-gram tuned mass damper at the panel's quarter-span point—a targeted solution that emerged directly from the partial fraction decomposition rather than requiring expensive iterative prototyping.
Scenario: Graduate Student Solving Laplace Transform Problem
Emma, a chemical engineering graduate student, is modeling the transient response of a continuous stirred-tank reactor (CSTR) with recycle stream. Her differential equation model, after Laplace transformation, yields the concentration profile C(s) = (5s + 17)/[s(s² + 6s + 8)]. Emma needs the time-domain solution c(t) to validate her model against experimental data collected every 30 seconds for 5 minutes. She enters her expression into the partial fractions calculator, which decomposes it into C(s) = 2.125/s - 0.625/(s + 2) - 1.500/(s + 4). Each term has a known inverse Laplace transform: the first term (2.125/s) represents the steady-state concentration of 2.125 mol/L, the second term gives a transient component -0.625e^(-2t) that decays with a 0.5-minute time constant, and the third term contributes -1.500e^(-4t) that vanishes within 0.25 minutes. Emma plots c(t) = 2.125 - 0.625e^(-2t) - 1.500e^(-4t) against her experimental data and achieves an R² correlation of 0.97, confirming her reactor model captures the essential dynamics. The partial fraction decomposition not only solved her mathematical problem but also revealed that the fast mode (4-minute⁻¹ decay rate) dominates early response while the slow mode (2-minute⁻¹) controls the approach to equilibrium—insights that inform her recommendations for optimal reactor startup procedures.
Frequently Asked Questions
What should I do if my rational function is improper (numerator degree ≥ denominator degree)? +
How do I handle complex conjugate roots in partial fraction decomposition? +
What is the cover-up method and when should I use it versus coefficient comparison? +
How many terms appear in the partial fraction decomposition for a repeated root of multiplicity n? +
Can partial fractions be used with functions other than polynomials in the denominator? +
What are the most common errors when performing partial fraction decomposition by hand? +
Free Engineering Calculators
Explore our complete library of free engineering and physics calculators.
Browse All Calculators →🔗 Explore More Free Engineering Calculators
About the Author
Robbie Dickson — Chief Engineer & Founder, FIRGELLI Automations
Robbie Dickson brings over two decades of engineering expertise to FIRGELLI Automations. With a distinguished career at Rolls-Royce, BMW, and Ford, he has deep expertise in mechanical systems, actuator technology, and precision engineering.
📹 Video Walkthrough — Partial Fractions Interactive Calculator
📹 Video Walkthrough — Partial Fractions Interactive Calculator
Need to implement these calculations?
Explore the precision-engineered motion control solutions used by top engineers.
