Prove the statement: $$a_n=\frac{(2^{n-1}+1)a_1-2^{n-1}+1}{2^{n-1}+1-(2^{n-1}-1)a_1}$$ for the given recursive sequence: $$a_{n+1}=\frac{3a_n-1}{3-a_n}. $$
My attempt:
Proof by induction:
(1) Base: $\tau(1)$. $$ \begin{align} a_2 &=\frac{(2+1)a_1-2+1}{2+1-(2-1)a_1} \\ &=\frac{3a_1-1}{3-a_1}. \end{align} $$
(2) Assumption: Let $$ a_n=\frac{(2^{n-1}+1)a_1-2^{n-1}+1}{2^{n-1}+1-(2^{n-1}-1)a_1} $$ hold for some $n\in\mathbb N$.
(3) Step $\tau(n+1)$: $$ \begin{align} a_n &= \frac{3\cdot\frac{(2^{n-1}+1)a_1-2^{n-1}+1}{2^{n-1}+1-(2^{n-1}-1)a_1}-1}{3-\frac{(2^{n-1}+1)a_1-2^{n-1}+1}{2^{n-1}+1-(2^{n-1}-1)a_1}} \\ &= \frac{3\cdot2^{n-1}a_1+3a_1-3\cdot2^{n-1}+3-2^{n-1}-1+2^{n-1}a_1-a_1}{3\cdot 2^{n-1}+3-3\cdot 2^{n-1}a_1+3a_1-2^{n-1}a_1-a_1+2^{n-1}-1} \\ &= \frac{(4\cdot 2^{n-1}+2)a_1-4\cdot 2^{n-1}+2}{4\cdot 2^{n-1}+2-(4\cdot 2^{n-1}+2)a_1} \\ &= \frac{(2^n+1)a_1-2^n+1}{2^n+1-(2^n+1)a_1} \end{align} $$
Is this correct and is there a more efficient method than induction?
Embed $a_n$ in a projective coordinate $[a_n,1]$ for each $n$. Then the given recursive relation may be written as
$$\left[\begin{array}{c}a_{n+1}\\ 1\end{array}\right]=\left[\begin{array}{cc}3&-1\\ -1&3\end{array}\right]\left[\begin{array}{c}a_n\\ 1\end{array}\right].$$
It follows by recursion that $$\left[\begin{array}{c}a_n\\ 1\end{array}\right]=A^{n-1}\left[\begin{array}{c}a_1\\ 1\end{array}\right],n\geq 1~\quad (1)$$ where $$A=\left[\begin{array}{cc}3&-1\\ -1&3\end{array}\right].$$ Diagonalize $A$ to get $$ \frac 12\left[\begin{array}{cc}1&1\\ 1&-1\end{array}\right]A\left[\begin{array}{cc}1&1\\ 1&-1\end{array}\right]=\left[\begin{array}{cc}2&0\\ 0&4\end{array}\right]$$ $$\Rightarrow A^{n-1}=\left[\begin{array}{cc}1&1\\ 1&-1\end{array}\right]\left[\begin{array}{cc}2^{n-1}&0\\ 0&2^{2n-2}\end{array}\right]\frac 12\left[\begin{array}{cc}1&1\\ 1&-1\end{array}\right]$$ $$=\left[\begin{array}{cc}2^{n-2}+2^{2n-3}&2^{n-2}-2^{2n-3}\\ 2^{n-2}-2^{2n-3}&2^{n-2}+2^{2n-3}\end{array}\right].$$ Substituting this back to (1), one gets that $$\left[\begin{array}{c}a_n\\ 1\end{array}\right]=\left[\begin{array}{cc}2^{n-2}+2^{2n-3}&2^{n-2}-2^{2n-3}\\ 2^{n-2}-2^{2n-3}&2^{n-2}+2^{2n-3}\end{array}\right]\left[\begin{array}{c}a_1\\ 1\end{array}\right]$$ $$\Leftrightarrow a_n=\frac{(2^{n-2}+2^{2n-3})a_1+(2^{n-2}-2^{2n-3})}{(2^{n-2}-2^{2n-3})a_1+(2^{n-2}+2^{2n-3})}=\frac{(2^{n-1}+1)a_1-2^{n-1}+1}{2^{n-1}+1-(2^{n-1}-1)a_1},$$ as required. QED