Encontré un posible enfoque numérico compatible con Gabriel Romon arriba. Calcular $a_n$ para $n=0,1,2,...,10$.
Definir $q_n=(3-a_n)6^{n}$.
Luego, $\ln(1-\frac{q_n}{q_{n+1}})$ puede ajustarse con una línea recta como función de $n$ teniendo $r^2=0.9983$.
Si se toma el argumento del logaritmo natural aquí como la razón de la primera derivada a la función misma, esto implica que $q(n)$ es una doble exponencial.
Es decir:
$$q_n=(3-a_n)6^n \approx q(n)= c_1e^{(\frac{e^b}{m}e^{mn})}$$
donde m y b son la pendiente y la intersección respectivamente para la mencionada línea, y $c_1$ es una constante.
En mis modelos, m=-1.8368, b = -0.4174.
Esto sugiere que $c_1$ es el límite en cuestión y ajustar $q(n)$ contra $q_n$ producirá $c_1$, aproximando el resultado deseado, sin embargo hasta ahora mis intentos dan un valor de aproximadamente 4.27.