¿Qué acerca de la cuadratura? Tenemos
$$ \left(\sum_{n\geq 0}x^{n^2}\right)^{2} = \sum_{n\geq 0}r_2(n)\,x^n $$
donde
$$ r_2(n) = \left|\left\{(a,b)\in\mathbb{N}^2: a^2+b^2 = n\right\}\right| $$
es una buena función aritmética. Por ejemplo, $\sum_{n=0}^{N}r_2(n)$ es el número de celosía puntos con los no-negativo coordenadas en el interior del círculo $x^2+y^2=N$, pero también el coeficiente de $x^N$ en
$$ \frac{1}{1-x}\left(\sum_{n\geq 0}x^{n^2}\right)^{2} = \sum_{N\geq 0}\left(\sum_{n=0}^{N}r_2(n)\right) x^N. $$
Por Gauss círculo problema con Voronoi obligado tenemos que $\sum_{n=0}^{N}r_2(n)=\frac{\pi N}{4}+O(N^{1/3})$, por lo tanto $x=1$ es un doble polo de la función anterior y el primer término de su Laurent de expansión alrededor de $x=1$ está dado por $\frac{\pi}{4(1-x)^2}$. Multiplicando por $(1-x)^2$ y tomando el límite cuando $x\to 1^-$ obtenemos:
$$ \lim_{x\to 1^-}\,(1-x)\left(\sum_{n\geq 0}x^{n^2}\right)^2 = \frac{\pi}{4} $$
y el reclamo fácilmente de la siguiente manera.
Es interesante señalar que el mismo enfoque muestra que, para cualquier $k\geq 2$,
$$ \lim_{x\to 1^-}(1-x)\left(\sum_{n\geq 0}x^{n^k}\right)^k $$
es exactamente la (hiper-)el volumen encerrado por $x_1\geq 0,x_2\geq 0,\ldots x_k\geq 0$$x_1^k+x_2^k+\ldots+x_k^k = 1$, que puede ser calculada de la siguiente manera: si asumimos que el $X_1,\ldots,X_k$ son yo.yo.d. variables aleatorias con distribución uniforme sobre $(0,1)$, el PDF de $X_i$ es compatible en $(0,1)$ y dado por $f(t)=\frac{1}{k} t^{1/k-1}$. Para calcular el volumen anterior, es suficiente para el cálculo de la función característica de a $X_i$, tome su $k$-ésima potencia, considere la inversa de la transformada de Fourier y la integran más de $(0,1)$ para obtener la probabilidad de que $X_1^k+\ldots+X_k^k\leq 1.$