I know the large majority of people on here might not know what difference equations are (or maybe you know them as recurrence relations, because that's what Wikipedia lists them as), but I got a question that maybe someone can help me out with.

I have a difference equation of the form:

y(k+2) + y(k+1) + y(k) = b * r^k

My particular equation is this:

y(k+2) + y(k+1) - 2y(k) = 2 * 1^k

Now, this is a non-homogeneous difference equation, so the first thing I did is create an associated homogeneous difference equation:

z(k+2) + z(k+1) - 2z(k) = 0

Now, if weassume that z(k) = Lambda^k, we get:

(I am going to use the letter L for lambda)

L^(k+2) + L^(k+1) - 2 * L^k = 0

Simplify:

L^2 + L - 2 = 0

We can then solve for lambda, and we get

L1 = 1

L2 = -2

Great. We have our lambda values. Now we need to find the general solution. We know that

y(k) = y*(k) + z(k)

where y*(k) is a particular solution to y(k) and z(k) is a general solution to our homogenous equation. Now, in a normal situation, I would do the following:

y*(k) = d * r^k

And then I would solve for "d".

My problem is that "r" is equal to one of our roots (lambda). In other words, r = Lambda_1 (because r is 1 in our very first equation up above).

So if I do y*(k) = dr^k, we end up solving it, and it turns out that 0 = 2, which we know is wrong.

How can I solve for y*(k)?

I was reading some notes, and I read that"when r is equal to a root, treat it as a repeated root of the extended polynomial."The notes gave some information on how to handle that in z(k)...but not in y*(k). Basically, right now I know that:

y(k) = y*(k) + c1 * L1^k + c2 * k * L2^2

But I need to find y*(k) and I cannot figure out how to do that. Anyone know how to do that?