How to solve initial value problems with Laplace transforms

Laplace transforms turn the derivatives in an initial value problem into expressions involving one unknown, Y(s). We can solve for that unknown with algebra, then take the inverse transform to find y(t). Let’s work through the steps on a second-order equation.

Transform the derivatives and use the initial values

Find y(t) for the following equation and initial conditions.

y′′-10y′+9y=5ty(0)=-1,y′(0)=2

Write Y(s)=ℒ(y(t)). The derivative transforms include the initial values, so we can substitute them as we transform each term.

ℒ(y′)=sY-y(0)ℒ(y′′)=s2Y-sy(0)-y′(0)ℒ(5t)=5s2

Here the first derivative transforms to sY+1 and the second transforms to s2Y+s-2. Bring the transformed terms together.

(s2Y+s-2)-10(sY+1)+9Y=5s2

Solve for the transformed function

Collect the terms containing the unknown transform, then factor the quadratic in the denominator.

(s2-10s+9)Y=5s2+12-sY(s)=5+12s2-s3s2(s-9)(s-1)

Decompose into partial fractions

The factor s occurs twice, so we need terms with denominators s and s2. The other factors each contribute one term.

Y(s)=As+Bs2+Cs-9+Ds-1

Multiply by s2(s-9)(s-1). Setting s=0 isolates B. Setting s=9 isolates C, and setting s=1 isolates D.

B=5(-9)(-1)=59C=24892(9-1)=3181D=1612(1-9)=-2

The coefficient of s3 in the original numerator is -1. In the combined fractions, the A, C, and D terms each contribute to that coefficient.

A+C+D=-1A=-1-3181+2=5081

Take the inverse Laplace transform

The inverse transforms of 1/s, 1/s2, and 1/(s-a) are 1, t, and eat, respectively. Apply those pairs to the four terms.

y(t)=5081+59t+3181e9t-2et

Check the initial values and the equation

y(0)=50+3181-2=-1y′(0)=59+319-2=2

The exponentials satisfy the homogeneous equation because 9 and 1 are roots of r2-10r+9. Applying y′′-10y′+9y to the remaining terms 50/81+5t/9 gives -50/9+50/9+5t=5t. This checks the forcing term.

Shift the time variable for initial values at another point

For initial values at t=t0, set τ=t-t0 and write u(τ)=y(τ+t0). The given values now specify u(0) and u′(0). Transform the equation in τ, then substitute τ=t-t0 into the solution.