Aquí hay otra solución. Después del cambio de variable $x = -y$$y \to 1^-$, tenemos
$$ f(x) = \sum_{n=0}^{\infty} y^{4n^2}(1 - y^{4n+1}). \tag{1}$$
Ahora la idea es identificar a $\text{(1)}$ como una especie de suma de Riemann para la integral $\int_{0}^{\infty} e^{-2t} \, dt = \frac{1}{2}$.
Para ello, vamos a la primera revisión arbitraria $\epsilon \in (0, 1)$ y elija $N = N(\epsilon, y) \in \mathbb{Z}_{> 0}$ tal que $y^{4N^2} \to \epsilon$$y \to 1^-$. Por ejemplo, elija $N$ a ser cualquier entero más cercano a $(\log \epsilon / 4 \log y)^{1/2}$. A continuación, mediante la alternancia de comportamiento de $\text{(1)}$ nos encontramos con que
$$ \left| \sum_{n=N+1}^{\infty} y^{4n^2}(1 - y^{4n+1}) \right| \leq y^{4N^2} \xrightarrow[y \to 1^-]{} \epsilon $$
que se le muestra al ser insignificante, ya que la $\epsilon \to 0^+$. A continuación, vamos a definir $t_n = t_n(y)$ por
$$t_n = \sum_{k=0}^{n} (4k+1)(1-y) = (2n^2+3n+1)(1-y) $$
y $\Delta t_n = \Delta t_n (y) := t_n - t_{n-1}$. Luego tenemos a $y = 1 - \frac{t_n}{2n^2+3n+1}$. Ahora, a partir de la siguiente desigualdad
$$ (4n+1)y^{4n+1}(1-y) \leq 1 - y^{4n+1} \leq (4n+1)(1-y), $$
nos encontramos con que
\begin{align*}
&\sum_{n=0}^{N} \exp\left\{ 4(n+1)^2 \log\left( 1 - \frac{t_n}{2n^2+3n+1} \right) \right\} \Delta t_n \tag{2}\\
&\qquad \leq \sum_{n=0}^{N} y^{4n^2}(1 - y^{4n+1}) \\
&\qquad \leq \sum_{n=0}^{N} \exp\left\{ 4n^2 \log\left( 1 - \frac{t_n}{2n^2+3n+1} \right) \right\} \Delta t_n \tag{3}
\end{align*}
Desde $t_N \to \frac{1}{2}\log\frac{1}{\epsilon}$$\max_{0 \leq n \leq N} \Delta t_n = (4N+1)(1-y) \to 0$$y \to 1^-$, podemos pensar en los límites $\text{(2)}$ $\text{(3)}$ como sumas de Riemann de $e^{-2t}$ sobre el inteval $[0, \frac{1}{2}\log\frac{1}{\epsilon}]$. Tomando como límite $y \to 1^-$, no es difícil demostrar (pero un poco molesto para dar un argumento) que ambos límites convergen en el mismo valor, que es
$$ \lim_{y \to 1^-} \sum_{n=0}^{N} y^{4n^2}(1 - y^{4n+1}) = \int_{0}^{\frac{1}{2}\log\frac{1}{\epsilon}} e^{-2t} \, dt = \frac{1 - \epsilon}{2}. $$
A continuación, por la costumbre de 3 epsilon argumento, tenemos
$$ \limsup_{y \1^-} \left| f(-y) - \frac{1}{2} \right| \leq \frac{3 \epsilon}{2}
\xrightarrow[\epsilon \to 0^+]{} 0 $$
y de ahí el reclamo de la siguiente manera.