Solve Difference Equation for c_n

Join the discussion
Ask a follow-up here, or get your own question answered by working scientists, mathematicians and engineers — people, not an autocomplete.
Real named experts · corrections over time · the nuance an AI answer skips
5 replies · 3K views
foxjwill
Messages
350
Reaction score
0

Homework Statement


Find a closed-form expression for [tex]c_n[/tex].

[tex]c_{n+1} = \frac{c_0 (3c_n + c_0)}{2c_n + c_0}[/tex]


Homework Equations





The Attempt at a Solution


Besides finding [tex]c_1, c_2, c_3, \ldots[/tex] and looking for a pattern, I have absolutely no idea.
 
Physics news on Phys.org
Every recurrence equation of the form
[tex]c_{n+1}=\frac{\kappa\,c_n+\lambda}{\mu\,c_n+\nu},\quad \kappa\,\nu-\mu\,\lambda\neq0[/tex]
can be brought into the linear form [itex]x_{n+1}=x_n+b_n[/itex] with [itex]b_n[/itex] known.
The transformation that does that, is
[tex]c_n=\frac{\alpha^{-n}}{x_n}+\beta[/tex]
Plug this to your equation and choose [itex]\alpha,\beta[/itex] in order to arrive to [itex]x_{n+1}=x_n+b_n[/itex]
 
Rainbow Child said:
Every recurrence equation of the form
[tex]c_{n+1}=\frac{\kappa\,c_n+\lambda}{\mu\,c_n+\nu},\quad \kappa\,\nu-\mu\,\lambda\neq0[/tex]
can be brought into the linear form [itex]x_{n+1}=x_n+b_n[/itex] with [itex]b_n[/itex] known.
The transformation that does that, is
[tex]c_n=\frac{\alpha^{-n}}{x_n}+\beta[/tex]
Plug this to your equation and choose [itex]\alpha,\beta[/itex] in order to arrive to [itex]x_{n+1}=x_n+b_n[/itex]

How do I get that to simplify into the required form? I tried plugging it in like this:

[tex]\frac{ \alpha ^{-n-1} }{x_{n+1} }+ \beta =\left( \frac{1}{x_{0} }+ \beta \right)\left[ \frac{3\left( \frac{ \alpha ^{-n} }{x_{n} }+ \beta \right)+\frac{1}{x_{0} }+ \beta }{2\left( \frac{ \alpha ^{-n} }{x_{n} }+ \beta \right)+\frac{1}{x_{0} }+ \beta } \right][/tex]

but I couldn't figure out how to simplify it without being left with a [tex]x_n x_{n+1}[/tex] term.
 
You don't change [itex]c_0[/itex] into [itex]\frac{1}{x_0}+\beta[/itex] because it is a constant. Just write

[tex]\frac{ \alpha ^{-n-1} }{x_{n+1} }+ \beta =c_0 \frac{3\left( \frac{ \alpha ^{-n} }{x_{n} }+ \beta \right)+c_0}{2\left( \frac{ \alpha ^{-n} }{x_{n} }+ \beta \right)+c_0}[/tex]

and calculate [itex]\alpha,\beta[/itex]. No term [itex]x_n\,x_{n+1}[/itex] must survive.
 
Here's what I have so far:

After making the substitution you suggested, I solved for [tex]\beta[/tex] by noticing that in order for the [tex]x_n x_{n+1}[/tex] term to disappear, either [tex]\alpha=0[/tex] or [tex]c_0^2 + 2c_0\beta -2\beta^2 = 0[/tex]. Solving, I got [tex]\beta=\pm \frac{c_0}{2} (1+\sqrt{3}).[/tex] For convenience, I only used [tex]\beta = \frac{c_0}{2} (1+\sqrt{3}).[/tex] Am I allowed to do that?

So, substituting in, I got
[tex]x_{n+1} = \left ( \frac{2+\sqrt{3}}{2-\sqrt{3}} \right ) \frac{x_n}{\alpha} + \frac{2}{c_0\alpha^{n+1} (2-\sqrt{3})}.[/tex]​

This is where I get "stuck". I figured out a way to solve for [tex]x_n[/tex] that's analagous to finding an integrating factor for a linear ODE, but it involves really nasty algebra. Is there a better way?
 
Last edited:
You are almost there! :smile:
Choose for
[tex]\alpha=\frac{2+\sqrt{3}}{2-\sqrt{3}}[/tex]
in order to eliminate the coefficient of [itex]x_n[/itex]. Can you solve the resulting equation?