Problems

---

Solution:

1. \(\,\) \begin{align*} && u_{n+2} - u_n &= 1 + \frac{1}{u_{n+1}} - \left (1 + \frac{1}{u_{n-1}} \right) \\ &&&= \frac{1}{1 + \frac1{u_n}} - \frac{1}{1 + \frac{1}{u_{n-2}}} \\ &&&= \frac{u_n}{u_n+1} - \frac{u_{n-2}}{1+u_{n-2}} \\ &&&= \frac{u_n(1+u_{n-2}) - u_{n-2}(1+u_n)}{(1+u_n)(1+u_{n-2})} \\ &&&= \frac{u_n - u_{n-2}}{(1+u_n)(1+u_{n-2})} \end{align*}
2. Claim: \(u_n \in [1,2]\) Proof: (By induction). Note that \(u_1 = 1, u_2 = 2\) so our claim is true for the first few terms. Note that if \(u_n \in [1,2]\), \(\frac{1}{u_n} \in [\tfrac12, 1]\) and \(1+\frac{1}{u_{n}} \in [\tfrac32,2] \subset [1,2]\). Therefore \(u_{n+1} \in [1,2]\). Therefore since \(u_1 \in [1,2]\) and \(u_n \in [1,2] \Rightarrow u_{n+1} \in [1,2]\) \(u_n \in [1,2]\) for all \(n \ge 1\)
3. First notice that \(u_3 = \frac32 > u_1\) and therefore by the recursion we found in the first part, \(u_{2n+1}-u_{2n-1} > 0\) so \(u_{2k+1}\) is increasing and bounded, and so by our theorem converges to a limit. Suppose this limit is \(L\), then we must have \(L = 1 + \frac1{L} \Rightarrow L^2 - L - 1 = 0 \Rightarrow L = \frac{1+\sqrt5}{2}\) since it must be in \([1,2]\). Similarly, not that \(u_4 = \frac{5}{3} < u_2\) and so \(u_{2k+2} - u_{2k} < 0\) and \(-u_{2k}\) is increasing and bounded above. Therefore it tends to a limit (and so does \(u_{2k}\)). By the same reasoning as before, it's the same limit, \(\frac{1+\sqrt5}{2}\) and therefore the sequence converges. If \(u_1 = 3, u_2 = \frac43 \in [1,2]\) so we have our sequence being bounded and all the same logic will follow through.

---

Solution: Let \(t = \tan \frac12 x\), then \begin{align*} \frac{\d t}{\d x} &= \tfrac12 \sec^2 \tfrac12 t \\ &= \tfrac12 (1 + \tan^2 \tfrac12 ) \\ &= \tfrac12 (1 + t^2) \\ \\ \sin x &= 2 \sin \tfrac12 x \cos \tfrac12 \\ &= \frac{2 \frac{\sin \tfrac12 x}{ \cos \tfrac12x}}{\frac{1}{\cos^2 \tfrac12 x}} \\ &= \frac{2 \tan \tfrac12 x}{\sec^2 \tfrac12 } \\ &= \frac{2t }{1+t^2} \end{align*} Now consider \begin{align*} t = \tan \tfrac12 x: && \int_0^{\frac12\pi} \frac 1{1+a \sin x}\, \d x &= \int_{t=0}^{t = 1} \frac{1}{1 + a \frac{2t}{1+t^2}} \frac{2}{1+t^2} \d t \\ &&&= \int_0^1 \frac{2}{1+2at+t^2} \d t \\ &&&= \int_0^1 \frac{2}{(t+a)^2 + 1-a^2} \d t \\ (1-a^2) > 0: &&&= \left [ \frac{2}{\sqrt{1-a^2}} \arctan \frac{t+a}{\sqrt{1-a^2}} \right]_0^1 \\ &&&= \frac{2}{\sqrt{1-a^2}} \left ( \arctan \frac{1+a}{\sqrt{1-a^2}} - \arctan \frac{a}{\sqrt{1-a^2}} \right) \\ &&&= \frac{2}{\sqrt{1-a^2}} \arctan \left ( \frac{\frac{1+a}{\sqrt{1-a^2}}-\frac{a}{\sqrt{1-a^2}}}{1+\frac{1+a}{\sqrt{1-a^2}}\frac{a}{\sqrt{1-a^2}}} \right) \\ &&&= \frac{2}{\sqrt{1-a^2}} \arctan \left ( \frac{\frac{1}{\sqrt{1-a^2}}}{\frac{1+a}{1-a^2}} \right) \\ &&&= \frac{2}{\sqrt{1-a^2}} \arctan \left ( \frac{\sqrt{1-a}}{\sqrt{1+a}} \right) \end{align*} as required. Let \[ I_n = \int_0^{\frac12\pi} \frac{ \sin ^nx}{2+\sin x} \, \d x \qquad \quad (n\ge0). \] and consider \begin{align*} I_{n+1} + 2I_n &= \int_0^{\frac12\pi} \frac{ \sin ^{n+1}x+2\sin^{n} x}{2+\sin x} \, \d x \\ &= \int_0^{\frac12\pi} \frac{ \sin^n x (2 + \sin x)}{2+\sin x} \, \d x \\ &= \int_0^{\frac12\pi} \sin^n x \d x \end{align*} Therefore we can compute \begin{align*} I_0 &= \int_0^{\pi/2} \frac{1}{2 + \sin x} \d x \\ &= \frac12 \int_0^{\pi/2} \frac{1}{1 + \frac12 \sin x} \d x \\ &= \frac{1}{\sqrt{3/4}} \arctan \frac{\sqrt{1/2}}{\sqrt{3/2}} \\ &= \frac{2}{\sqrt{3}} \arctan \frac{1}{\sqrt{3}} \\ &= \frac{\pi}{3\sqrt{3}} \\ \\ I_1 &= \int_0^{\pi/2} 1 \d x - 2 I_0 \\ &= \frac{\pi}{2} - \frac{2\pi}{3\sqrt{3}} \\ I_2 &= \int_0^{\pi/2} \sin x \d x - 2I_1 \\ &= 1 - \pi + \frac{4\pi}{3\sqrt{3}} \\ I_3 &= \int_0^{\pi/2} \sin^2 x \d x - 2I_2 \\ &= \frac12 \int_0^{\pi/2} \sin^2 + \cos^2 x \d x - 2I_2 \\ &= \frac{\pi}{4} - 2 + 2\pi - \frac{8\pi}{3\sqrt{3}} \\ &= -2 + \frac{9\pi}{4} - \frac{8\pi}{3\sqrt{3}} \end{align*}

---

Solution: \begin{align*} && y &= \frac{\arcsin x}{\sqrt{1-x^2}} \\ && \frac{\d y}{\d x} &= \frac{(1-x^2)^{-1/2} \cdot (1-x^2)^{1/2}-\arcsin x \cdot (-x)(1-x^2)^{-1/2}}{1-x^2} \\ &&&= \frac{1+ xy}{1-x^2} \\ \Rightarrow && 0 &= (1-x^2) \frac{\d y}{\d x} -xy-1\\ \\ \frac{\d^n}{\d x^{n+1}}: && 0 &= \left ( (1-x^2) y' \right)^{(n+1)} - (xy)^{(n+1)} \\ \Rightarrow && 0 &= (1-x^2)y^{(n+2)} + \binom{n+1}{1}(1-x^2)^{(1)}y^{(n+1)}+\binom{n+1}{2} (1-x^2)^{(2)}y^{(n)} - (xy^{(n+1)} +\binom{n+1}{1} y^{(n)} ) \\ &&&= (1-x^2)y^{(n+2)}+\left ( (n+1)\cdot(-2x)-x \right)y^{(n+1)} + \left ( \frac{(n+1)n}{2} \cdot (-2)-(n+1) \right)y^{(n)} \\ &&&= (1-x^2)y^{(n+2)}-\left ( 2n+3 \right)xy^{(n+1)} - \left ( (n+1)n+(n+1)\right)y^{(n)} \\ &&&= (1-x^2)y^{(n+2)}-\left ( 2n+3 \right)xy^{(n+1)} - \left ( n+1\right)^2y^{(n)} \\ \end{align*} Since \(y(0) = 0, y'(0) = 1\) we can look at the recursion: \(y^{(n+2)} - (n+1)^2y^{(n)}\) for larger terms, ie \(y^{(2k)}(0) = 0\) \(y^{(1)}(0) = 1, y^{(3)}(0) = (1+1)^2 \cdot 1 = 2^2, y^{(5)}(0) = (3+1)^2 y^{(3)} = 4^2 \cdot 2^2\) and \(y^{(2k+1)}(0) = (2k)^2 \cdot (2k-2)^2 \cdots 2^2 \cdot 1^2 = 2^{2k} \cdot (k!)^2\). Therefore \begin{align*} && \frac{\arcsin x}{\sqrt{1-x^2}} &= \sum_{k=0}^{\infty} \frac{2^{2k} \cdot (k!)^2}{(2k+1)!} x^{2k+1} \\ \\ \Rightarrow && \frac{\arcsin \frac12}{\sqrt{1-\left (\frac12 \right)^2}} &= \sum_{k=0}^{\infty} \frac{2^{2k} \cdot (k!)^2}{(2k+1)!} 2^{-2k-1}\\ &&&= \frac12 \sum_{k=0}^{\infty} \frac{ (k!)^2}{(2k+1)!} \\ &&&= \frac12 \left ( 1 + \frac1{3!} + \frac{2^2}{5!} + \cdots+ \right) \\ \Rightarrow&& S &= 2 \frac{2\frac{\pi}{6}}{\sqrt{3}} = \frac{2\pi}{3\sqrt{6}} \end{align*}

---

Solution:

1. Suppose \(t_n = x\) for all \(n\), then we must have \begin{align*} && x &= p x + q x \\ \Leftrightarrow && 1 &= p+q \end{align*} and this clearly works for any value of \(x\).
2. Suppose \(t_{2n} = x, t_{2n+1} = y\) for all \(n\), then \begin{align*} && x &= py + q x \\ && y &= px + q y \\ \Rightarrow && 0 &= py + (q-1) x \\ && 0 &= px + (q-1) y \\ \Rightarrow && p &= (q-1) \frac{x}{y} = (q-1) \frac{y}{x} \\ \Rightarrow && \frac{y}{x} = \pm 1 & \text{ or } q = 1, p = 0 \\ \Rightarrow && y = \pm x & \text{ or } (p,q) = (0,1) \\ \end{align*}
3. Suppose \(t_{3n} = x\), \(t_{3n+1}=y\) and \(t_{3n+2}=z\) , so \begin{align*} && x &= pz + qy \\ && y & = px + qz \\ && z &= py + qx \\ \\ && z &= p(px+qz) + q(pz + qy) \\ &&&= p^2x + 2pqz + q^2 y \\ &&&= p^2(pz+qy) + 2pqz + q^2(px+qz) \\ &&&= p^3 z + p^2qy + 2pqz + q^2p x + q^3 z \\ &&&= (p^3+q^3+2pq)z + pq(py+qx) \\ &&&= (p^3 + q^3 + 2pq)z + pq z \\ &&&= (p^3 + q^3 + 3pq)z \\ \Rightarrow && 0 &= p^3 + q^3 + 3pq- 1 \\ &&&= (p+q-1)(p^2+q^2+1+p+q-pq) \\ &&&= \tfrac12(p+q-1)((p-q)^2+(p+1)^2+(q+1)^2) \end{align*} Therefore \(p+q = 1\) or \((p-q)^2+(p+1)^2+(q+1)^2 = 0 \Rightarrow p = q = -1\). If \(p+q = 1\), then \(z = py + (1-p)x\) and \(x = p(py+(1-p)x) + (1-p)y \Rightarrow (1-p+p^2)x = (1-p+p^2)y \Rightarrow x = y \Rightarrow x= y = z\). If \(p = q = -1\) then adding all the equations we get \(x + y + z = -2(x+y+z) \Rightarrow x + y + z = 0\)

Note that what is actually going on here is that solutions must be of the form \(t_n = \lambda^n\) so the only way to be constant is for \(\lambda = 1\) to be a root, the only way for it to be \(2\)-periodic is for \(\lambda = -1\) to be a root, and the only way for it to be \(3\)-periodic is for \(\lambda = 1, \omega, \omega^2\) to be the roots (although we see this via the classic \(x^3 + y^3 + z^3 - 3xyz = (x+y+z)(x + \omega y + \omega^2 z)(x+\omega^2 y +\omega z)\) which is because of the real constraint in the question.

---

Solution: \begin{align*} && \beta - \alpha &> q \\ \Rightarrow &&(\beta - \alpha)^2 &> q^2 \\ \Rightarrow && \beta^2 +\alpha^2 - 2\beta \alpha &> q^2 \\ \Rightarrow && \alpha^2+\beta^2-q^2 -2 \beta \alpha &> 0 \\ \Rightarrow && \frac{\alpha^2+\beta^2-q^2}{\alpha\beta} - 2 &> 0 \end{align*} \begin{align*} && u_n(u_n+u_{n+2}) &= u_n \cdot \left (u_n + \frac {u_{n+1}^2 -q^2}{u_{n}}\right) \\ &&&= u_n^2 + u_{n+1}^2-q^2 \\ &&&= u_n^2 + u_{n+1}^2 - (u_n^2-u_{n-1}u_{n+1}) \\ &&&= u_{n+1}^2 + u_{n+1}u_{n-1} \\ &&&= u_{n+1}(u_{n-1}+u_{n+1}) \\ \\ && u_{n+1}-pu_n+u_{n-1} &= -pu_n+\frac{u_{n}(u_{n-2}+u_n)}{u_{n-1}} \\ &&&= \frac{u_n(u_{n}-pu_{n-1}+u_{n-2})}{u_{n-1}} \end{align*} Therefore if \(u_2 -pu_1 + u_0 = 0\) it is always zero, ie if \begin{align*} && u_2 &= p\beta - \alpha \\ && u_2 &= \frac{\beta^2-q^2}{\alpha} \\ \Rightarrow && \frac{\beta^2-q^2}{\alpha} &= p\beta - \alpha \\ \Rightarrow && p &= \frac{\alpha^2+\beta^2-q^2}{\alpha\beta} \end{align*} If \(\beta > \alpha + q\) we must have that \(p > 2\), and so \(u_{n+1}-u_n = (p-1)u_n - u_{n-1} > u_n-u_{n-1} > 0\), therefore the sequence is strictly increasing. If \(\beta = \alpha + q\) the sequence follows \(u_{n+1} - 2u_n + u_{n-1} =0\) and so \(u_{n+1}-u_n = u_n - u_{n-1}\) for all \(n\) (which is still increasing - it's an arithmetic progression with common difference \(\beta - \alpha\)).

---

Solution:

1. \(\,\) \begin{align*} && 0 &= a+b \tag{1}\\ && 1 &= a\lambda -a\mu \tag{2} \\ && 1 &= a\lambda^2 -a\mu^2 \tag{3} \\ && 2 &= a\lambda^3 - a\mu^3 \tag{4} \\ (4) \div (2): && 2 & = \lambda^2+\lambda \mu + \mu^2 \\ (3) \div (2): && 1 &= \lambda + \mu \\ \Rightarrow && 2 &= \lambda^2 + \lambda(1-\lambda) + (1-\lambda)^2 \\ &&&= \lambda^2-\lambda+1\\ \Rightarrow && \lambda, \mu &= \frac{1 \pm \sqrt{5}}{2} \\ \Rightarrow && a &= \frac{1}{\lambda - \mu} = \frac{1}{\sqrt{5}} \\ \Rightarrow && b &= -\frac{1}{\sqrt{5}} \end{align*} (NB: This is Binet's formula)
2. \(\,\) \begin{align*} F_6 &= \frac{1}{\sqrt{5}} \left ( \left ( \frac{1 +\sqrt{5}}{2} \right)^6- \left ( \frac{1 -\sqrt{5}}{2} \right)^6 \right) \\ &= \frac{1}{2^6 \sqrt{5}} \left ( (1+\sqrt{5})^6-(1-\sqrt{5})^6 \right) \\ &= \frac{1}{2^5 \cdot \sqrt{5}} \left (6 \sqrt{5} +\binom{6}{3} (\sqrt{5})^3+\binom{6}{5}(\sqrt{5})^5 \right)\\ &= \frac{1}{2^5} \left (6 +20\cdot 5+6\cdot 5^2 \right)\\ &= \frac{1}{2^5} 256 = 2^3 = 8 \end{align*} (way more painful than just computing it by adding terms!)
3. \(\,\) \begin{align*} && \sum_{n=0}^{\infty} \frac{F_n}{2^{n+1}} &= \sum_{n=0}^{\infty} \frac{a\lambda^n + b\mu^n}{2^{n+1}} \\ &&&= \frac12 \left ( \frac{a}{1-\frac{\lambda}2} + \frac{b}{1-\frac{\mu}2} \right) \\ &&&= \frac12 \left ( \frac{2a}{2-\lambda} + \frac{2b}{2-\mu}\right) \\ &&&= \frac{2a}{4-2\lambda} + \frac{2b}{4-2\mu}\\ &&&= \frac{2a}{4-(1+\sqrt{5})} - \frac{2a}{4-(1-\sqrt{5})} \\ &&&= \frac{2}{3\sqrt{5}-5} - \frac{2}{3\sqrt{5}+5} \\ &&&= \frac{6\sqrt{5}+10-6\sqrt{5}+10}{45-25} \\ &&&= 1 \end{align*}

---

Solution: \begin{align*} && y &= \cos(m \arcsin x) \\ && y' &= -m \sin (m \arcsin x) \cdot (1-x^2)^{-\frac12} \\ && y'' &= -m^2 \cos(m \arcsin x) \cdot (1-x^2)^{-1} -m \sin(m \arcsin x) \cdot (1-x^2)^{-\frac32} \cdot (-x) \\ &&&= -m^2 y (1-x^2)^{-1} + x(1-x^2)^{-1} y' \\ \Rightarrow && 0 &= (1-x^2)y'' - x y' + m^2y \\ \\ && 0 &= (1-x^2)y^{(3)} -2xy'' - xy''-y' + m^2y' \\ &&&= (1-x^2)y^{(3)} - 3xy'' + (m^2-1)y' \\ \\ && 0 &= (1-x^2)y^{(4)} - 2xy^{(3)} - 3xy^{(3)} - 3y^{(2)} + (m^2-1)y^{(2)} \\ &&&= (1-x^2)y^{(4)}- 5xy^{(3)} - (m^2-4)y^{(2)} \end{align*} Claim: \(0 = (1-x^2)y^{(n+2)} - (2n+1)y^{(n+1)} + (m^2-n^2)y^{(n)}\) Proof: (By induction) Clearly the first few base cases are true. Suppose it is true for some \(n\), then \begin{align*} && 0 &= (1-x^2)y^{(n+2)} - (2n+1)xy^{(n+1)} + (m^2-n^2)y^{(n)} \\ \Rightarrow && 0 &= (1-x^2)y^{(n+3)} - 2xy^{(n+2)} - (2n+1)xy^{(n+2)} - (2n+1)y^{(n+1)} + (m^2-n^2)y^{(n+1)} \\ &&&= (1-x^2)y^{(n+3)} - (2n+3)xy^{(n+2)} + (m^2-n^2-2n-1)y^{(n+1)} \\ &&&= (1-x^2)y^{(n+1+2)} - (2(n+1)+1)xy^{(n+1+1)} +(m^2-(n+1)^2)y^{(n)} \end{align*} And so we can conclude the result by induction. Notice that \begin{align*} && y(0) &= \cos(m 0) = 1 \\ && y'(0) &= -m\sin(m 0) = 0 \\ && y''(0) &= -m^2 y(0) = -m^2\\ \end{align*} Notice that \(y^{(n+2)}(0) + (m^2-n^2)y^{(n)} = 0\) so in particular all the odd terms will be \(0\) and the even terms will be \(1, -m^2, m^2(m^2-2^2), \cdots\), therefore \begin{align*} && \cos (m \arcsin x) &= 1 -\frac{m^2}{2!} x^2 + \frac{m^2(m^2-2^2)}{4!}x^4 - \cdots \\ \Rightarrow && \cos(m \theta) &= 1 - \frac{m^2}{2!} \sin^2 \theta + \frac{m^2(m^2-2^2)}{4!} \sin^4 \theta \end{align*} Notice that if \(m\) is even, then at some point we will have \(m^2-m^2\) appearing in our expansion and all remaining terms will be zero. Therefore we will end up with a polynomial series, of degree \(m\) in \(\sin \theta\)

---

Solution: \begin{array}{c|r} n & F_n \\ \hline 1 & 1 \\ 2 & 1 \\ 3 & 2 \\ 4 & 3 \\ 5 & 5 \\ 6 & 8 \\ 7 & 13 \\ 8 & 21 \\ 9 & 34 \\ 10 & 55 \end{array} \begin{questionparts} \item Claim: \(\frac1{F_i} > \frac1{2F_{i-1}}\) for \(i \geq 4\). Proof: Since \(F_i = F_{i-1}+F_{i-2}\) and \(F_i > 1\) for \(i \geq 1\) we have \(F_i > F_{i-1}\) for \(i \geq 3\). In particular we have \(F_i = F_{i-1}+F_{i-2} < 2F_{i-1}\) for \(i -1 \geq 3\) or \(i \geq 4\). Therefore \(\frac{1}{F_i} > \frac1{2F_{i-1}}\)

---

Solution:

1. Let \(\displaystyle y= \sum_{n=0}^\infty a_n x^n\,\) then \begin{align*} y' &= \frac{\d}{\d x} \l \sum_{n=0}^\infty a_n x^n \r \\ &= \sum_{n=0}^\infty \frac{\d}{\d x} \l a_n x^n \r \\ &= \sum_{n=0}^\infty n a_n x^{n-1} \\ &= \sum_{n=1}^\infty n a_n x^{n-1} \\ \\ y'' &= \frac{\d}{\d x} \l\sum_{n=1}^\infty n a_n x^{n-1} \r \\ &= \sum_{n=1}^\infty \frac{\d}{\d x} \l n a_n x^{n-1} \r \\ &= \sum_{n=1}^\infty n(n-1) a_n x^{n-2} \\ &= \sum_{n=2}^\infty n(n-1) a_n x^{n-2} \\ \\ y &= a_0 + a_1 x+ a_2x^2 + a_3x^3 + \cdots \\ y'&= a_1 + 2a_2x+3a_3x^2 + \cdots \\ y'' &= 2a_2 + 6a_3x + \cdots \end{align*}
2. \begin{align*} && 0 &= xy''-y'+4x^3y \\ &&&= x\sum_{n=2}^\infty n(n-1) a_n x^{n-2} - \sum_{n=1}^\infty n a_n x^{n-1} + 4x^3 \sum_{n=0}^\infty a_n x^n \\ &&&= \sum_{n=2}^\infty n(n-1) a_n x^{n-1} - \sum_{n=1}^\infty n a_n x^{n-1} + \sum_{n=0}^\infty 4a_n x^{n+3} \\ &&&= \sum_{n=2}^\infty n(n-1) a_n x^{n-1} - \sum_{n=1}^\infty n a_n x^{n-1} + \sum_{n=4}^\infty 4a_{n-4} x^{n-1} \\ &&&= \sum_{n=4}^{\infty} \l n(n-1) a_n- n a_n +4a_{n-4} \r x^{n-1} + 2a_2x + 6a_3x^2-a_1-2a_2x-3a_3x^2 \\ &&&= \sum_{n=4}^{\infty} \l n(n-2) a_n +4a_{n-4} \r x^{n-1}+ 3a_3x^2-a_1 \\ \end{align*} Therefore since all coefficients are \(0\), \(a_1 = 0\), \(a_3 = 0\) and \(\displaystyle a_n = -\frac{4}{n(n-2)}a_{n-4}\). If \(a_0 = 1, a_2 = 0\), and since \(a_1 = 0, a_3 = 0\) the only values which will take non-zero value are \(a_{4k}\). We can compute these values as: \(a_{4k} = -\frac{4}{(4k)(4k-2)} a_{4k-4} = \frac{1}{2k(2k-1)}a_{4k-r}\) so \(a_{4k} = \frac{(-1)^k}{(2k)!}\), which are precisely the coefficients in the expansion \(\cos x^2\). If \(a_0 = 0, a_2 = 1\) then since \(a_1 = 0, a_3 = 0\) the only values which take non-zero values are \(a_{4k+2}\) we can compute these values as: \(a_{4k+2} = -\frac{4}{(4k+2)(4k)}a_{4k-2} = -\frac{1}{(2k+1)2k}a_{4k-2}\) so we can see that \(a_{4k+2}= \frac{(-1)^k}{(2k+1)!}\) precisely the coefficients of \(\sin x^2\)

---

Solution: \begin{align*} && (x+y+z)^2 &= x^2 + y^2 + z^2 + 2(xy+yz+zx) \\ \Rightarrow && 1^2 &= 2 + 2(xy+yz+zx) \\ \Rightarrow && xy+yz+zx &= -\frac12 \end{align*} \begin{align*} && 1 \cdot 2 &= (x+y+z)(x^2+y^2+z^2) \\ &&&= x^3 + y^3 + z^3 + x^2y+x^2z+y^2z+y^2x+z^2x+z^2y \\ &&&= 3 + x^2y+x^2z+y^2z+y^2x+z^2x+z^2y\\ \Rightarrow && -1 &= x^2y+x^2z+y^2z+y^2x+z^2x+z^2y \end{align*} \begin{align*} && (x+y+z)^3 &= x^3 + y^3 + z^3 + \\ &&&\quad \quad 3xy^2 + 3xz^2 + \cdots + 3zx^2 + 3zy^2 + \\ &&&\quad \quad \quad 6xyz \\ \Rightarrow && 1 &= 3 + 3(-1) + 6xyz \\ \Rightarrow && xyz &= \frac16 \end{align*} Since we have \(f(t) = (t-x)(t-y)(t-z) = t^3-t^2-\frac12 t - \frac16\) is zero for \(x,y,z\) we can notice that: \(t^{n+1} = t^n +\frac12 t^{n-1} + \frac16 t^{n-2}\) is also true for \(x,y,z\) (by multiplying by \(t^{n-2}\). Therefore: \(S_{n+1} = S_n + \frac12 S_{n-1} + \frac16 S_{n-2}\)

---

Solution: This is a second order recurrence relation, so we need to find \(m\) and \(n\) such that; \begin{align*} &&\frac15 &= m\frac13 + n \\ &&\frac17 &= m \frac15 + n\frac13 \\ \Rightarrow && m,n &= \frac67, - \frac{3}{35} \end{align*} So we now need to solve the characteristic equation: \(\lambda^2 - \frac67 \lambda + \frac{3}{35} = 0\) So \(x,y = \frac{15 \pm 2 \sqrt{30}}{35}\). We need, \begin{align*} && 1 &= a+ b \\ && \frac13 &= a \frac{15 + 2 \sqrt{30}}{35} + b \frac{15 - 2 \sqrt{30}}{35} \\ && \frac13 &= \frac{15}{35} + \frac{2 \sqrt{30}}{35}(a-b) \\ \Rightarrow && -\frac{\sqrt{30}}{18} &= a-b \\ \Rightarrow && a &= \frac{18-\sqrt{30}}{36} \\ && b &= \frac{18+\sqrt{30}}{38} \end{align*} So our two answers are: \[ (a,b,x,y) = \left (\frac{18\pm\sqrt{30}}{36} ,\frac{18\mp\sqrt{30}}{36},\frac{15 \pm 2 \sqrt{30}}{35},\frac{15 \mp 2 \sqrt{30}}{35}, \right)\]