Laplace Transforms

Study Sheet

Laplace Transforms

Turning calculus into algebra

The Transform Method

Tip
The transform and the short table
ty1234560.5f = t
L{t}\mathcal L\{t\} at s=1s = 1: the falling weight e−te^{-t} (dashed) times f=tf = t gives te−tte^{-t}, whose shaded area is exactly 11.

L{f}(s)=∫0∞f(t)e−stdt\mathcal L\{f\}(s) = \int_0^\infty f(t)e^{-st}dt. The working table: L{eat}=1s−a\mathcal L\{e^{at}\} = \tfrac{1}{s-a}, L{tn}=n!sn+1\mathcal L\{t^n\} = \tfrac{n!}{s^{n+1}}, L{sin⁡bt}=bs2+b2\mathcal L\{\sin bt\} = \tfrac{b}{s^2+b^2}, L{cos⁡bt}=ss2+b2\mathcal L\{\cos bt\} = \tfrac{s}{s^2+b^2}.

Each entry is one integral. For eate^{at}, ∫0∞e−(s−a)t dt=1s−a\int_0^\infty e^{-(s-a)t}\,dt = \tfrac{1}{s - a} when s>as > a (proved below). The trigonometric entries follow from the same integral through Euler's formula: cos⁡bt\cos bt and sin⁡bt\sin bt are the real and imaginary parts of eibte^{ibt}, whose transform is 1s−ib=s+ibs2+b2\tfrac{1}{s - ib} = \tfrac{s + ib}{s^2 + b^2}. The transform weighs f(t)f(t) by the decaying factor e−ste^{-st} and adds up the result; the figure shows this for f(t)=tf(t) = t at s=1s = 1, where the shaded area under te−tte^{-t} is L{t}(1)=112=1\mathcal L\{t\}(1) = \tfrac{1}{1^2} = 1.

Tip
Derivatives become multiplication

L{y′}=s L{y}−y(0)\mathcal L\{y'\} = s\,\mathcal L\{y\} - y(0): differentiation turns into multiplying by ss, with the initial condition built in automatically. An initial-value problem becomes algebra in ss — solve for L{y}\mathcal L\{y\}, split by partial fractions, read the table backward.

Why multiplication by ss: integrating by parts moves the derivative off yy and onto e−ste^{-st}, and differentiating e−ste^{-st} brings down the factor −s-s; the boundary term at t=0t = 0 leaves −y(0)-y(0) behind (see the proof below). So every derivative in the equation becomes a power of ss times YY plus known initial values.

Example
Reading the table backward

F(s)=6s4F(s) = \tfrac{6}{s^4}: since L{t3}=3!s4=6s4\mathcal L\{t^3\} = \tfrac{3!}{s^4} = \tfrac{6}{s^4}, the inverse transform is f(t)=t3f(t) = t^3. Recognizing table entries inside partial-fraction pieces is the entire skill.

Reading backward is legitimate because two continuous functions with the same transform are equal (Lerch's theorem), so each transform has only one continuous inverse. Constants must match exactly: 6s4\tfrac{6}{s^4} is L{t3}\mathcal L\{t^3\} because 3!=63! = 6, while 1s4\tfrac{1}{s^4} would be L{t36}\mathcal L\left\{\tfrac{t^3}{6}\right\}.

Side note
Significance: the engineer’s transform

Control theory runs on Laplace: a system's transfer function H(s)H(s) is the transform of its impulse response, and stability is "poles in the left half-plane" — statements about ss, not tt. Circuit analysis, signal processing, and PID controller design all speak this language natively.

Try it
Try it: one transform, no tables open

(a) Compute L{5e−2t+3t}\mathcal L\{5e^{-2t} + 3t\}. (b) Which function has transform 4s2+4\dfrac{4}{s^2 + 4}?

Show the answer

(a) The transform is linear, so transform term by term: L{e−2t}=1s+2\mathcal L\{e^{-2t}\} = \tfrac{1}{s + 2} and L{t}=1s2\mathcal L\{t\} = \tfrac{1}{s^2}, giving 5s+2+3s2\dfrac{5}{s + 2} + \dfrac{3}{s^2} (valid for s>0s > 0). (b) The sine entry bs2+b2\tfrac{b}{s^2 + b^2} with b=2b = 2 is 2s2+4\tfrac{2}{s^2 + 4}, so 4s2+4=2⋅2s2+4\tfrac{4}{s^2 + 4} = 2\cdot\tfrac{2}{s^2 + 4} is the transform of 2sin⁡2t2\sin 2t. The cosine entry ss2+4\tfrac{s}{s^2 + 4} has ss in the numerator, so it does not match.

Partial Fractions and the Shift Rule

Tip
The s-shift, and the three kinds of pieces

L{eatf(t)}=F(s−a)\mathcal L\{e^{at}f(t)\} = F(s - a): a transform in which ss appears only as s−as - a inverts to eate^{at} times the inverse of the unshifted function. With partial fractions, the pieces are

  • distinct linear factors: As−a+Bs−b\dfrac{A}{s - a} + \dfrac{B}{s - b}, which invert to Aeat+BebtAe^{at} + Be^{bt};
  • a repeated factor: 1(s−a)2\dfrac{1}{(s - a)^2}, which inverts to teatte^{at};
  • an irreducible quadratic: complete the square, e.g. s2+2s+5=(s+1)2+4s^2 + 2s + 5 = (s + 1)^2 + 4, and aim at s+1(s+1)2+4→e−tcos⁡2t\dfrac{s + 1}{(s + 1)^2 + 4} \to e^{-t}\cos 2t and 2(s+1)2+4→e−tsin⁡2t\dfrac{2}{(s + 1)^2 + 4} \to e^{-t}\sin 2t.

The shift rule is one line: ∫0∞e−steatf(t) dt=∫0∞e−(s−a)tf(t) dt=F(s−a)\int_0^\infty e^{-st}e^{at}f(t)\,dt = \int_0^\infty e^{-(s - a)t}f(t)\,dt = F(s - a). Multiplying by an exponential in tt slides the transform in ss. That is why completed squares appear: the transform of the damped oscillation e−tcos⁡2te^{-t}\cos 2t is the transform of cos⁡2t\cos 2t with ss replaced by s+1s + 1. The repeated-factor entry is the shift rule applied to L{t}=1s2\mathcal L\{t\} = \tfrac1{s^2}.

Example
Worked: distinct linear factors

Invert Y(s)=3s+1(s−1)(s+2)Y(s) = \dfrac{3s + 1}{(s - 1)(s + 2)}.

Step 1 — write 3s+1(s−1)(s+2)=As−1+Bs+2\dfrac{3s + 1}{(s - 1)(s + 2)} = \dfrac{A}{s - 1} + \dfrac{B}{s + 2} and clear denominators: 3s+1=A(s+2)+B(s−1)3s + 1 = A(s + 2) + B(s - 1).

Step 2 — choose ss to remove one unknown at a time: s=1s = 1 gives 4=3A4 = 3A, and s=−2s = -2 gives −5=−3B-5 = -3B.

Step 3 — A=43A = \tfrac43 and B=53B = \tfrac53.

Step 4 — invert with the eate^{at} entry: y(t)=43et+53e−2ty(t) = \tfrac43e^{t} + \tfrac53e^{-2t}.

Step 5 — check: y(0)=3y(0) = 3, which matches lim⁡s→∞sY(s)=3\lim_{s\to\infty} sY(s) = 3 (the initial-value theorem).

Example
Worked: completing the square
ty123456−11
The inverse e−t(cos⁡2t+sin⁡2t)e^{-t}(\cos 2t + \sin 2t) (solid) is the unshifted inverse cos⁡2t+sin⁡2t\cos 2t + \sin 2t (dashed) damped by e−te^{-t}.

Invert Y(s)=s+3s2+2s+5Y(s) = \dfrac{s + 3}{s^2 + 2s + 5}.

Step 1 — the denominator has no real roots (4−20<04 - 20 < 0), so complete the square: s2+2s+5=(s+1)2+4s^2 + 2s + 5 = (s + 1)^2 + 4.

Step 2 — rewrite the numerator in terms of s+1s + 1: s+3=(s+1)+2s + 3 = (s + 1) + 2.

Step 3 — split: Y=s+1(s+1)2+4+2(s+1)2+4Y = \dfrac{s + 1}{(s + 1)^2 + 4} + \dfrac{2}{(s + 1)^2 + 4}.

Step 4 — each piece is a table entry with ss replaced by s+1s + 1, i.e. the shift a=−1a = -1: y(t)=e−tcos⁡2t+e−tsin⁡2t=e−t(cos⁡2t+sin⁡2t)y(t) = e^{-t}\cos 2t + e^{-t}\sin 2t = e^{-t}(\cos 2t + \sin 2t). Without the shift the inverse would be cos⁡2t+sin⁡2t\cos 2t + \sin 2t; the shift multiplies it by e−te^{-t}.

Example
Worked: a repeated factor

Invert Y(s)=1s2(s+1)Y(s) = \dfrac{1}{s^2(s + 1)}.

Step 1 — the factor s2s^2 is repeated, so the form needs both powers: As+Bs2+Cs+1\dfrac As + \dfrac{B}{s^2} + \dfrac{C}{s + 1}.

Step 2 — clear denominators: 1=As(s+1)+B(s+1)+Cs21 = As(s + 1) + B(s + 1) + Cs^2. Setting s=0s = 0 gives B=1B = 1; setting s=−1s = -1 gives C=1C = 1.

Step 3 — compare the s2s^2 coefficients: 0=A+C0 = A + C, so A=−1A = -1.

Step 4 — invert: y(t)=−1+t+e−ty(t) = -1 + t + e^{-t}.

Step 5 — check by transforming back: −1s+1s2+1s+1=−s(s+1)+(s+1)+s2s2(s+1)=1s2(s+1)-\tfrac1s + \tfrac{1}{s^2} + \tfrac{1}{s + 1} = \dfrac{-s(s + 1) + (s + 1) + s^2}{s^2(s + 1)} = \dfrac{1}{s^2(s + 1)} ✓.

Side note
Where it goes wrong: ignoring the linear term

Faced with 1s2+2s+5\dfrac{1}{s^2 + 2s + 5}, it is tempting to match bs2+b2\dfrac{b}{s^2 + b^2} with b2=5b^2 = 5 and answer sin⁡5 t5\dfrac{\sin\sqrt5\,t}{\sqrt5}. But L{15sin⁡5 t}=1s2+5\mathcal L\left\{\tfrac{1}{\sqrt5}\sin\sqrt5\,t\right\} = \dfrac{1}{s^2 + 5}, which has no 2s2s. Completing the square, 1(s+1)2+4=12⋅2(s+1)2+4\dfrac{1}{(s + 1)^2 + 4} = \dfrac12\cdot\dfrac{2}{(s + 1)^2 + 4}, gives 12e−tsin⁡2t\tfrac12e^{-t}\sin 2t. The two are different functions: at t=1t = 1, sin⁡55≈0.35\tfrac{\sin\sqrt5}{\sqrt5} \approx 0.35 while 12e−1sin⁡2≈0.17\tfrac12e^{-1}\sin 2 \approx 0.17.

Try it
Try it: invert by partial fractions

Invert Y(s)=2s−1s2−s−2Y(s) = \dfrac{2s - 1}{s^2 - s - 2}.

Show the answer

Factor: s2−s−2=(s−2)(s+1)s^2 - s - 2 = (s - 2)(s + 1). Write 2s−1(s−2)(s+1)=As−2+Bs+1\dfrac{2s - 1}{(s - 2)(s + 1)} = \dfrac{A}{s - 2} + \dfrac{B}{s + 1}, so 2s−1=A(s+1)+B(s−2)2s - 1 = A(s + 1) + B(s - 2). At s=2s = 2: 3=3A3 = 3A, so A=1A = 1. At s=−1s = -1: −3=−3B-3 = -3B, so B=1B = 1. Hence y(t)=e2t+e−ty(t) = e^{2t} + e^{-t}. Check: L{e2t+e−t}=(s+1)+(s−2)(s−2)(s+1)=2s−1s2−s−2\mathcal L\{e^{2t} + e^{-t}\} = \dfrac{(s + 1) + (s - 2)}{(s - 2)(s + 1)} = \dfrac{2s - 1}{s^2 - s - 2} ✓.

Second-Order Problems and Step Responses

Tip
The pipeline for a second-order IVP

For ay′′+by′+cy=f(t)ay'' + by' + cy = f(t) with y(0)=y0y(0) = y_0 and y′(0)=y1y'(0) = y_1, transforming gives

(as2+bs+c) Y(s)=F(s)+a(sy0+y1)+by0,(as^2 + bs + c)\,Y(s) = F(s) + a(sy_0 + y_1) + by_0,

so Y=F(s)+a(sy0+y1)+by0as2+bs+cY = \dfrac{F(s) + a(sy_0 + y_1) + by_0}{as^2 + bs + c}: the characteristic polynomial sits in the denominator and the initial data in the numerator.

Substitute L{y′′}=s2Y−sy0−y1\mathcal L\{y''\} = s^2Y - sy_0 - y_1 and L{y′}=sY−y0\mathcal L\{y'\} = sY - y_0 and collect the YY terms: a(s2Y−sy0−y1)+b(sY−y0)+cY=Fa(s^2Y - sy_0 - y_1) + b(sY - y_0) + cY = F. The roots of as2+bs+cas^2 + bs + c are the characteristic roots, so the partial-fraction pieces of YY over them reproduce the homogeneous solutions, while the poles of FF produce the particular part. A constant force KK switched on at t=0t = 0 has F=KsF = \tfrac Ks, and its pole at s=0s = 0 is the new equilibrium.

Example
Worked: an unforced second-order IVP

Solve y′′−y′−2y=0y'' - y' - 2y = 0, y(0)=1y(0) = 1, y′(0)=0y'(0) = 0.

Step 1 — transform, with the initial values in place: (s2Y−s)−(sY−1)−2Y=0(s^2Y - s) - (sY - 1) - 2Y = 0.

Step 2 — collect: (s2−s−2)Y=s−1(s^2 - s - 2)Y = s - 1, so Y=s−1(s−2)(s+1)Y = \dfrac{s - 1}{(s - 2)(s + 1)}.

Step 3 — partial fractions: s−1=A(s+1)+B(s−2)s - 1 = A(s + 1) + B(s - 2); s=2s = 2 gives A=13A = \tfrac13 and s=−1s = -1 gives B=23B = \tfrac23.

Step 4 — y(t)=13e2t+23e−ty(t) = \tfrac13e^{2t} + \tfrac23e^{-t}.

Step 5 — check: y(0)=1y(0) = 1 ✓ and y′(0)=23−23=0y'(0) = \tfrac23 - \tfrac23 = 0 ✓, and both exponentials are homogeneous solutions, since 22 and −1-1 are the characteristic roots.

Example
Worked: a constant force switched on

A system at rest obeys y′′+4y′+3y=6y'' + 4y' + 3y = 6 for t≥0t \ge 0, with y(0)=y′(0)=0y(0) = y'(0) = 0.

Step 1 — transform: (s2+4s+3)Y=6s(s^2 + 4s + 3)Y = \dfrac6s, so Y=6s(s+1)(s+3)Y = \dfrac{6}{s(s + 1)(s + 3)}.

Step 2 — partial fractions by covering up each factor and evaluating the rest at its root: at s=0s = 0, 6(1)(3)=2\dfrac{6}{(1)(3)} = 2; at s=−1s = -1, 6(−1)(2)=−3\dfrac{6}{(-1)(2)} = -3; at s=−3s = -3, 6(−3)(−2)=1\dfrac{6}{(-3)(-2)} = 1.

Step 3 — Y=2s−3s+1+1s+3Y = \dfrac2s - \dfrac{3}{s + 1} + \dfrac{1}{s + 3}, so y(t)=2−3e−t+e−3ty(t) = 2 - 3e^{-t} + e^{-3t}.

Step 4 — check: y(0)=2−3+1=0y(0) = 2 - 3 + 1 = 0 ✓, y′(0)=3−3=0y'(0) = 3 - 3 = 0 ✓, and y′′+4y′+3y=(−3e−t+9e−3t)+4(3e−t−3e−3t)+3(2−3e−t+e−3t)=6y'' + 4y' + 3y = (-3e^{-t} + 9e^{-3t}) + 4(3e^{-t} - 3e^{-3t}) + 3(2 - 3e^{-t} + e^{-3t}) = 6 ✓.

Step 5 — the system settles at the new equilibrium y=2y = 2 (where 3y=63y = 6) without overshooting, because both characteristic roots are real.

Example
Worked: a step response that overshoots
ty1234560.51t = π
The response to the constant force: it overshoots the final value 11 (dashed) at t=πt = \pi, then settles.

Solve y′′+2y′+2y=2y'' + 2y' + 2y = 2 for t≥0t \ge 0, from rest.

Step 1 — (s2+2s+2)Y=2s(s^2 + 2s + 2)Y = \dfrac2s, so Y=2s(s2+2s+2)Y = \dfrac{2}{s(s^2 + 2s + 2)}.

Step 2 — partial fractions: 2s(s2+2s+2)=As+Bs+Cs2+2s+2\dfrac{2}{s(s^2 + 2s + 2)} = \dfrac As + \dfrac{Bs + C}{s^2 + 2s + 2} with A=22=1A = \tfrac22 = 1 (cover-up at s=0s = 0); comparing the s2s^2 and ss coefficients in 2=A(s2+2s+2)+(Bs+C)s2 = A(s^2 + 2s + 2) + (Bs + C)s gives B=−1B = -1 and C=−2C = -2.

Step 3 — complete the square: s+2(s+1)2+1=s+1(s+1)2+1+1(s+1)2+1\dfrac{s + 2}{(s + 1)^2 + 1} = \dfrac{s + 1}{(s + 1)^2 + 1} + \dfrac{1}{(s + 1)^2 + 1}.

Step 4 — invert with the shift rule: y(t)=1−e−t(cos⁡t+sin⁡t)y(t) = 1 - e^{-t}(\cos t + \sin t).

Step 5 — the overshoot: y′(t)=2e−tsin⁡ty'(t) = 2e^{-t}\sin t first vanishes again at t=πt = \pi, where y(π)=1+e−π≈1.043y(\pi) = 1 + e^{-\pi} \approx 1.043, about 4%4\% above the final value 11. The complex roots −1±i-1 \pm i are what make the response overshoot.

Side note
Where it goes wrong: dropping the initial values

For y′′+y=0y'' + y = 0, y(0)=1y(0) = 1, y′(0)=0y'(0) = 0, writing L{y′′}=s2Y\mathcal L\{y''\} = s^2Y gives (s2+1)Y=0(s^2 + 1)Y = 0, so Y=0Y = 0 and y≡0y \equiv 0, which has y(0)=0y(0) = 0, not 11. The correct transform s2Y−s⋅1−0s^2Y - s\cdot 1 - 0 gives (s2+1)Y=s(s^2 + 1)Y = s, Y=ss2+1Y = \dfrac{s}{s^2 + 1} and y=cos⁡ty = \cos t. In the transform method the initial values are the only place the data enter; drop them and the solution is lost.

Try it
Try it: a second-order IVP

Solve y′′+5y′+6y=0y'' + 5y' + 6y = 0, y(0)=1y(0) = 1, y′(0)=0y'(0) = 0, by Laplace transform.

Show the answer

Transform: (s2Y−s)+5(sY−1)+6Y=0(s^2Y - s) + 5(sY - 1) + 6Y = 0, so (s2+5s+6)Y=s+5(s^2 + 5s + 6)Y = s + 5 and Y=s+5(s+2)(s+3)Y = \dfrac{s + 5}{(s + 2)(s + 3)}. Cover-up: at s=−2s = -2, 31=3\dfrac{3}{1} = 3; at s=−3s = -3, 2−1=−2\dfrac{2}{-1} = -2. So Y=3s+2−2s+3Y = \dfrac{3}{s + 2} - \dfrac{2}{s + 3} and y=3e−2t−2e−3ty = 3e^{-2t} - 2e^{-3t}. Check: y(0)=1y(0) = 1 ✓ and y′(0)=−6+6=0y'(0) = -6 + 6 = 0 ✓.

Proofs & Why It Matters

Tip
Proof: Le^at = 1/(s−a)

Directly: ∫0∞eate−stdt=∫0∞e−(s−a)tdt=[e−(s−a)t−(s−a)]0∞=1s−a\int_0^\infty e^{at}e^{-st}dt = \int_0^\infty e^{-(s-a)t}dt = \left[\dfrac{e^{-(s-a)t}}{-(s-a)}\right]_0^\infty = \dfrac{1}{s-a} for s>as > a (the boundary term at ∞\infty dies because s−a>0s - a > 0). Setting a=0a = 0 gives L{1}=1s\mathcal L\{1\} = \tfrac1s as a corollary. ■\blacksquare

Tip
Proof: Ly′ = sLy − y(0)

Integrate by parts: ∫0∞y′(t)e−stdt=[y(t)e−st]0∞+s∫0∞y(t)e−stdt=−y(0)+s L{y}\int_0^\infty y'(t)e^{-st}dt = \left[y(t)e^{-st}\right]_0^\infty + s\int_0^\infty y(t)e^{-st}dt = -y(0) + s\,\mathcal L\{y\}, assuming yy grows slower than este^{st} so the boundary term at infinity vanishes.

Iterating gives L{y′′}=s2L{y}−sy(0)−y′(0)\mathcal L\{y''\} = s^2\mathcal L\{y\} - sy(0) - y'(0) — initial conditions enter automatically, which is the method's superpower. ■\blacksquare

Tip
Proof: Ltⁿ = n!/s^n+1

Induction with parts: L{tn}=∫0∞tne−stdt=ns∫0∞tn−1e−stdt=ns L{tn−1}\mathcal L\{t^n\} = \int_0^\infty t^ne^{-st}dt = \tfrac ns\int_0^\infty t^{n-1}e^{-st}dt = \tfrac ns\,\mathcal L\{t^{n-1}\} (boundary terms vanish). Unwinding from L{1}=1s\mathcal L\{1\} = \tfrac1s: L{tn}=ns⋅n−1s⋯1s⋅1s=n!sn+1\mathcal L\{t^n\} = \tfrac{n}{s}\cdot\tfrac{n-1}{s}\cdots\tfrac{1}{s}\cdot\tfrac1s = \tfrac{n!}{s^{n+1}}. ■\blacksquare

Going Deeper: Explanations & Worked Problems

Concept
The transform pipeline, and where each rule enters

Solving an IVP by Laplace is a fixed pipeline. (1) Transform the equation term by term: L{y′}=sY−y(0)\mathcal L\{y'\} = sY - y(0) and L{y′′}=s2Y−sy(0)−y′(0)\mathcal L\{y''\} = s^2Y - sy(0) - y'(0) fold the initial conditions in at the start — no constants to find later.

(2) Solve the resulting algebraic equation for Y(s)Y(s) — always possible, it is linear in YY. (3) Partial-fraction YY into table-sized pieces: 1(s−a)(s−b)=1a−b(1s−a−1s−b)\tfrac{1}{(s-a)(s-b)} = \tfrac{1}{a-b}\left(\tfrac{1}{s-a} - \tfrac{1}{s-b}\right), quadratic pieces aimed at the sine/cosine entries. (4) Invert with the table. The bookkeeping wins where classical methods strain: piecewise forcing (step functions), impulses, and long chains of initial conditions all ride through the algebra unchanged. The cost is fluency with a small table — and partial fractions, which do the heavy lifting in every problem.

Example
Worked: an IVP through the full pipeline

Solve y′+y=1y' + y = 1, y(0)=0y(0) = 0.

Step 1 — transform both sides: (sY−y(0))+Y=1s(sY - y(0)) + Y = \tfrac1s, so (s+1)Y=1s(s + 1)Y = \tfrac1s.

Step 2 — algebra: Y=1s(s+1)Y = \dfrac{1}{s(s+1)}.

Step 3 — partial fractions: write 1s(s+1)=As+Bs+1\dfrac{1}{s(s+1)} = \dfrac{A}{s} + \dfrac{B}{s+1}; multiplying out, 1=A(s+1)+Bs1 = A(s+1) + Bs; set s=0s = 0: A=1A = 1; set s=−1s = -1: B=−1B = -1.

So Y=1s−1s+1Y = \tfrac1s - \tfrac{1}{s+1}.

Step 4 — invert with the table: y(t)=1−e−ty(t) = 1 - e^{-t}.

Step 5 — verify: y′=e−ty' = e^{-t} and y′+y=e−t+1−e−t=1y' + y = e^{-t} + 1 - e^{-t} = 1 ✓, y(0)=0y(0) = 0 ✓. Compare with the integrating-factor route: same answer, but here the initial condition entered at Step 1 and no constant CC ever appeared.

Example
Worked: a second-order IVP with a sine, by table and partial fractions

Solve y′′+y=0y'' + y = 0, y(0)=0y(0) = 0, y′(0)=3y'(0) = 3 — deliberately simple, to watch the machinery.

Step 1 — transform: s2Y−s⋅0−3+Y=0s^2Y - s\cdot0 - 3 + Y = 0.

Step 2 — solve: Y=3s2+1Y = \dfrac{3}{s^2 + 1}.

Step 3 — no partial fractions needed: this is 33 times the table entry L{sin⁡t}=1s2+1\mathcal L\{\sin t\} = \tfrac{1}{s^2+1}.

Step 4 — invert: y=3sin⁡ty = 3\sin t.

Step 5 — verify: y′′=−3sin⁡t=−yy'' = -3\sin t = -y ✓, y(0)=0y(0) = 0 ✓, y′(0)=3cos⁡0=3y'(0) = 3\cos 0 = 3 ✓. Now the same skeleton with forcing, damping, or a step function only lengthens

Step 3 — the pipeline itself never changes shape, which is exactly its virtue.