There is no simple formula for $u_n$ in terms of $a$. Indeed, $u_n$ is a polynomial expression on $a$ of degree $2^{n-1}$.
If you seek a simple formula for $u_n$ to understand its convergence, then this is hopeless because convergence is explained by Julia sets, which are quite complicated. In this case, the Julia set of $f(z)=z^2+z$:
There is no simple formula for $u_n$ in terms of $a$. Indeed, $u_n$ is a polynomial expression on $a$ of degree $2^{n-1}$.
If you seek a simple formula for $u_n$ to understand its convergence, then this is hopeless because convergence is explained by Julia sets, which are quite complicated. In this case, the Julia set of $f(z)=z^2+z$:
(Image computed with WA)