Colin Mitchell
Now that we have formed our model,

with
, we can solve it using the method of separation of
variables. So we want to seek a solution
of the form
.
So inserting
![]()
into the governing equation gives us

Since the left hand side is dependent only on x, and the right hand side is dependent only on t, then both sides must be equal to a constant. So we can form ordinary differential equations in both functions,

Let
, so that the equations reduce,
.
Since
can be any number
, we notice that different cases will give different
solutions to each ODE. Basic intuition
allows us to understand that our model will have a sinusoidal shape, and as
such we should seek a solution of sines and cosines. For
, we will get solutions of the form

which will satisfy the governing
PDE, but will not give a proper model. For
, we will get solutions of the form
![]()
The boundary conditions will give a trivial solution, and
this form will also not work. So we must
assume that
. First, let us solve
for X. The solution for
is of the form
![]()
Using the left boundary condition
, we get
![]()
Using the right boundary condition
, we get
![]()
Since we do not want the trivial solution
, we must have
. This then tells us that
, so
. So we reach
.
This is a sequence of solutions, each of which satisfies the
PDE and the boundary conditions. We
notice that we will have constants in both X
and T, and so we will omit the constant
of amplification in X, and T will define the amplitude of the
vibrations. So we can write summation of
each
as the overall form of
X,
![]()
Now we will
focus on T. Similarly, the solution to
will take the form
. Since we know that
, we will replace it in this equation to get

It is also apparent here that T is also a sequence of solutions, and so we can combine it with X to reach the solution
.
Our
final step is to fit the solution to the initial conditions
and
. Using the first
condition, we get

Now, we want to find the coefficients
. If we take
, multiply by
, and integrate from 0 to L,
we get
.
Using the integral
,
we can reduce this to
.
Now we must
consider the initial condition
. Inserting the
solution into this, we get

Again, we can multiply by
and integrate to get
,
and reducing with the integral identity, we get
