Limit of sequence given by $x_{n} = \frac{x_{n-2}+x_{n-3}}{2}$? Let $x_1,x_2,x_3\in\mathbb{R}$ be three distinct real numbers. I am interested in the convergence of the sequence
$$
x_{n} = \frac{x_{n-2}+x_{n-3}}{2},\quad n = 4,5,\ldots
$$
ie,
$$
x_{4} = \frac{x_{2}+x_{1}}{2},\quad x_{5} = \frac{x_{3}+x_{2}}{2},\quad x_{6} = \frac{x_{4}+x_{3}}{2},\quad\cdots
$$
Please note that, although related to, this is not a recursive averaging sequence with two initial values, as described and studied elsewhere --- see, for instance, Smith, Scott G. "Recursive Averaging". The Mathematics Teacher, Vol. 108, No. 7 (March 2015), pp. 553-557.
In the present case, convergence is easily verified by noting that
$$
\big|\,x_{n}-x_{n-1}\,\big| = \big|\frac{x_{n-2}+x_{n-3}}{2}-\frac{x_{n-3}+x_{n-4}}{2}\big| = \big|\,\frac{1}{2} (x_{n-2}-x_{n-4})\,\big| = \big|\,\frac{1}{2^{\,n-4}}\,\big|\cdot\big|\, x_{3}-x_{1}\,\big| 
$$
Therefore the sequence is Cauchy and hence it converges in $\mathbb{R}$.
However, I was not able to find the limit to where the sequence converges. Computer experiments carried out with several different sets of initial values suggest that this limit should be a linear combination of the those initial values, as occurs in the case mentioned above. Also, the results of these experiments point to the existence of (how many?) sub-sequences, all converging to the very same limit.
I would appreciate some help in the analysis of this problem.

The characteristic equation  of the linear recurrence is
$$\lambda^3 - \frac{1}{2}\lambda-\frac{1}{2} = 0$$ with roots $\lambda_1=1$, $\lambda_{2,3} = \frac{1}{\sqrt{2}}\cdot(-\frac{1}{\sqrt{2}} \pm \frac{i}{\sqrt{2}})$. There exist unique $a$, $b$, $c$ such that for all $n \ge 0$ we have
$$x_n = a \cdot 1^n + b \cdot \lambda_2^n + c \cdot \lambda_3^n$$
Since $|\lambda_{2,3}|= \frac{1}{\sqrt{2}}< 1$, we have
$$\lambda_2^n \to 0, \ \ \lambda_3^n \to 0$$
for $n \to \infty$. We conclude that $x_n \to a$. Now we only need to obtain $a$ from the equalities
$$x_0 = a + b + c \\
x_1 = a + b \lambda_2 + c \lambda_3\\
x_2 = a + b \lambda_2^2 + c \lambda_3^3$$
We get  with Cramer's rule
$$ a = \frac{\left| \begin{matrix} x_0 & 1 & 1 \\ x_1 & -\frac{1}{2} + \frac{i}{2} & -\frac{1}{2} - \frac{i}{2}\\ x_2 & (-\frac{1}{2} + \frac{i}{2})^2 & (-\frac{1}{2} - \frac{i}{2})^2 \end{matrix} \right |} {\left| \begin{matrix} 1 & 1 & 1 \\ 1 & -\frac{1}{2} + \frac{i}{2} & -\frac{1}{2} - \frac{i}{2}\\ 1 & (-\frac{1}{2} + \frac{i}{2})^2 & (-\frac{1}{2} - \frac{i}{2})^2 \end{matrix} \right |} = \frac{1}{5}(x_0 + 2 x_1 + 2 x_2)$$
$\bf{Added:}$ Since we are dealing with linear recurrences it's worth getting comfortable with determinants of the form
$$D(x_0, x_1, \ldots x_{n-1}) \colon = \left| \begin{matrix} x_0 & 1 & \ldots &1 \\
                        x_1 & \mu_2& \ldots& \mu_{n} \\
                        x_2 & \mu_2^2& \ldots & \mu_{n}^2 \\
                        \ldots& & & \ldots \\
                        x_{n-1} & \mu_2^{n-1} & \ldots & \mu_{n}^{n-1} \end{matrix} \right |$$
Note that when $x_k = x^k$ for $0\le k \le n-1$ we have a Vandermonde determinant, which equals
$$V(x, \mu_2, \ldots, \mu_{n}) = V( \mu_2, \ldots, \mu_{n}) \cdot \prod_{i=2}^{n}( \mu_i - x)= \\ =(-1)^{n-1} V( \mu_2, \ldots, \mu_{n})\cdot (x-\mu_2)\cdot (x-\mu_{n}) = (-1)^{n-1} V( \mu_2, \ldots, \mu_{n})\cdot ( x^{n-1} - s_1 x^{n-2} + s_2 x^{n-3} - \cdots)$$
It follows that the determinant with first column $(x_0, x_1, \ldots x_{n-1})$ equals
$$(-1)^{n-1} V(\mu_2, \ldots, \mu_{n+1}) \cdot ( x_{n-1} - s_1 x_{n-2} + s_2 x_{n-3} - \cdots)$$
Now when solving the system we have a quotient of such determinants $D$, the extra factor $(-1)^{n-1} \cdot V$ can be omitted. We get: the coefficient in front of a power $\mu_1^n$ equals
$$\frac{P(x)/(x-\mu_1)}{P'(\mu_1)}$$
where the numerator is considered "in the umbral sense".
Example with our problem at hand:
The characteristic polynomial equals $P(x)=x^3 - \frac{1}{2} x - \frac{1}{2}$.

*

*For the numerator : divide $P(x)$ by $x-\mu_1= x-1$ and get
$\frac{1}{2}(2 x^2 + 2 x + 1)$. This umbrally is $\frac{1}{2}(2 x_2 + 2x_1 + x_0)$.


*For the denominator we have $P'(x) = 3 x^2 - \frac{1}{2}$ so $P'(1) = \frac{5}{2}$


*Set up the fraction
$$a_1 = a =\frac{\frac{1}{2}(2 x_2 + 2 x_1 + x_0)}{\frac{5}{2}} = \frac{2 x_2 + 2 x_1 + x_0}{5}$$
