# Laplace Transformation

This tutorial contains many *Mathematica* scripts. You, as the user, are free to use all codes for your needs, and have the
right to distribute this tutorial and refer to this tutorial as long as
this tutorial is accredited appropriately. Any comments and/or contributions for this tutorial are welcome; you can send your remarks to <Vladimir_Dobrushkin@brown.edu>

Return to computing page for the second course APMA0330

Return to Mathematica tutorial for the first course

Return to Mathematica tutorial for the second course

Return to the main page for the course APMA0330

Return to the main page for the course APMA0340

# Brief History

The Laplace Transform is named after the French mathematician and astronomer Pierre-Simon Laplace (1749--1827). However, he did not actually invent what we now call the Laplace transform. Indeed, Laplace himself, a notoriously vain and selfish
person in spite of his scientific genius, was careful to credit Leonhard Euler (1707--1783) with the basic formula. Well before the
work of Laplace, however, mathematical genius Leonhard Euler had studied differential equations. One of his many
noteworthy contributions in this field was the idea of transforming a function *X(x)* into a new function *z* via
the equation

*Methodus aequationes differentiales altiorum graduum integrandi ulterius promota*-- it’s a good thing mathematicians don’t use Latin any more…), Euler used methods based on this transform to give a systematic method of solving second order linear differential equations. Later in his career, he further clarified the method and introduced the definite integral form

*Institutiones Calculi Integralis*, vol. II, where he used it to solve the equation

However, Euler did not pursue this topic very far. Joseph Louis Lagrange (1736--1813), born as Giuseppe Lodovico Lagrangia in Turin, Italy, who succeded Euler (since Leonhard returned to Russia) as the director of mathematics at the Prussian Academy of Sciences in Berlin, began to study integrals in the form \( \int_0^{\infty} f(t)\,e^{-at}\,\mathrm{d}t \) in connection with his work on integrating probability density functions. Laplace was the next person to seriously work on this topic, and took a critical step forward by applying the idea of a "transformation" rather than just looking for a solution in the form of an integral. He looked for solutions with the following equation: \( \int x^{s}\,\phi (x)\,\mathrm{d}x = F(s) .\) In 1809, Laplace extended his transform to find solutions that diffused indefinitely -- giving us our popular Laplace transformation. Laplace appeared to have quickly understood the importance of his discovery, as he went on to use Laplace transforms numerous times in his later work and generalize his integrals to create Fourier and Mellin transforms as well. So while Laplace may not have invented his transforms, he certainly deserves credit for producing a systematic body of theory that went far beyond anything created by his predecessors.

Although the results had been published for at least 70 years, the transformation was not given a true physical and mathematical meaning until Oliver Heaviside (1850--1925) came up with completely new ideas on his own in the 1880s. His predecessors used (what we call now the Laplace) integral as an anlytical tool and never even tried to establish the inverse Laplace transform---without this tool the theory is not complete at all. Oliver Heaviside did not actually use (and most likely was unfamiliar with) the Laplace integral in his derivations because it was not needed in his pioneering work in establishing the operator methods.

Let us pause explanations of Heaviside's achievements and focus on his biography and historical circumstances that surrounded his science breakthrough. Oliver Heaviside was born on May 18, 1850 in Camden Town, London, England. He caught scarlet fever when he was a young child and this affected his hearing. This was to have a major effect on his life, making his childhood unhappy, with relations between himself and other children difficult. However, his school results were rather good, and in 1865 he was placed fifth from five hundred pupils. Academic subjects seemed to hold little attraction for Heaviside, however, and at age 16 he left school. Perhaps he was more disillusioned with school than with learning since he continued to study (Oliver was completely self-taught) after leaving school, in particular he learnt Morse code, studied electricity, and studied further languages, in particular Danish and German. He was aiming at a career as a telegrapher, in this he was advised and helped by his uncle Charles Wheatstone (the piece of electrical apparatus known as the Wheatstone bridge is named after him).

In 1868 Heaviside went to Denmark and became a telegrapher. He
progressed quickly in his profession and returned to England in 1871
to take up a post in Newcastle upon Tyne in the office of Great
Northern Telegraph Company, which dealt with overseas
traffic. Heaviside became increasingly deaf but he worked on his own
research into electricity. While still working as chief operator
in Newcastle he began to publish papers on electricity, the first in
1872, and then the second in 1873, which was of sufficient interest to
James Clerk Maxwell (1831--1879) that he mentioned the results in the
second edition of his
*Treatise on Electricity and Magnetism*. Despite this hatred of
rigour, Heaviside was able to greatly simplify Maxwell's twenty
equations with twenty variables, replacing them by four equations with two
vector variables (the electric field **E** and the magnetic
field **B**). Today we call these 'Maxwell's equations' forgetting that
they are in fact 'Heaviside's equations.'

Heaviside went on to achieve further advances in knowledge, again
receiving less than his just desserts. In a 1887 paper Heaviside gave,
for the first time, the conditions necessary to transmit a signal
without distortion. In *Electromagnetic Theory* (1893--1912), he
postulated that an electric charge would increase in mass as its
velocity increases, an anticipation of an aspect of Einstein’s special
theory of relativity. Heaviside was elected a Fellow of the Royal
Society in 1891, perhaps the greatest honour he ever received. In 1902
Heaviside predicted that there was a conducting layer in the
atmosphere which allowed radio waves to follow the Earth's
curvature. This layer in the atmosphere, the Heaviside layer, is named
after him. Its existence was proved in 1923 when radio pulses were
transmitted vertically upward and the returning pulses from the
reflecting layer were received.

On the eve of the twentieth century, new physical ideas required new mathematical tools to utilize these theories. It is not an accident that many new mathematical notations and definitions were introduced by physicists, but not mathematicians. We cannot not to mention the famous Dirac's delta-function as well as Einstein summation abbreviation. The latest example gave us the Internet, including the HTML protocol. Therefore, Heaviside's invention of operational calculus has naturally evolved to quantum mechanics and then as a byproduct, a useful technique for solving initial value problems for differential equations.

The revolutionary idea of Heaviside's operational calculus consisted in establishing the spectral decomposition of the (unbounded) derivative operator \( \texttt{D} = {\text d}/{\text d} t \) acting in the space of functions on half line \( [0, \infty ) . \) This means that the operator was replaced by simple multiplication. We will meet this approach in the second part of this tutorial when spectral decomposition will be applied to square matrices.

Due to his research, Heaviside dealt with multi-degree differential
equations for electrical systems in 1880s. To solve corresponding
initial value problems, Heaviside substituted the derivative operator by a
letter *p* (we use contemporary notation \(
\texttt{D} \) instead), which yields an algebraic
equations. After that he quickly solved the algebraic equation and proceeded with a solution to the
original differential equations. These
kinds of equations (usually with 10 or more derivatives of dependent
variable) would usually take days or week for most people to solve.
Heaviside was able to solve these kinds of equations within hours.
Single-handedly, between 1880--1887, he invented operational calculus
-- a new method for solving the differential equations.

Unfortunately, Oliver Heaviside failed to explain how he had derived the solutions to the initial value problems, and as often happens the work of geniuses is not often understood. It took many years for academia to accept his results because of a lack of a rigorous proof for his methods. He replied to this criticism with the famous statement "Mathematics is an experimental science, and definitions do not come first, but later on," claiming that "I do not refuse my dinner simply because I do not understand the process of digestion." As a result, mathematicians did not accept his technique and tried to explain what is going on behind his manipulations. After many years of work, it became clear that Heaviside used the integral transformation

*s*in the complex plane such that

*s*is greater than the real part of all singularities of

*F(λ)*and

*F(λ)*is bounded on the line of integration. This integral is usually referred to as the Bromwich integral. To this day, the Inverse Laplace Transform is the most difficult process to understand when solving differential equations in this method. In this tutorial, we will not use the above formal definition of the inverse Laplace transform; instead, we apply the residue method based on the novel approach by Vladimir Dobrushkin.

# Concept Overview

Essentially, we will be doing almost the same thing that Heaviside had done-- we will take a differential equation, subject to the initial conditions

*n*and

*P(λ)*is the polynomial of degree less than

*n*in λ depending on the initial conditions; finally, we apply the inverse Laplace transform to obtain the required solution. We again will consider Initial Value Problems (IVP) with specified initial conditions that have unique solutions. Thus, our answer will also be unique when using the Laplace Transform method -- we will have no arbitrary constants. Something to note, however, is that this method is ONLY valid for "

*t*" (our independent variable) when "

*t*" is > 0. We will be using a Heaviside function that will allow our answer to maintain these conditions (the Heaviside function will be more well defined in the section labled "Heaviside Function")