Given a sequence of functions $\{f_k\}$, suppose for any $k\geq 4$,
$$ f_k=f_1f_{k-1}+\frac{1}{2}(f_2-f_1^2)f_{k-2}+(\frac{1}{6}f_1^3-\frac{1}{2}f_1f_2+\frac{1}{3}f_3)f_{k-3}. $$
I want to obtain the general explicit expression of $f_k$ in terms of $f_1, f_2, f_3$. Is it possible?
I have tried the method that $$ f_k-Af_{k-1}-Bf_{k-2}=\lambda(f_{k-1}-Af_{k-2}-Bf_{k-3}) $$ but find that the degree 3 equation of $\lambda$ cannot be solved (always with some square root or cubic root)... Hence I failed.
This is a classical linear recurrence, and the characteristic polynomial is indeed of the third degree.
Unless some simplification occurs in the coefficients, by luck, you can't escape the complexity of Cardano's formulas. You should try to put the polynomial in the depressed form.