We have:
$$\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$$
and since:
$$ e^z = \sum_{j=0}^{+\infty}\frac{z^j}{j!} $$
we have:
$$\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*}$$
(see this question about an identity for the incomplete Beta function)
By summing the terms up to $j=5$ in the last series we get:
$$1.469\ldots = \frac{73}{135}e\geq \int_{0}^{1} e^{x^2}\,dx \geq \frac{207}{385}e = 1.461\ldots. $$