Tenemos:
$$\int_{0}^{1}e^{x^2}\,dx = e\cdot\int_{0}^{1}e^{1-x^2}\,dx = e\cdot\int_{0}^{1}e^{x(2-x)}\,dx=2e\cdot\int_{0}^{1/2}e^{4x(1-x)}\,dx$$
y puesto que:
$$ e^z = \sum_{j=0}^{+\infty}\frac{z^j}{j!} $$
tenemos:
$$\begin{eqnarray*}\int_{0}^{1}e^{x^2}\,dx &=& 2e\cdot\int_{0}^{1/2}\sum_{j=0}^{+\infty}\frac{(-1)^j(4x(1-x))^j}{j!}\,dx =2e\cdot\sum_{j=0}^{+\infty}\frac{(-4)^j B(1/2,j+1,j+1)}{j!}\\&=&\frac{e}{2}\sum_{j=1}^{+\infty}\frac{(-1)^{j+1}4^j}{j!\binom{2j}{j}}=\frac{e\sqrt{\pi}}{2}\sum_{j=0}^{+\infty}\frac{(-1)^j}{\Gamma(j+3/2)}.\end{eqnarray*}$$
(ver esta pregunta acerca de una identidad propia de la función Beta incompleta)
Sumando los términos de la a a $j=5$ en la última serie, se obtiene:
$$1.469\ldots = \frac{73}{135}e\geq \int_{0}^{1} e^{x^2}\,dx \geq \frac{207}{385}e = 1.461\ldots. $$