Uno podría tenga en cuenta que:
$$\frac{n!}{(2n)!}=\frac{\sqrt\pi}{4^n\Gamma(n+1/2)}$$
De hecho, esto hace que su serie es un caso especial de la función siguiente:
$$f_\alpha(x)=\sum_{n=0}^\infty\frac{x^n}{\Gamma(n+\alpha)}$$
Con $S=\sqrt\pi f_{1/2}(1/4)$.
Maravillosamente, este problema sirve como una buena introducción para las fracciones de los derivados. Vamos a definir aquí los siguientes:
$$\frac{\mathrm d^\alpha}{\mathrm dx^\alpha}f(x)=\frac1{\Gamma(1-\{\alpha\})}\frac{\mathrm d^{\lceil\alpha\rceil}}{\mathrm dx^{\lceil\alpha\rceil}}\int_0^x(x-t)^{-\{\alpha\}}f(t)~\mathrm dt$$
donde $\lceil\alpha\rceil$ es el límite máximo de la función y $\{\alpha\}=\alpha-\lfloor\alpha\rfloor$ es la parte fraccionaria de $\alpha$. De la definición anterior, se puede deducir que
$$\frac{\mathrm d^\alpha}{\mathrm dx^\alpha}x^\beta=\frac{x^{\beta-\alpha}\Gamma(\beta+1)}{\Gamma(\beta-\alpha+1)}$$
A partir de esto, uno puede deducir que
$$\frac{x^n}{\Gamma(n+\alpha)}=x^{1-\alpha}\frac{\mathrm d^{1-\alpha}}{\mathrm dx^{1-\alpha}}\frac{x^n}{n!}$$
Por lo tanto, nos encontramos con que
$$f_\alpha(x)=x^{1-\alpha}\frac{\mathrm d^{1-\alpha}}{\mathrm dx^{1-\alpha}}e^x$$
En el caso particular de $\alpha=1/2$,
$$\frac{\mathrm d^{1/2}}{\mathrm dx^{1/2}}e^x=\frac1{\Gamma(1/2)}\int_0^x(x-t)^{-1/2}e^t~\mathrm dt$$
Deje $x-t=u^2$ para obtener
$$\frac{\mathrm d^{1/2}}{\mathrm dx^{1/2}}e^x=\frac1{\Gamma(3/2)}\int_0^{\sqrt x}u^2e^{x-u^2}~\mathrm du$$
Con lo cual la relación con la función de error es inmediata.
El caso general para $f_\alpha(x)$ queda como crédito adicional para el lector con la sugerencia adicional que Gamma incompleta funciones son recomendados.