Solution of the Ordinary Differential Equations

   

The Differential Equation in t

The ordinary differential equation in time is an old friend.

Clear[w]

timeDiffEqn2

Its solution is

Clear[wSolution,t]

wSolution[t_] := c1 Sin[c k t] + c2 Cos[c k t]

We can satisfy the initial condition, du/dt (r,0) = 0 by setting

initialVelocity = wSolution'[0] == 0

We know that neither the k's nor c are zero, so we must set c1 = 0. So now

Clear[wSolution,t]

wSolution[t_] := c2 Cos[c k t]

Substituting the values of k we found in the previous section gives us a list of solutions to the wave equation, w[t] v[r].

solutionList = Thread[Times[Array[c,Length[j0zeroes]],

wSolution[t]/c2 BesselJ[0,k r] /. j0zeroes]]

Each of these solutions satisfies the boundary conditions u(0,t) bounded and u(1,t) = 0, and the initial condition du/dt (r,0) = 0.

Go up to Vibrating Drumhead