Homogeneous Linear ODEs with Constant Coefficients
As I was finalizing my study of the method of separation of variables, I ran into the inconsistency that I had forgotten how to solve homogeneous linear ODEs with constant coefficients, that is, equations of the form Basically, we perform the following lambda substitution: and solve the resulting quadratic equation for : There are three cases:
- If
, then we have a solution of the form - If
, then we have a solution of the form - If
, then we have a solution of the form
The last one is derived from Euler's formula , and the solution is taken to have the form .
Example To find the complementary function , we perform a lambda substitution and solve: The particular solution has the form . Differentiating and substituting accordingly gives Therefore, by the superposition principle, the general solution to the ODE is
The Method of Separation of Variables
We want to solve the PDE with initial and boundary conditions To do that, let's assume the solution has the form Don't ask, just accept it. There's a reason behind picking such form, but it wouldn't help to explain it now. With that said, let's plug this solution into our PDE to obtain Now, let's separate variables! If you look at what we just obtained, you're eventually going to ask yourself: "how's it possible for a function in terms of to be equal to a function in terms of ? This makes no sense at all!" Well, of course it doesn't, unless they're both equal to the same constant! Let's call this constant : Again, the reason we separated our equation this way and the reason we put a negative sign before are for the purpose of convenience only. It'll all make sense down the road, so stop scratching your head and keep reading! By the way, notice that our partial derivatives turned into normal ones because we're dealing with functions with one variable.
Now, let's solve the first ODE Our famous toolkit methods lead to the solution With that out of the way, we can now solve the second ODE (it's a BVP) Performing a lambda substitution would be confusing, so let's call it an substitution. It'll yield Again we have three cases:
Now, let's solve the first ODE
- If
, then From the boundary conditions, To avoid the trivial solution, we must have that where . As a result, we're left with - If
, then From the boundary conditions, Therefore, is not an eigenvalue for this problem.
- If
, then where . Let's use the hyperbolic functions to represent this From the boundary conditions, because and is never zero for a positive argument. Therefore, is not an eigenvalue for this problem, either.
In summary, we obtained product solutions of the heat equation Finally, by the principle of superposition, we claim that the following infinite series is the solution of our heat conduction problem: Phew! I guess I could now do the same thing for the multidimensional case, lol.