Chapter 7

Laplace Transform
The Laplace transform can be used to solve differential equations. Besides being a different and efficient alternative to variation of parameters and undetermined coefficients, the Laplace method is particularly
advantageous for input terms that are piecewise-defined, periodic or impulsive.
The direct Laplace transform or the Laplace integral of a function
f (t) defined for 0 t < is the ordinary calculus integration problem
Z

f (t)est dt,

succinctly denoted L(f (t)) in science and engineering literature. The

Lnotation recognizes that integration always proceeds over t = 0 to
t = and that the integral involves an integrator est dt instead of the
usual dt. These minor differences distinguish Laplace integrals from
the ordinary integrals found on the inside covers of calculus texts.

7.1 Introduction to the Laplace Method

The foundation of Laplace theory is Lerchs cancellation law
R
0

(1)

y(t)est dt =

L(y(t) = L(f (t))

f (t)est dt

implies
or
implies

y(t) = f (t),
y(t) = f (t).

In differential equation applications, y(t) is the sought-after unknown

while f (t) is an explicit expression taken from integral tables.
Below, we illustrate Laplaces method by solving the initial value problem
y 0 = 1, y(0) = 0.
The method obtains a relation L(y(t)) = L(t), whence Lerchs cancellation law implies the solution is y(t) = t.
The Laplace method is advertised as a table lookup method, in which
the solution y(t) to a differential equation is found by looking up the
answer in a special integral table.

7.1 Introduction to the Laplace Method

247

g(t)est dt is called the Laplace

R
integral of the function g(t). It is defined by limN 0N g(t)est dt and
depends on variable s. The ideas will be illustrated for g(t) = 1, g(t) = t
and g(t) = t2 , producing the integral formulas in Table 1.

Laplace Integral. The integral

R
0

R
0

= 1/s

(t)est dt =

=
0

R
0

(t2 )est dt

=
=
=

R
0

Laplace integral of g(t) = t.

(1)est dt

Use

d
st )dt
0 ds (te
d R
st dt
ds
0 (t)e
d
(1/s2 )
ds
2/s3

1
s

d
ds F (t, s)dt

d
ds

F (t, s)dt.

Differentiate.

(1)est dt =

Use L(1) = 1/s.

Table 1. The Laplace integral

Assumed s > 0.

d
(1/s)
ds
1/s2

Laplace integral of g(t) = 1.

d
(est )dt
ds

d
= ds

=
R

t=

(1)est dt = (1/s)est t=0

R
0

R
0

Laplace integral of g(t) = t2 .

Use L(t) = 1/s2 .

g(t)est dt for g(t) = 1, t and t2 .

(t)est dt =

1
s2

In summary, L(tn ) =

R
0

n!

(t2 )est dt =

2
s3

s1+n

An Illustration. The ideas of the Laplace method will be illustrated for the solution y(t) = t of the problem y 0 = 1, y(0) = 0. The
method, entirely different from variation of parameters or undetermined
coefficients, uses basic calculus and college algebra; see Table 2.
Table 2. Laplace method details for the illustration y 0 = 1, y(0) = 0.

y 0 (t)est = est

Multiply y 0 = 1 by est .

Integrate t = 0 to t = .

R
0
st dt = est dt
0 y (t)e
0
R 0
st dt = 1/s
y
(t)e
0
R
s 0 y(t)est dt y(0) = 1/s
R
st dt = 1/s2
0 y(t)e
R
R
st dt = (t)est dt
0 y(t)e
0

y(t) = t

Use Table 1.
Integrate by parts on the left.
Use y(0) = 0 and divide.
Use Table 1.
Apply Lerchs cancellation law.

248

Laplace Transform

In Lerchs law, the formal rule of erasing the integral signs is valid provided the integrals are equal for large s and certain conditions hold on y
and f see Theorem 2. The illustration in Table 2 shows that Laplace
theory requires an in-depth study of a special integral table, a table
which is a true extension of the usual table found on the inside covers
of calculus books. Some entries for the special integral table appear in
Table 1 and also in section 7.2, Table 4.
The L-notation for the direct Laplace transform produces briefer details,
as witnessed by the translation of Table 2 into Table 3 below. The reader
is advised to move from Laplace integral notation to the Lnotation as
soon as possible, in order to clarify the ideas of the transform method.
Table 3. Laplace method L-notation details for y 0 = 1, y(0) = 0
translated from Table 2.

L(y 0 (t)) = L(1)

Apply L across y 0 = 1, or multiply y 0 =

1 by est , integrate t = 0 to t = .

L(y 0 (t)) = 1/s

Use Table 1.

sL(y(t)) y(0) = 1/s

Integrate by parts on the left.

L(y(t)) =

1/s2

Use y(0) = 0 and divide.

L(y(t)) = L(t)

Apply Table 1.

y(t) = t

Invoke Lerchs cancellation law.

Some Transform Rules. The formal properties of calculus integrals

plus the integration by parts formula used in Tables 2 and 3 leads to these
rules for the Laplace transform:
L(f (t) + g(t)) = L(f (t)) + L(g(t))

The integral of a sum is the

sum of the integrals.

L(cf (t)) = cL(f (t))

Constants c pass through the

integral sign.

L(y 0 (t)) = sL(y(t)) y(0)

The t-derivative rule, or integration by parts. See Theorem 3.

Lerchs cancellation law. See
Theorem 2.

L(y(t)) = L(f (t)) implies y(t) = f (t)

1 Example (Laplace method) Solve by Laplaces method the initial value

problem y 0 = 5 2t, y(0) = 1.
Solution: Laplaces method is outlined in Tables 2 and 3. The L-notation of
Table 3 will be used to find the solution y(t) = 1 + 5t t2 .

7.1 Introduction to the Laplace Method

L(y 0 (t)) = L(5 2t)
2
5
L(y 0 (t)) = 2
s s
5
2
sL(y(t)) y(0) = 2
s s
1
5
2
L(y(t)) = + 2 3
s s
s
L(y(t)) = L(1) + 5L(t) L(t2 )
2

= L(1 + 5t t )

y(t) = 1 + 5t t

249

Apply L across y 0 = 5 2t.

Use Table 1.
Apply the t-derivative rule, page 248.
Use y(0) = 1 and divide.
Apply Table 1, backwards.
Linearity, page 248.
Invoke Lerchs cancellation law.

2 Example (Laplace method) Solve by Laplaces method the initial value

problem y 00 = 10, y(0) = y 0 (0) = 0.
Solution: The L-notation of Table 3 will be used to find the solution y(t) = 5t2 .
L(y 00 (t)) = L(10)

Apply L across y 00 = 10.

sL(y 0 (t)) y 0 (0) = L(10)

s[sL(y(t)) y(0)] y 0 (0) = L(10)

s2 L(y(t)) = L(10)
10
L(y(t)) = 3
s
L(y(t)) = L(5t2 )

Apply the t-derivative rule to y 0 , that is,

replace y by y 0 on page 248.
Repeat the t-derivative rule, on y.
Use y(0) = y 0 (0) = 0.
Use Table 1. Then divide.
Apply Table 1, backwards.

y(t) = 5t2

Invoke Lerchs cancellation law.

est f (t) dt
is known to exist in the sense of the improper integral definition1

Existence of the Transform. The Laplace integral

Z

g(t)dt = lim

N 0

g(t)dt

provided f (t) belongs to a class of functions known in the literature as

functions of exponential order. For this class of functions the relation
(2)

f (t)
=0
t eat
lim

is required to hold for some real number a, or equivalently, for some

constants M and ,
(3)
|f (t)| M et .
In addition, f (t) is required to be piecewise continuous on each finite
subinterval of 0 t < , a term defined as follows.
1
An advanced calculus background is assumed for the Laplace transform existence
proof. Applications of Laplace theory require only a calculus background.

250

Laplace Transform

Definition 1 (piecewise continuous)

A function f (t) is piecewise continuous on a finite interval [a, b] provided there exists a partition a = t0 < < tn = b of the interval [a, b]
and functions f1 , f2 , . . . , fn continuous on (, ) such that for t not
a partition point

(4)

f1 (t)

..
f (t) =
.

fn (t)

t0

< t < t1 ,
..
.

tn1 < t < tn .

The values of f at partition points are undecided by equation (4). In

particular, equation (4) implies that f (t) has one-sided limits at each
point of a < t < b and appropriate one-sided limits at the endpoints.
Therefore, f has at worst a jump discontinuity at each partition point.
3 Example (Exponential order) Show that f (t) = et cos t + t is of exponential order, that is, show that f (t) is piecewise continuous and find > 0
such that limt f (t)/et = 0.
Solution: Already, f (t) is continuous, hence piecewise continuous. From
LHospitals rule in calculus, limt p(t)/et = 0 for any polynomial p and
any > 0. Choose = 2, then
lim

cos t
t
f (t)
= lim
+ lim 2t = 0.
2t
t
t
t
e
e
e

Theorem 1 (Existence of L(f ))

Let f (t) be piecewise continuous on every finite interval in t 0 and satisfy
|f (t)| M et for some constants M and . Then L(f (t)) exists for s >
and lims L(f (t)) = 0.
Proof: It has to be shown that the Laplace integral of f is finite for s > .
Advanced calculus implies that it is sufficient to show that the integrand is absolutely bounded above by an integrable function g(t). Take g(t) = M e (s)t .
Then g(t) 0. Furthermore, g is integrable, because
Z
M
g(t)dt =
.
s

0
Inequality |f (t)| M et implies the absolute value of the Laplace transform
integrand f (t)est is estimated by


f (t)est  M et est = g(t).

M
, because the
s
right side of this inequality has limit zero at s = . The proof is complete.
The limit statement follows from |L(f (t))|

R
0

g(t)dt =

7.1 Introduction to the Laplace Method

251

Theorem 2 (Lerch)
R
If f (t) and f2 (t) are continuous, of exponential order and 0 f1 (t)est dt =
R 1
st dt for all s > s , then f (t) = f (t) for t 0.
0
1
2
0 f2 (t)e

Proof: See Widder [?].

Theorem 3 (t-Derivative Rule)

If f (t) is continuous, lim f (t)est = 0 for all large values of s and f 0 (t)
t

is piecewise continuous, then L(f 0 (t)) exists for all large s and L(f 0 (t)) =
sL(f (t)) f (0).
Proof: See page 276.

Exercises 7.1
PN
Solve the given 18. f (t) =
n=1 cn sin(nt), for any
choice of the constants c1 , . . . , cN .
initial value problem using Laplaces
method.
Existence of transforms. Let f (t) =
1. y 0 = 2, y(0) = 0.
2
2
tet sin(et ). Establish these results.
2. y 0 = 1, y(0) = 0.
19. The function f (t) is not of expo-

Laplace method.

nential order.

3. y 0 = t, y(0) = 0.

20. The
integral of f (t),
R Laplace
st
f
(t)e
dt,
converges for all
0
s > 0.

4. y 0 = t, y(0) = 0.
5. y 0 = 1 t, y(0) = 0.
6. y 0 = 1 + t, y(0) = 0.
7. y 0 = 3 2t, y(0) = 0.
8. y 0 = 3 + 2t, y(0) = 0.
00

9. y = 2, y(0) = y (0) = 0.
00

10. y = 1, y(0) = y (0) = 0.

00

11. y = 1 t, y(0) = y (0) = 0.

12. y 00 = 1 + t, y(0) = y 0 (0) = 0.
00

00

13. y = 3 2t, y(0) = y (0) = 0.

14. y = 3 + 2t, y(0) = y (0) = 0.

Jump Magnitude. For f piecewise

continuous, define the jump at t by
J(t) = lim f (t + h) lim f (t h).
h0+

h0+

Compute J(t) for the following f .

21. f (t) = 1 for t 0, else f (t) = 0
22. f (t) = 1 for t 1/2, else f (t) = 0
23. f (t) = t/|t| for t 6= 0, f (0) = 0
24. f (t) = sin t/| sin t| for t 6= n,
f (n) = (1)n

.
P series
Exponential order. Show that f (t) Taylor
n

The series relation

P

n
L( n=0 cn t ) =
n=0 cn L(t ) often
is of exponential order, by finding a
holds, in which case the result L(tn ) =
constant 0 in each case such that
n!s1n can be employed to find a
f (t)
= 0.
lim
series representation of the Laplace
t et
transform. Use this idea on the fol15. f (t) = 1 + t
lowing to find a series formula for
L(f (t)).
16. f (t) = et sin(t)
P
2t
n
PN
17. f (t) = n=0 cn xn , for any choice 25. f (t) = e = n=0 (2t) /n!
P
of the constants c0 , . . . , cN .
26. f (t) = et = n=0 (t)n /n!

252

Laplace Transform

7.2 Laplace Integral Table

The objective in developing a table of Laplace integrals, e.g., Tables 4
and 5, is to keep the table size small. Table manipulation rules appearing in Table 6, page 257, effectively increase the table size manyfold,
making it possible to solve typical differential equations from electrical
and mechanical problems. The combination of Laplace tables plus the
table manipulation rules is called the Laplace transform calculus.
Table 4 is considered to be a table of minimum size to be memorized.
Table 5 adds a number of special-use entries. For instance, the Heaviside
entry in Table 5 is memorized, but usually not the others.
Derivations are postponed to page 270. The theory of the gamma function (x) appears below on page 255.
Table 4. A minimal Laplace integral table with L-notation
R
0

R
0

(tn )est dt =
(eat )est dt =

n!
1
sa

(cos bt)est dt =

s1+n

L(eat ) =

s
+ b2
R
b
st
dt = 2
0 (sin bt)e
s + b2
R

n!

L(tn ) =

s1+n

1
sa

s
+ b2
b
L(sin bt) = 2
s + b2
L(cos bt) =

s2

s2

Table 5. Laplace integral table extension

L(H(t a)) =

eas
(a 0)
s

L((t a)) = eas

L(floor(t/a)) =

eas
s(1 eas )

L(sqw(t/a)) =

1
tanh(as/2)
s

L(a trw(t/a)) =

1
tanh(as/2)
s2

(1 + )
s1+
r

1/2
L(t
)=
s

L(t ) =

Heaviside
 unit step, defined by
1
for t 0,
H(t) =
0
otherwise.
Dirac delta, (t) = dH(t).
Special usage rules apply.
Staircase function,
floor(x) = greatest integer x.
Square wave,
sqw(x) = (1)floor(x) .
TriangularRwave,
x
trw(x) = 0 sqw(r)dr.

Generalized power
R function,
(1 + ) = 0 ex x dx.

Because (1/2) =

7.2 Laplace Integral Table

253

4 Example (Laplace transform) Let f (t) = t(t 1) sin 2t + e3t . Compute

L(f (t)) using the basic Laplace table and transform linearity properties.
Solution:
L(f (t)) = L(t2 5t sin 2t + e3t )
2

Expand t(t 5).

3t

= L(t ) 5L(t) L(sin 2t) + L(e )

2
5
2
1
= 3 2 2
+
s
s
s +4 s3

Linearity applied.
Table lookup.

5 Example (Inverse Laplace transform) Use the basic Laplace table backwards plus transform linearity properties to solve for f (t) in the equation
L(f (t)) =

s2

s
2
s+1
+
+ 3 .
+ 16 s 3
s

Solution:
1
1
1 2
s
+2
+ 2+
+ 16
s3 s
2 s3
3t
= L(cos 4t) + 2L(e ) + L(t) + 12 L(t2 )

L(f (t)) =

s2

3t

= L(cos 4t + 2e + t +
3t

f (t) = cos 4t + 2e + t +

1 2
2t )

1 2
2t

Convert to table entries.

Laplace table (backwards).
Linearity applied.
Lerchs cancellation law.

6 Example (Heaviside) Find the Laplace transform of f (t) in Figure 1.

5

1
1

Figure 1. A piecewise defined

function f (t) on 0 t < : f (t) = 0
except for 1 t < 2 and 3 t < 4.

Solution: The details require the use of the Heaviside function formula
H(t a) H(t b) =
The formula for f (t):


1 1 t < 2,
1
5
3

t
<
4,
f (t) =
=
0

0 otherwise

1 a t < b,
0 otherwise.

1 t < 2,
+5
otherwise

1
0

3 t < 4,
otherwise

Then f (t) = f1 (t) + 5f2 (t) where f1 (t) = H(t 1) H(t 2) and f2 (t) =
H(t 3) H(t 4). The extended table gives
L(f (t)) = L(f1 (t)) + 5L(f2 (t))

= L(H(t 1)) L(H(t 2)) + 5L(f2 (t))

Linearity.
Substitute for f1 .

254

Laplace Transform
es e2s
+ 5L(f2 (t))
s
es e2s + 5e3s 5e4s
=
s

Extended table used.

Similarly for f2 .

7 Example (Dirac delta) A machine shop tool that repeatedly hammers a

P
die is modeled by the Dirac impulse model f (t) = N
n=1 (t n). Show
PN
ns
that L(f (t)) = n=1 e .
Solution:

P

N
L(f (t)) = L
(t

n)
n=1
PN
= n=1 L((t n))
P
ns
= N
n=1 e

Linearity.
Extended Laplace table.

8 Example (Square wave) A periodic camshaft force f (t) applied to a mechanical system has the idealized graph shown in Figure 2. Show that
f (t) = 1 + sqw(t) and L(f (t)) = 1s (1 + tanh(s/2)).
2

Figure 2. A periodic force f (t) applied

to a mechanical system.

Solution:
1 + sqw(t)

1+1
11

2n t < 2n + 1, n = 0, 1, . . .,
2n + 1 t < 2n + 2, n = 0, 1, . . .,

2
0
= f (t).

2n t < 2n + 1, n = 0, 1, . . .,
otherwise,

By the extended Laplace table, L(f (t)) = L(1) + L(sqw(t)) =

1 tanh(s/2)
+
.
s
s

9 Example (Sawtooth wave) Express the P -periodic sawtooth wave represented in Figure 3 as f (t) = ct/P c floor(t/P ) and obtain the formula
L(f (t)) =

ceP s
c

.
P s2 s seP s

4P

Figure 3. A P -periodic sawtooth

wave f (t) of height c > 0.

7.2 Laplace Integral Table

255

Solution: The representation originates from geometry, because the periodic

function f can be viewed as derived from ct/P by subtracting the correct constant from each of intervals [P, 2P ], [2P, 3P ], etc.
The technique used to verify the identity is to define g(t) = ct/P c floor(t/P )
and then show that g is P -periodic and f (t) = g(t) on 0 t < P . Two P periodic functions equal on the base interval 0 t < P have to be identical,
hence the representation follows.
The fine details: for 0 t < P , floor(t/P ) = 0 and floor(t/P + k) = k. Hence
g(t + kP ) = ct/P + ck c floor(k) = ct/P = g(t), which implies that g is
P -periodic and g(t) = f (t) for 0 t < P .
c
L(t) cL(floor(t/P ))
P
c
ceP s
=

P s2
s seP s

L(f (t)) =

Linearity.
Basic and extended table applied.

10 Example (Triangular wave) Express the triangular wave f of Figure 4 in

5
terms of the square wave sqw and obtain L(f (t)) = 2 tanh(s/2).
s
5

Figure 4. A 2-periodic triangular

wave f (t) of height 5.

Solution: The representation of f in terms of sqw is f (t) = 5

R t/
0

sqw(x)dx.

Details: A 2-periodic triangular wave of height 1 is obtained by integrating

the square wave of period 2. A wave of height c and period 2 is given by
Rt
R 2t/P
c trw(t) = c 0 sqw(x)dx. Then f (t) = c trw(2t/P ) = c 0
sqw(x)dx where
c = 5 and P = 2.

Laplace transform details: Use the extended Laplace table as follows.

L(f (t)) =

5
5
L( trw(t/)) = 2 tanh(s/2).

Gamma Function. In mathematical physics, the Gamma function or the generalized factorial function is given by the identity
(1)

(x) =

et tx1 dt,

x > 0.

This function is tabulated and available in computer languages like Fortran, C and C++. It is also available in computer algebra systems and
numerical laboratories. Some useful properties of (x):
(2)

(1 + x) = x(x)

(3)

(1 + n) = n! for integers

n 1.

256

Laplace Transform
R

et dt = 1, which gives
(1) = 1. Use this identity and successively relation (2) to obtain relation (3).
To prove identity (2), integration by parts is applied, as follows:
R
(1 + x) = 0 et tx dt
Definition.
R
t=

= tx et |t=0 + 0 et xtx1 dt
Use u = tx , dv = et dt.
R t x1
=x 0 e t
dt
Boundary terms are zero
for x > 0.
= x(x).

Details for relations (2) and (3): Start with

Exercises 7.2
Laplace

transform.
Compute Inverse Laplace transform. Solve
L(f (t)) using the basic Laplace table the given equation for the function
and the linearity properties of the f (t). Use the basic table and linearity
transform. Do not use the direct properties of the Laplace transform.
Laplace transform!
21. L(f (t)) = s2
1. L(2t)
22. L(f (t)) = 4s2
2. L(4t)
23. L(f (t)) = 1/s + 2/s2 + 3/s3
2
3. L(1 + 2t + t )
24. L(f (t)) = 1/s3 + 1/s
2
4. L(t 3t + 10)
25. L(f (t)) = 2/(s2 + 4)
5. L(sin 2t)
26. L(f (t)) = s/(s2 + 4)
6. L(cos 2t)
27. L(f (t)) = 1/(s 3)
7. L(e2t )
28. L(f (t)) = 1/(s + 3)
8. L(e2t )
29. L(f (t)) = 1/s + s/(s2 + 4)
9. L(t + sin 2t)
30. L(f (t)) = 2/s 2/(s2 + 4)
10. L(t cos 2t)
31. L(f (t)) = 1/s + 1/(s 3)
11. L(t + e2t )
32. L(f (t)) = 1/s 3/(s 2)
12. L(t 3e2t )
33. L(f (t)) = (2 + s)2 /s3
13. L((t + 1)2 )
34. L(f (t)) = (s + 1)/s2
14. L((t + 2)2 )
35. L(f (t)) = s(1/s2 + 2/s3 )
15. L(t(t + 1))
36. L(f (t)) = (s + 1)(s 1)/s3
16. L((t + 1)(t + 2))
P10
37. L(f (t)) = n=0 n!/s1+n
P10 n
17. L( n=0 t /n!)
P10
38. L(f (t)) = n=0 n!/s2+n
P10 n+1
18. L( n=0 t
/n!)
P10
n
39. L(f (t)) = n=1 2
P10
s + n2
19. L( n=1 sin nt)
P10
s
P10
40. L(f (t)) = n=0 2
20. L( n=0 cos nt)
s + n2

7.3 Laplace Transform Rules

257

7.3 Laplace Transform Rules

In Table 6, the basic table manipulation rules are summarized. Full
statements and proofs of the rules appear in section 7.7, page 275.
The rules are applied here to several key examples. Partial fraction
expansions do not appear here, but in section 7.4, in connection with
Heavisides coverup method.
Table 6. Laplace transform rules
L(f (t) + g(t)) = L(f (t)) + L(g(t))

Linearity.

L(cf (t)) = cL(f (t))

Linearity.

L(y 0 (t)) = sL(y(t)) y(0)

The t-derivative rule.

 1
L
g(x)dx = L(g(t))
s
d
L(tf (t)) = L(f (t))
ds
R

t
0

The Laplace of a sum is the sum of the Laplaces.

Constants move through the L-symbol.

Derivatives L(y 0 ) are replaced in transformed equations.

The t-integral rule.

The s-differentiation rule.
Multiplying f by t applies d/ds to the transform of f .

L(eat f (t)) = L(f (t))|s(sa)

First shifting rule.

L(f (t a)H(t a)) = eas L(f (t)),

L(g(t)H(t a)) = eas L(g(t + a))
RP
f (t)est dt
L(f (t)) = 0
1 eP s

Second shifting rule.

L(f (t))L(g(t)) = L((f g)(t))

Convolution rule.
Rt

Multiplying f by eat replaces s by s a.

First and second forms.

Rule for P -periodic functions.

Assumed here is f (t + P ) = f (t).

Define (f g)(t) =

f (x)g(t x)dx.

11 Example (Harmonic oscillator) Solve by Laplaces method the initial value

problem x00 + x = 0, x(0) = 0, x0 (0) = 1.
Solution: The solution is x(t) = sin t. The details:
L(x00 ) + L(x) = L(0)
0

sL(x ) x (0) + L(x) = 0

0

s[sL(x) x(0)] x (0) + L(x) = 0

2

(s + 1)L(x) = 1
1
L(x) = 2
s +1
= L(sin t)
x(t) = sin t

Apply L across the equation.

Use the t-derivative rule.

Use again the t-derivative rule.

Use x(0) = 0, x0 (0) = 1.
Divide.
Basic Laplace table.
Invoke Lerchs cancellation law.

258

Laplace Transform
2
.
(s 5)3

12 Example (s-differentiation rule) Show the steps for L(t2 e5t ) =

Solution:




d
L(t e ) =

L(e5t )
ds


1
2 d d
= (1)
ds ds s 5


d
1
=
ds (s 5)2
2
=
(s 5)3
2 5t

ds

Apply s-differentiation.
Basic Laplace table.
Calculus power rule.
Identity verified.

13 Example (First shifting rule) Show the steps for L(t2 e3t ) =

2
.
(s + 3)3

Solution:

L(t2 e3t ) = L(t2 )ss(3)



2

=

2+1
s
ss(3)
=

First shifting rule.

Basic Laplace table.

2
(s + 3)3

Identity verified.

14 Example (Second shifting rule) Show the steps for

L(sin t H(t )) =

es
.
s2 + 1

Solution: The second shifting rule is applied as follows.

L(sin t H(t )) = L(g(t)H(t a)
= eas L(g(t + a)
=e

=e

= es

Choose g(t) = sin t, a = .

Second form, second shifting theorem.

L(sin(t + )) Substitute a = .

L( sin t)
1
s2 + 1

Sum rule sin(a + b) = sin a cos b +

sin b cos a plus sin = 0, cos = 1.
Basic Laplace table. Identity verified.

15 Example (Trigonometric formulas) Show the steps used to obtain these

Laplace identities:
3
2
s 2 a2
2 cos at) = 2(s 3sa )
(a) L(t cos at) = 2
(c)
L(t
(s + a2 )2
(s2 + a2 )3
2
3
2sa
2 sin at) = 6s a a
(d)
L(t
(b) L(t sin at) = 2
(s + a2 )2
(s2 + a2 )3

7.3 Laplace Transform Rules

259

Solution: The details for (a):

L(t cos at) = (d/ds)L(cos at)


s
d
=
ds s2 + a2
=

s2 a 2
(s2 + a2 )2

Use s-differentiation.
Basic Laplace table.
Calculus quotient rule.

The details for (c):

L(t2 cos at) = (d/ds)L((t) cos at)


s2 a 2
d
2
=
ds
(s + a2 )2
=

2s3 6sa2 )
(s2 + a2 )3

Use s-differentiation.
Result of (a).
Calculus quotient rule.

The similar details for (b) and (d) are left as exercises.

16 Example (Exponentials) Show the steps used to obtain these Laplace

identities:
2
2
sa
at cos bt) = (s a) b
(c)
L(te
(a) L(eat cos bt) =
(s a)2 + b2
((s a)2 + b2 )2
2b(s a)
b
(d) L(teat sin bt) =
(b) L(eat sin bt) =
(s a)2 + b2
((s a)2 + b2 )2
Solution: Details for (a):
L(eat cos bt) = L(cos bt)|ssa



s

=
s2 + b2 ssa
sa
=
(s a)2 + b2

First shifting rule.

Basic Laplace table.
Verified (a).

Details for (c):

L(teat cos bt) = L(t cos bt)|ssa



d
= L(cos bt) 
ds
ssa




s
d

=

2
2
ds s + b
ssa
 2


s b2

=
(s2 + b2 )2 ssa
=

(s a)2 b2
((s a)2 + b2 )2

Left as exercises are (b) and (d).

First shifting rule.

Apply s-differentiation.
Basic Laplace table.
Calculus quotient rule.
Verified (c).

260

Laplace Transform

17 Example (Hyperbolic functions) Establish these Laplace transform facts

about cosh u = (eu + eu )/2 and sinh u = (eu eu )/2.
s 2 + a2
(s2 a2 )2
2as
(d) L(t sinh at) = 2
(s a2 )2

s
2
s a2
a
(b) L(sinh at) = 2
s a2

(c) L(t cosh at) =

(a) L(cosh at) =

Solution: The details for (a):

L(cosh at) = 12 (L(eat ) + L(eat ))


1
1
1
+
=
2 sa s+a
s
= 2
s a2

Definition plus linearity of L.

Basic Laplace table.
Identity (a) verified.

The details for (d):



d
a
L(t sinh at) =
ds s2 a2
a(2s)
= 2
(s a2 )2

Apply the s-differentiation rule.

Calculus power rule; (d) verified.

Left as exercises are (b) and (c).

18 Example (s-differentiation) Solve L(f (t)) =

(s2

2s
for f (t).
+ 1)2

Solution: The solution is f (t) = t sin t. The details:

2s
(s2 + 1)2


d
1
=
ds s2 + 1
d
= (L(sin t))
ds
= L(t sin t)

L(f (t)) =

f (t) = t sin t

Calculus power rule (un )0 = nun1 u0 .

Basic Laplace table.
Apply the s-differentiation rule.
Lerchs cancellation law.

19 Example (First shift rule) Solve L(f (t)) =

s+2
for f (t).
22 + 2s + 2

Solution: The answer is f (t) = et cos t + et sin t. The details:

L(f (t)) =
=

s+2
s2 + 2s + 2
s+2
(s + 1)2 + 1

Signal for this method: the denominator has complex roots.

Complete the square, denominator.

7.3 Laplace Transform Rules

261

S+1
S2 + 1
1
S
+
= 2
S + 1 S2 + 1
= L(cos t) + L(sin t)|sS=s+1
=

= L(e

f (t) = e

cos t) + L(e

cos t + e

sin t)

sin t

Substitute S for s + 1.
Split into Laplace table entries.
Basic Laplace table.
First shift rule.
Invoke Lerchs cancellation law.

20 Example (Damped oscillator) Solve by Laplaces method the initial value

problem x00 + 2x0 + 2x = 0, x(0) = 1, x0 (0) = 1.
Solution: The solution is x(t) = et cos t. The details:
L(x00 ) + 2L(x0 ) + 2L(x) = L(0)

Apply L across the equation.

s[sL(x) x(0)] x0 (0)

+2[L(x) x(0)] + 2L(x) = 0

The t-derivative rule on x.

The t-derivative rule on x0 .

sL(x ) x (0) + 2L(x ) + 2L(x) = 0

(s2 + 2s + 2)L(x) = 1 + s
s+1
L(x) = 2
s + 2s + 2
s+1
=
(s + 1)2 + 1

Use x(0) = 1, x0 (0) = 1.

Divide.
Complete the square in the denominator.
Basic Laplace table.

= L(cos t)|ss+1
= L(et cos t)

First shifting rule.

x(t) = et cos t

Invoke Lerchs cancellation law.

21 Example (Rectified sine wave) Compute the Laplace transform of the

rectified sine wave f (t) = | sin t| and show it can be expressed in the
form
s

coth 2
.
L(| sin t|) =
s2 + 2
Solution: The periodic function formula will be applied
with period P =
R
P

2/. The calculation reduces to the evaluation of J = 0 f (t)est dt. Because

sin t 0 on / t 2/, integral J can be written as J = J1 + J2 , where
J1 =

sin t est dt,

J2 =

2/
/

sin t est dt.

Integral tables give the result

Z
est cos(t) sest sin(t)
sin t est dt =

.
s2 + 2
s2 + 2
Then
J1 =

(es/ + 1)
,
s2 + 2

J2 =

(e2s/ + es/ )
,
s2 + 2

262

Laplace Transform

J=

(es/ + 1)2
.
s2 + 2

The remaining challenge is to write the answer for L(f (t)) in terms of coth.
The details:
L(f (t)) =
=
=

J
1 eP s
(1

Periodic function formula.

J

eP s/2 )(1

Apply 1 x2 = (1 x)(1 + x),

x = eP s/2 .

+ eP s/2 )

(1 + eP s/2 )
(1 eP s/2 )(s2 + 2 )

Cancel factor 1 + eP s/2 .

eP s/4 + eP s/4

eP s/4 eP s/4 s2 + 2
2 cosh(P s/4)

=
2
2 sinh(P s/4) s + 2
coth(P s/4)
=
s2 + 2


s
coth 2
=
s2 + 2

Factor out eP s/4 , then cancel.

Apply cosh, sinh identities.

Use coth u = cosh u/ sinh u.
Identity verified.

22 Example (Halfwave rectification) Compute the Laplace transform of the

halfwave rectification of sin t, denoted g(t), in which the negative cycles
of sin t have been canceled to create g(t). Show in particular that
1
s
L(g(t)) =
1 + coth
2
2
2s +
2




Solution: The halfwave rectification of sin t is g(t) = (sin t + | sin t|)/2.

Therefore, the basic Laplace table plus the result of Example 21 give
L(2g(t)) = L(sin t) + L(| sin t|)

cosh(s/(2))
= 2
+
2
s +
s2 + 2

(1 + cosh(s/(2))
= 2
s + 2
Dividing by 2 produces the identity.

23 Example (Shifting rules) Solve L(f (t)) = e3s

s2

s+1
for f (t).
+ 2s + 2

Solution: The answer is f (t) = e3t cos(t 3)H(t 3). The details:
s+1
(s + 1)2 + 1
S
= e3s 2
S +1
= e3S+3 (L(cos t))|sS=s+1

L(f (t)) = e3s

Complete the square.

Replace s + 1 by S.
Basic Laplace table.

7.3 Laplace Transform Rules


= e3 e3s L(cos t) sS=s+1

= e3 (L(cos(t 3)H(t 3)))|sS=s+1

3

= e L(e

f (t) = e

3t

cos(t 3)H(t 3))

cos(t 3)H(t 3)

24 Example () Solve L(f (t) =

s2

263
Regroup factor e3S .
Second shifting rule.
First shifting rule.
Lerchs cancellation law.

s+7
for f (t).
+ 4s + 8

Solution: The answer is f (t) = e2t (cos 2t + 52 sin 2t). The details:
s+7
(s + 2)2 + 4
S+5
= 2
S +4
5 2
S
+
= 2
S + 4 2 S2 + 4

s
5 2 
= 2
+
s + 4 2 s2 + 4 

Complete the square.

L(f (t)) =

Replace s + 2 by S.
Split into table entries.
Prepare for shifting rule.
sS=s+2


= L(cos 2t) + 25 L(sin 2t)sS=s+2

= L(e

f (t) = e

2t

2t

(cos 2t +

(cos 2t +

5
2

5
2

sin 2t))

sin 2t)

Basic Laplace table.

First shifting rule.
Lerchs cancellation law.

264

Laplace Transform

7.4 Heavisides Method

This practical method was popularized by the English electrical engineer
Oliver Heaviside (18501925). A typical application of the method is to
solve
2s
= L(f (t))
(s + 1)(s2 + 1)
for the t-expression f (t) = et + cos t + sin t. The details in Heavisides
method involve a sequence of easy-to-learn college algebra steps.
More precisely, Heavisides method systematically converts a polynomial quotient
a0 + a1 s + + an sn
(1)
b0 + b1 s + + bm sm
into the form L(f (t)) for some expression f (t). It is assumed that
a0 , .., an , b0 , . . . , bm are constants and the polynomial quotient (1) has
limit zero at s = .

Partial Fraction Theory

In college algebra, it is shown that a rational function (1) can be expressed as the sum of terms of the form
(2)

A
(s s0 )k

where A is a real or complex constant and (s s0 )k divides the denominator in (1). In particular, s0 is a root of the denominator in (1).
Assume fraction (1) has real coefficients. If s0 in (2) is real, then A is
real. If s0 = + i in (2) is complex, then (s s0 )k also appears, where
s0 = i is the complex conjugate of s0 . The corresponding terms
in (2) turn out to be complex conjugates of one another, which can be
combined in terms of real numbers B and C as
(3)

B+Cs
A
A
+
=
.
k
k
(s s0 )
(s s0 )
((s )2 + 2 )k

Simple Roots. Assume that (1) has real coefficients and the denominator of the fraction (1) has distinct real roots s1 , . . . , sN and distinct
complex roots 1 + i1 , . . . , M + iM . The partial fraction expansion
of (1) is a sum given in terms of real constants Ap , Bq , Cq by
(4)

N
M
X
X
Ap
Bq + Cq (s q )
a0 + a1 s + + an sn
=
+
.
m
b0 + b1 s + + bm s
s sp q=1 (s q )2 + q2
p=1

7.4 Heavisides Method

265

Multiple Roots. Assume (1) has real coefficients and the denominator of the fraction (1) has possibly multiple roots. Let Np be the
multiplicity of real root sp and let Mq be the multiplicity of complex root
q + iq , 1 p N , 1 q M . The partial fraction expansion of (1)
is given in terms of real constants Ap,k , Bq,k , Cq,k by
(5)

N
X

p=1 1kNp

M
X
X Bq,k + Cq,k (s q )
Ap,k
+
.
k
(s sp )
((s q )2 + q2 )k
q=1 1kM
q

Heavisides Coverup Method

The method applies only to the case of distinct roots of the denominator
in (1). Extensions to multiple-root cases can be made; see page 266.
To illustrate Oliver Heavisides ideas, consider the problem details
(6)

2s + 1
s(s 1)(s + 1)

A
B
C
+
+
s
s1 s+1

= L(A) + L(Bet ) + L(Cet )

= L(A + Bet + Cet )
The first line (6) uses college algebra partial fractions. The second and
third lines use the Laplace integral table and properties of L.

Heavisides mysterious method. Oliver Heaviside proposed to

find in (6) the constant C =

1
2

by a coverup method:






2s + 1
s(s 1)

s+1 =0

The instructions are to coverup the matching factors (s + 1) on the left

, then evaluate on the left at the root s which
and right with box
makes the contents of the box zero. The other terms on the right are
replaced by zero.
To justify Heavisides coverup method, multiply (6) by the denominator
s + 1 of partial fraction C/(s + 1):
(2s + 1) (s + 1)
s(s 1) (s + 1)

A (s + 1)
s

B (s + 1)
s1

C (s + 1)
(s + 1)

Set (s + 1) = 0 in the display. Cancellations left and right plus annihilation of two terms on the right gives Heavisides prescription
2s + 1 
= C.
s(s 1) s+1=0


266

Laplace Transform

The factor (s + 1) in (6) is by no means special: the same procedure

applies to find A and B. The method works for denominators with
simple roots, that is, no repeated roots are allowed.

Extension to Multiple Roots. An extension of Heavisides method

is possible for the case of repeated roots. The basic idea is to factorout
the repeats. To illustrate, consider the partial fraction expansion details
R=

1
(s + 1)2 (s + 2)

A sample rational function having

repeated roots.

1
1
=
s + 1 (s + 1)(s + 2)


1
1
1
+
=
s+1 s+1 s+2


1
1
+
2
(s + 1)
(s + 1)(s + 2)
1
1
1
=
+
+
(s + 1)2
s+1 s+2
=

Factorout the repeats.

Apply the coverup method to the
simple root fraction.
Multiply.
Apply the coverup method to the
last fraction on the right.

Terms with only one root in the denominator are already partial fractions. Thus the work centers on expansion of quotients in which the
denominator has two or more roots.

Special Methods. Heavisides method has a useful extension for the

case of roots of multiplicity two. To illustrate, consider these details:
1
(s + 1)2 (s + 2)
C
A
B
+
=
+
2
s + 1 (s + 1)
s+2
1
A
1
+
=
+
2
s + 1 (s + 1)
s+2

R=

1
1
1
+
+
s + 1 (s + 1)2 s + 2

A fraction with multiple roots.

See equation (5).
Find B and C by Heavisides cover
up method.
Multiply by s+1. Set s = . Then
0 = A + 1.

The illustration works for one root of multiplicity two, because s =

will resolve the coefficient not found by the coverup method.
In general, if the denominator in (1) has a root s0 of multiplicity k, then
the partial fraction expansion contains terms
A2
Ak
A1
+
+ +
.
s s0 (s s0 )2
(s s0 )k
Heavisides coverup method directly finds Ak , but not A1 to Ak1 .

7.5 Heaviside Step and Dirac Delta

267

7.5 Heaviside Step and Dirac Delta

Heaviside Function. The unit step function or Heaviside function is defined by
H(x) =

1
0

for x 0,
for x < 0.

The most oftenused formula involving the Heaviside function is the

characteristic function of the interval a t < b, given by
(1)

H(t a) H(t b) =

1 a t < b,
0 t < a, t b.

To illustrate, a square wave sqw(t) = (1)floor(t) can be written in the

series form

(1)n (H(t n) H(t n 1)).

n=0

Dirac Delta. A precise mathematical definition of the Dirac delta,

denoted , is not possible to give here. Following its inventor P. Dirac,
the definition should be
(t) = dH(t).
The latter is nonsensical, because the unit step does not have a calculus derivative at t = 0. However, dH(t) could have the meaning of
a Riemann-Stieltjes integrator, which restrains dH(t) to have meaning
only under an integral sign. It is in this sense that the Dirac delta is
defined.
What do we mean by the differential equation
x00 + 16x = 5(t t0 )?
The equation x00 + 16x = f (t) represents a spring-mass system without
damping having Hookes constant 16, subject to external force f (t). In
a mechanical context, the Dirac delta term 5(t t0 ) is an idealization
of a hammer-hit at time t = t0 > 0 with impulse 5.
More precisely, the forcing term f (t) can be formally written as a RiemannStieltjes integrator 5dH(tt0 ) where H is Heavisides unit step function.
The Dirac delta or derivative of the Heaviside unit step, nonsensical
as it may appear, is realized in applications via the two-sided or central
difference quotient
H(t + h) H(t h)
dH(t).
2h

268

Laplace Transform

Therefore, the force f (t) in the idealization 5(t t0 ) is given for h > 0
very small by the approximation
f (t) 5

H(t t0 + h) H(t t0 h)
.
2h

The impulse2 of the approximated force over a large interval [a, b] is

computed from
Z

f (t)dt 5

h
h

H(t t0 + h) H(t t0 h)
dt = 5,
2h

due to the integrand being 1/(2h) on |t t0 | < h and otherwise 0.

Modeling Impulses. One argument for the Dirac delta idealization

is that an infinity of choices exist for modeling an impulse. There are in
addition to the central difference quotient two other popular difference
quotients, the forward quotient (H(t + h) H(t))/h and the backward
quotient (H(t) H(t h))/h (h > 0 assumed). In reality, h is unknown
in any application, and the impulsive force of a hammer hit is hardly
constant, as is supposed by this naive modeling.
The modeling logic often applied for the Dirac delta is that the external
force f (t) is used in the model in a limited manner, in which only the
momentum p = mv is important. More
precisely, only the change in
R
momentum or impulse is important, ab f (t)dt = p = mv(b) mv(a).
The precise force f (t) is replaced during the modeling by a simplistic
piecewise-defined force that has exactly the same impulse p. The replacement is justified by arguing that if only the impulse is important,
and not the actual details of the force, then both models should give
similar results.

Function or Operator? The work of physics Nobel prize winner P.

Dirac (19021984) proceeded for about 20 years before the mathematical
community developed a sound mathematical theory for his impulsive
force representations. A systematic theory was developed in 1936 by
the soviet mathematician S. Sobolev. The French mathematician L.
Schwartz further developed the theory in 1945. He observed that the
idealization is not a function but an operator or linear functional, in
particular, maps or associates to each function (t) its value at t = 0, in
short, () = (0). This fact was observed early on by Dirac and others,
during the replacement of simplistic forces by . In Laplace theory, there
is a natural encounter with the ideas, because L(f (t)) routinely appears
on the right of the equation after transformation. This term, in the case
2

Momentum is defined to be mass times velocity. RIf the force f is given by Newtons
b
d
(mv(t)) and v(t) is velocity, then a f (t)dt = mv(b) mv(a) is the
law as f (t) = dt
net momentum or impulse.

7.5 Heaviside Step and Dirac Delta

269

of an impulsive force f (t) = c(H(tt0 h)H(tt0 +h))/(2h), evaluates

for t0 > 0 and t0 h > 0 as follows:

c
(H(t t0 h) H(t t0 + h))est dt
2h
0
Z t0 +h
c st
=
e dt
t0 h 2h
!
sh esh
e
= cest0
2sh

L(f (t)) =

esh esh
The factor
is approximately 1 for h > 0 small, because of
2sh
LHospitals rule. The immediate conclusion is that we should replace
the impulsive force f by an equivalent one f such that
L(f (t)) = cest0 .
Well, there is no such function f !
The apparent mathematical flaw in this idea was resolved by the work
of L. Schwartz on distributions. In short, there is a solid foundation
for introducing f , but unfortunately the mathematics involved is not
elementary nor especially accessible to those readers whose background
is just calculus.
Practising engineers and scientists might be able to ignore the vast literature on distributions, citing the example of physicist P. Dirac, who
succeeded in applying impulsive force ideas without the distribution theory developed by S. Sobolev and L. Schwartz. This will not be the case
for those who wish to read current literature on partial differential equations, because the work on distributions has forever changed the required
background for that topic.

270

Laplace Transform

7.6 Laplace Table Derivations

Verified here are two Laplace tables, the minimal Laplace Table 7.2-4
and its extension Table 7.2-5. Largely, this section is for reading, as it is
designed to enrich lectures and to aid readers who study alone.
Derivation of Laplace integral formulas in Table 7.2-4, page 252.

Proof of L(tn ) = n!/s1+n:

The first step is to evaluate L(tn ) for n = 0.

R
L(1) = 0 (1)est dt
Laplace integral of f (t) = 1.
t=

= (1/s)est |t=0
= 1/s

Evaluate the integral.

Assumed s > 0 to evaluate limt est .

The value of L(tn ) for n = 1 can be obtained by s-differentiation of the relation

L(1) = 1/s, as follows.
R
d
d
st
dt
Laplace integral for f (t) = 1.
ds L(1) = ds 0 (1)e
Rb
R b dF
R d st
d
Used ds
= 0 ds (e ) dt
a F dt = a ds dt.
R
= 0 (t)est dt
Calculus rule (eu )0 = u0 eu .
= L(t)

Definition of L(t).

Then
d
L(t) = ds
L(1)

d
ds
(1/s)
2

= 1/s

Rewrite last display.

Use L(1) = 1/s.
Differentiate.

d
L(t) and hence L(t2 ) = 2/s3 .
This idea can be repeated to give L(t2 ) = ds
d
n
n1
The pattern is L(t ) = ds L(t
) which gives L(tn ) = n!/s1+n .

Proof of L(eat ) = 1/(s a):

The result follows from L(1) = 1/s,

R
L(eat ) = 0 eat est dt
R
= 0 e(sa)t dt
R
= 0 eSt dt
= 1/S

= 1/(s a)

as follows.
Direct Laplace transform.
Use eA eB = eA+B .
Substitute S = s a.

Apply L(1) = 1/s.

Back-substitute S = s a.

Proof of L(cos bt) = s/(ss + b2 ) and L(sin bt) = b/(ss + b2 ):

Use will be made of Eulers formula ei = cos + i sin , usually first introduced

in trigonometry. In this formula, is a real number (in radians) and i = 1

is the complex unit.

7.6 Laplace Table Derivations

271

eibt est = (cos bt)est + i(sin bt)est

Substitute = bt into Eulers

formula and multiply by est .

Integrate t = 0 to t = . Use
properties of integrals.

eibt est dt =

R
(cos bt)est dt
0 R

+ i 0 (sin bt)est dt

R
1
= 0 (cos bt)est dt
s ib
R
+ i 0 (sin bt)est dt

1
= L(cos bt) + iL(sin bt)
s ib

Evaluate the left side using

L(eat ) = 1/(s a), a = ib.
Direct Laplace transform definition.

s + ib
= L(cos bt) + iL(sin bt)
s2 + b 2

Use complex rule 1/z = z/|z|2 ,

z= A + iB, z = A iB, |z| =
A2 + B 2 .

s
= L(cos bt)
s2 + b 2
b
= L(sin bt)
2
s + b2

Extract the real part.

Extract the imaginary part.

Derivation of Laplace integral formulas in Table 7.2-5, page 252.

Proof of the Heaviside formula L(H(t a)) = eas /s.

L(H(t a)) =

=
=
=

R
0

R
a

H(t a)est dt Direct Laplace transform. Assume a 0.

(1)est dt

s(x+a)

dx
0 (1)e
R
as
sx
e
(1)e
dx
0
as

=e

(1/s)

Because H(t a) = 0 for 0 t < a.

Change variables t = x + a.

Constant eas moves outside integral.

Apply L(1) = 1/s.

Proof of the Dirac delta formula L((t a)) = eas .

The definition of the delta function is a formal one, in which every occurrence of
(t a)dt under an integrand is replaced by dH(t a). The differential symbol
dH(t a) is taken in the sense of the Riemann-Stieltjes integral. This integral
is defined in [?] for monotonic integrators (x) as the limit
Z

f (x)d(x) = lim

N
X

n=1

f (xn )((xn ) (xn1 ))

where x0 = a, xN = b and x0 < x1 < < xN forms a partition of [a, b] whose

mesh approaches zero as N .

The steps in computing the Laplace integral of the delta function appear below.
Admittedly, the proof requires advanced calculus skills and a certain level of
mathematical maturity. The reward is a fuller understanding of the Dirac
symbol (x).
R
L((t a)) = 0 est (t a)dt
Laplace integral, a > 0 assumed.
R st
= 0 e dH(t a)
Replace (t a)dt by dH(t a).
R M st
= limM 0 e dH(t a) Definition of improper integral.

272

Laplace Transform
= esa

Explained below.

To explain the last step, apply the definition of the Riemann-Stieltjes integral:
Z

est dH(t a) = lim

N
1
X
n=0

estn (H(tn a) H(tn1 a))

where 0 = t0 < t1 < < tN = M is a partition of [0, M ] whose mesh

max1nN (tn tn1 ) approaches zero as N . Given a partition, if tn1 <
a tn , then H(tn a)H(tn1 a) = 1, otherwise this factor is zero. Therefore,
the sum reduces to a single term estn . This term approaches esa as N ,
because tn must approach a.
as
Proof of L(floor(t/a)) = s(1 e eas ) :
The library function floor present in computer languages C and Fortran is
defined by floor(x) = greatest whole integer x, e.g., floor(5.2) = 5 and
floor(1.9) = 2. The computation of the Laplace integral of floor(t) requires
ideas from infinite series, as follows.
R
F (s) = 0 floor(t)est dt
Laplace integral definition.
P R n+1
st
On n t < n + 1, floor(t) = n.
= n=0 n (n)e dt
P n ns
= n=0 (e
enss )
Evaluate each integral.
s
1 es P
sn
=
Common factor removed.
n=0 ne
s
x(1 x) P
n1
Define x = es .
=
n=0 nx
s
x(1 x) d P n
=
x
Term-by-term differentiation.
s
dx n=0
x(1 x) d 1
Geometric series sum.
=
s
dx 1 x
x
=
Compute the derivative, simplify.
s(1 x)
es
=
Substitute x = es .
s(1 es )
To evaluate the Laplace integral of floor(t/a), a change of variables is made.
R
L(floor(t/a)) = 0 floor(t/a)est dt
Laplace integral definition.
R
asr
= a 0 floor(r)e
dr
Change variables t = ar.
= aF (as)
eas
=
s(1 eas )

Apply the formula for F (s).

Simplify.

Proof of L(sqw(t/a)) = 1s tanh(as/2):

The square wave defined by sqw(x) = (1)floor(x) is periodic of period 2 and
R2
piecewise-defined. Let P = 0 sqw(t)est dt.

7.6 Laplace Table Derivations

P=
=

R1

273

R2
sqw(t)est dt + 1 sqw(t)est dt
R2
1 est dt

Apply

0
R 1 st
e dt
0

Rb
a

Rc
a

Rb
c

Use sqw(x) = 1 on 0 x < 1 and

sqw(x) = 1 on 1 x < 2.

1
1
(1 es ) + (e2s es )
s
s
1
= (1 es )2
s
=

Evaluate each integral.

Collect terms.

An intermediate step is to compute the Laplace integral of sqw(t):

R2
sqw(t)est dt
L(sqw(t)) = 0
Periodic function formula, page 275.
1 e2s
1
1
= (1 es )2
. Use the computation of P above.
s
1 e2s
1 1 es
=
.
Factor 1 e2s = (1 es )(1 + es ).
s 1 + es
1 es/2 es/2
.
Multiply the fraction by es/2 /es/2 .
=
s es/2 + es/2
1 sinh(s/2)
=
.
Use sinh u = (eu eu )/2,
s cosh(s/2)
cosh u = (eu + eu )/2.
1
= tanh(s/2).
Use tanh u = sinh u/ cosh u.
s
To complete the computation of L(sqw(t/a)), a change of variables is made:
R
L(sqw(t/a)) = 0 sqw(t/a)est dt
Direct transform.
R
asr
= 0 sqw(r)e
(a)dr
Change variables r = t/a.
a
=
tanh(as/2)
See L(sqw(t)) above.
as
1
= tanh(as/2)
s

Proof of L(a trw(t/a)) = s12 tanh(as/2):

The triangular wave is defined by trw(t) =
L(a trw(t/a)) =

1
(f (0) + L(f 0 (t))
s

1
L(sqw(t/a))
s

1
tanh(as/2)
s2

Rt
0

sqw(x)dx.

Let f (t) = a trw(t/a). Use L(f 0 (t)) =

sL(f (t)) f (0), page 251.
R t/a
Use f (0) = 0, (a 0 sqw(x)dx)0 =
sqw(t/a).
Table entry for sqw.

+ )
Proof of L(t ) = (1s1+
:
L(t ) =

R0
0

t est dt

Direct Laplace transform.

(u/s) eu du/s

Change variables u = st, du = sdt.

274

Laplace Transform
1 R u
u e du
s1+ 0
1
= 1+ (1 + ).
s
=

Where (x) =
definition.

R
0

ux1 eu du, by

The generalized factorial function (x) is defined for x > 0 and it agrees with
the classical factorial n! = (1)(2) (n) in case x = n + 1 is an integer. In
literature, ! means (1 + ). For more details about the Gamma function, see
Abramowitz and Stegun [?], or maple documentation.
r
1/2
Proof of L(t ) = s :
(1 + (1/2))
s11/2

=
s

L(t1/2 ) =

Apply the previous formula.

Use (1/2) =

7.7 Transform Properties

275

7.7 Transform Properties

Collected here are the major theorems and their proofs for the manipulation of Laplace transform tables.
Theorem 4 (Linearity)
The Laplace transform has these inherited integral properties:
(a)
(b)

L(f (t) + g(t)) = L(f (t)) + L(g(t)),

L(cf (t)) = cL(f (t)).

Theorem 5 (The t-Derivative Rule)

Let y(t) be continuous, of exponential order and let f 0 (t) be piecewise
continuous on t 0. Then L(y 0 (t)) exists and
L(y 0 (t)) = sL(y(t)) y(0).
Theorem 6 (The t-Integral Rule)
Let g(t) be of exponential order and continuous for t 0. Then
L

R
t
0


1
g(x) dx = L(g(t)).
s

Theorem 7 (The s-Differentiation Rule)

Let f (t) be of exponential order. Then
d
L(f (t)).
ds

L(tf (t)) =

Theorem 8 (First Shifting Rule)

Let f (t) be of exponential order and < a < . Then
L(eat f (t)) = L(f (t))|s(sa) .
Theorem 9 (Second Shifting Rule)
Let f (t) and g(t) be of exponential order and assume a 0. Then
(a)
(b)

L(f (t a)H(t a)) = eas L(f (t)),

L(g(t)H(t a)) = eas L(g(t + a)).

Theorem 10 (Periodic Function Rule)

Let f (t) be of exponential order and satisfy f (t + P ) = f (t). Then
L(f (t)) =

RP
0

f (t)est dt
.
1 eP s

Theorem 11 (Convolution Rule)

Let f (t) and g(t) be of exponential order. Then
L(f (t))L(g(t)) = L

Z

f (x)g(t x)dx .

276

Laplace Transform

Proof of Theorem 4 (linearity):

LHS = L(f (t) + g(t))
R
= 0 (f (t) + g(t))est dt
R
R
= 0 f (t)est dt + 0 g(t)est dt
= L(f (t)) + L(g(t))

LHS = L(cf (t))

R
= 0 cf (t)est dt
R
= c 0 f (t)est dt

Left side of the identity in (a).

Direct transform.
Calculus integral rule.
Equals RHS; identity (a) verified.
Left side of the identity in (b).
Direct transform.
Calculus integral rule.

= cL(f (t))

Equals RHS; identity (b) verified.

Proof of Theorem 5 (t-derivative rule): Already L(f (t)) exists, because

f is of exponential order and continuous. On an interval [a, b] where f 0 is
continuous, integration by parts using u = est , dv = f 0 (t)dt gives
Rb
a

t=b

f 0 (t)est dt = f (t)est |t=a

Rb

f (t)(s)est dt
Rb
= f (a)esa + f (b)esb + s a f (t)est dt.
a

On any interval [0, N ], there are finitely many intervals [a, b] on each of which
f 0 is continuous. Add the above equality across these finitely many intervals
[a, b]. The boundary values on adjacent intervals match and the integrals add
to give
Z

f 0 (t)est dt = f (0)e0 + f (N )esN + s

f (t)est dt.

Take the limit across this equality as N . Then the right side has limit
f (0) + sL(f (t)), because of the existence of L(f (t)) and limt f (t)est = 0
for large s. Therefore, the left side has a limit, and by definition L(f 0 (t)) exists
and L(f 0 (t)) = f (0) + sL(f (t)).

Proof of Theorem 6 (t-Integral rule): Let f (t) =

Rt

0 g(x)dx. Then f is of
exponential order and continuous. The details:
Rt
L( 0 g(x)dx) = L(f (t))
By definition.
1
= L(f 0 (t)) Because f (0) = 0 implies L(f 0 (t)) = sL(f (t)).
s
1
= L(g(t))
Because f 0 = g by the Fundamental theorem of
s
calculus.

Proof of Theorem 7 (s-differentiation): We prove the equivalent relation

L((t)f (t)) = (d/ds)L(f (t)). If f is of exponential order, then so is (t)f (t),
therefore L((t)f (t)) exists. It remains to show the s-derivative exists and
satisfies the given equality.
The proof below is based in part upon the calculus inequality
 x

e + x 1 x2 , x 0.
(1)

7.7 Transform Properties

277

The inequality is obtained from two applications of the mean value theorem
g(b)g(a) = g 0 (x)(ba), which gives ex +x1 = xxex1 with 0 x1 x x.

In addition, the existence of L(t2 |f (t)|) is used to define s0 > 0 such that
L(t2 |f (t)|) 1 for s > s0 . This follows from the transform existence theorem
for functions of exponential order, where it is shown that the transform has
limit zero at s = .
Consider h 6= 0 and the Newton quotient Q(s, h) = (F (s + h) F (s))/h for the
s-derivative of the Laplace integral. We have to show that
lim |Q(s, h) L((t)f (t))| = 0.

h0

This will be accomplished by proving for s > s0 and s + h > s0 the inequality
|Q(s, h) L((t)f (t))| |h|.
For h 6= 0,
Q(s, h) L((t)f (t)) =

f (t)

estht est + thest

dt.
h

Assume h > 0. Due to the exponential rule eA+B = eA eB , the quotient in the
integrand simplifies to give
 ht

Z
e
+ th 1
st
Q(s, h) L((t)f (t)) =
f (t)e
dt.
h
0
Inequality (1) applies with x = ht 0, giving
Z
|Q(s, h) L((t)f (t))| |h|

t2 |f (t)|est dt.

The right side is |h|L(t2 |f (t)|), which for s > s0 is bounded by |h|, completing
the proof for h > 0. If h < 0, then a similar calculation is made to obtain
Z
|Q(s, h) L((t)f (t))| |h|
t2 |f (t)estht dt.
0

The right side is |h|L(t2 |f (t)|) evaluated at s + h instead of s. If s + h > s0 ,

then the right side is bounded by |h|, completing the proof for h < 0.

Proof of Theorem 8 (first shifting rule): The left side LHS of the equality
can be written because of the exponential rule eA eB = eA+B as
Z
LHS =
f (t)e(sa)t dt.
0

This integral is L(f (t)) with s replaced by s a, which is precisely the meaning
of the right side RHS of the equality. Therefore, LHS = RHS.

Proof of Theorem 9 (second shifting rule): The details for (a) are
LHS = L(H(t a)f (t a))
R
= 0 H(t a)f (t a)est dt Direct transform.

278

Laplace Transform
=
=
=

H(t a)f (t a)est dt

a
R
s(x+a)
dx
0 H(x)f (x)e
R
sa
sx
e
dx
0 f (x)e
sa

=e

= RHS

L(f (t))

Because a 0 and H(x) = 0 for x < 0.

Change variables x = t a, dx = dt.
Use H(x) = 1 for x 0.

Direct transform.

Identity (a) verified.

In the details for (b), let f (t) = g(t + a), then

LHS = L(H(t a)g(t))

= L(H(t a)f (t a))

Use f (t a) = g(t a + a) = g(t).

= esa L(g(t + a))

Because f (t) = g(t + a).

= esa L(f (t))

Apply (a).

= RHS

Identity (b) verified.

Proof of Theorem 10 (periodic function rule):

LHS = L(f (t))
R
= 0 f (t)est dt
P R nP +P
= n=0 nP
f (t)est dt
P R P
= n=0 0 f (x + nP )esxnP s dx
RP
P
= n=0 enP s 0 f (x)esx dx
=

RP
0

RP
0

RP

f (x)esx dx

f (x)esx dx

f (x)esx dx
1 eP s
= RHS

n=0

1
1r

rn

Direct transform.
Additivity of the integral.
Change variables t = x + nP .
Because f is P -periodic and
eA eB = eA+B .
Common factor in summation.
Define r = eP s .
Sum the geometric series.
Substitute r = eP s .
Periodic function identity verified.

Left unmentioned here is the convergence of the infinite series on line 3 of the
proof, which follows from f of exponential order.

Proof of Theorem 11 (convolution rule): The details use Fubinis integration interchange theorem for a planar unbounded region, and therefore
this proof involves advanced calculus methods that may be outside the background of the reader. Modern calculus texts contain a less general version of
Fubinis theorem for finite regions, usually referenced as iterated integrals. The
unbounded planar region is written in two ways:
D = {(r, t) : t r < , 0 t < },

D = {(r, t) : 0 r < , 0 r t}.

Readers should pause here and verify that D = D.

7.7 Transform Properties

279

The change of variable r = x + t, dr = dx is applied for fixed t 0 to obtain

the identity
R
R
est 0 g(x)esx dx = 0 g(x)esxst dx
(2)
R
= t g(r t)ers dr.
The left side of the convolution identity is expanded as follows:
LHS = L(f (t))L(g(t))
R
R
= 0 f (t)est dt 0 g(x)esx dx
R
R
= 0 f (t) t g(r t)ers drdt
R
= D f (t)g(r t)ers drdt
R
= D f (t)g(r t)ers drdt
RRr
= 0 0 f (t)g(r t)dters dr

Direct transform.
Apply identity (2).
Fubinis theorem applied.
Descriptions D and D are the same.
Fubinis theorem applied.

Then

R

t
RHS = L 0 f (u)g(t u)du
RRt
= 0 0 f (u)g(t u)duest dt
RRr
= 0 0 f (u)g(r u)duesr dr
RRr
= 0 0 f (t)g(r t)dt esr dr
= LHS

Direct transform.
Change variable names r t.

Change variable names u t.

Convolution identity verified.

280

Laplace Transform

7.8 More on the Laplace Transform

Model conversion and engineering. A differential equation model
for a physical system can be subjected to the Laplace transform in order
to produce an algebraic model in the transform variable s. Lerchs theorem says that both models are equivalent, that is, the solution of one
model gives the solution to the other model.
In electrical and computer engineering it is commonplace to deal only
with the Laplace algebraic model. Engineers are in fact capable of having hour-long modeling conversations, during which differential equations
are never referenced! Terminology for such modeling is necessarily specialized, which gives rise to new contextual meanings to the terms input
and output. For example, an RLC-circuit would be discussed with input
F (s) =

s2

,
+ 2

and the listener must know that this expression is the Laplace transform
of the t-expression sin t. Hence the RLC-circuit is driven by a sinusoindal input of natural frequency . During the modeling discourse, it
could be that the output is
X(s) =

10
1
.
+ 2
s + 1 s + 2

Lerchs equivalence says that X(s) is the Laplace transform of et +

10 sin t, but that is extra work, if all that is needed from the model is a
statement about the transient and steady-state responses to the input.