How do I evaluate $\int \limits_{-\infty}^{a} e^{−t^2}dt$? I know that
$$I \equiv \int \limits_{-\infty}^\infty e^{−t^2} \, dt=\sqrt{\pi},\text{ and }\int \limits_{-\infty}^0 e^{−t^2} \, dt=\frac{\sqrt{\pi}}{2}.$$
However, I don't understand if (or how) I can find a similar solution for $\int \limits_{-\infty}^{a}e^{−t^2}dt, a \neq 0$, given that the error function actually does not yield closed form solutions.
Any help is greatly appreciated!

Since $e^{-t^2}$ is an even function, and since you know $\int_{-\infty}^{\infty}$, if you were able to find $\int_{-\infty}^b$ then you would be able to find $\int_a^b$ for any $a,b \in \mathbb{R}$. The fact that $e^{-t^2}$ does not have an elementary primative should suggest to you that you probably can't find $\int_a^b$ explicitly. 
As an A-level student (UK 16-18 pre-university), I always wondered why we had to use a table of normal distribution values to get approximate probabilities. Now I know that the probability density function is basically a stretched and translated version of $e^{-t^2}$.
Food for thought: the natural logarithm, $\ln x$, has an integral definition:
$$ \ln x := \int_1^x \frac{dt}{t}  $$
because $t^{-1}$ has no known primitive, yet we call $\ln x$ an elementary function, while $\text{erf} \, x$ is not.
