Section 3.4 : Repeated Roots
In this section we will be looking at the last case for the constant coefficient, linear, homogeneous second order differential equations. In this case we want solutions to
ay″+by′+cy=0where solutions to the characteristic equation
ar2+br+c=0are double roots r1=r2=r.
This leads to a problem however. Recall that the solutions are
y1(t)=er1t=erty2(t)=er2t=ertThese are the same solution and will NOT be “nice enough” to form a general solution. We do promise that we’ll define “nice enough” eventually! So, we can use the first solution, but we’re going to need a second solution.
Before finding this second solution let’s take a little side trip. The reason for the side trip will be clear eventually. From the quadratic formula we know that the roots to the characteristic equation are,
r1,2=−b±√b2−4ac2aIn this case, since we have double roots we must have
b2−4ac=0This is the only way that we can get double roots and in this case the roots will be
r1,2=−b2aSo, the one solution that we’ve got is
y1(t)=e−bt2aTo find a second solution we will use the fact that a constant times a solution to a linear homogeneous differential equation is also a solution. If this is true then maybe we’ll get lucky and the following will also be a solution
y2(t)=v(t)y1(t)=v(t)e−bt2awith a proper choice of v(t). To determine if this in fact can be done, let’s plug this back into the differential equation and see what we get. We’ll first need a couple of derivatives.
y′2(t)=v′e−bt2a−b2ave−bt2ay″2(t)=v″e−bt2a−b2av′e−bt2a−b2av′e−bt2a+b24a2ve−bt2a=v″e−bt2a−bav′e−bt2a+b24a2ve−bt2aWe dropped the (t) part on the v to simplify things a little for the writing out of the derivatives. Now, plug these into the differential equation.
a(v″e−bt2a−bav′e−bt2a+b24a2ve−bt2a)+b(v′e−bt2a−b2ave−bt2a)+c(ve−bt2a)=0We can factor an exponential out of all the terms so let’s do that. We’ll also collect all the coefficients of v and its derivatives.
e−bt2a(av″+(−b+b)v′+(b24a−b22a+c)v)=0e−bt2a(av″+(−b24a+c)v)=0e−bt2a(av″−14a(b2−4ac)v)=0Now, because we are working with a double root we know that that the second term will be zero. Also exponentials are never zero. Therefore, (1) will be a solution to the differential equation provided v(t) is a function that satisfies the following differential equation.
av″=0ORv″=0We can drop the a because we know that it can’t be zero. If it were we wouldn’t have a second order differential equation! So, we can now determine the most general possible form that is allowable for v(t).
v′=∫v″dt=cv(t)=∫v′dt=ct+kThis is actually more complicated than we need and in fact we can drop both of the constants from this. To see why this is let’s go ahead and use this to get the second solution. The two solutions are then
y1(t)=e−bt2ay2(t)=(ct+k)e−bt2aEventually you will be able to show that these two solutions are “nice enough” to form a general solution. The general solution would then be the following.
y(t)=c1e−bt2a+c2(ct+k)e−bt2a=c1e−bt2a+(c2ct+c2k)e−bt2a=(c1+c2k)e−bt2a+c2cte−bt2aNotice that we rearranged things a little. Now, c, k, c1, and c2 are all unknown constants so any combination of them will also be unknown constants. In particular, c1+c2k and c2c are unknown constants so we’ll just rewrite them as follows.
y(t)=c1e−bt2a+c2te−bt2aSo, if we go back to the most general form for v(t) we can take c=1 and k=0 and we will arrive at the same general solution.
Let’s recap. If the roots of the characteristic equation are r1=r2=r, then the general solution is then
y(t)=c1ert+c2tertNow, let’s work a couple of examples.
The characteristic equation and its roots are.
r2−4r+4=(r−2)2=0r1,2=2The general solution and its derivative are
y(t)=c1e2t+c2te2ty′(t)=2c1e2t+c2e2t+2c2te2tDon’t forget to product rule the second term! Plugging in the initial conditions gives the following system.
12=y(0)=c1−3=y′(0)=2c1+c2This system is easily solved to get c1=12 and c2=−27. The actual solution to the IVP is then.
y(t)=12e2t−27te2tThe characteristic equation and its roots are.
16r2−40r+25=(4r−5)2=0r1,2=54The general solution and its derivative are
y(t)=c1e5t4+c2te5t4y′(t)=54c1e5t4+c2e5t4+54c2te5t4Don’t forget to product rule the second term! Plugging in the initial conditions gives the following system.
3=y(0)=c1−94=y′(0)=54c1+c2This system is easily solve to get c1=3 and c2=−6. The actual solution to the IVP is then.
y(t)=3e5t4−6te5t4The characteristic equation and its roots are.
r2+14r+49=(r+7)2=0r1,2=−7The general solution and its derivative are
y(t)=c1e−7t+c2te−7ty′(t)=−7c1e−7t+c2e−7t−7c2te−7tPlugging in the initial conditions gives the following system of equations.
−1=y(−4)=c1e28−4c2e285=y′(−4)=−7c1e28+c2e28+28c2e28=−7c1e28+29c2e28Solving this system gives the following constants.
c1=−9e−28c2=−2e−28The actual solution to the IVP is then.
y(t)=−9e−28e−7t−2te−28e−7ty(t)=−9e−7(t+4)−2te−7(t+4)