Consider this picture. Integrating over this infinitesimal box gives the following equivalencies:
$$\int_{\Delta V} d^3r~{\rm div} \vec{E}(\vec{r}) = \int_{S(\Delta V)} d\vec{f} \cdot \vec{E}(\vec{r}),$$
which for $\Delta x \rightarrow 0$ gives $\Delta F \vec{n} \cdot (\vec{E_a} - \vec{E_i})$. Since this integral has to be equal to $\frac{1}{\epsilon_0}\sigma \Delta F$ we can deduce that $$\vec{n} \cdot (\vec{E_a}-\vec{E_i}) = \frac{\sigma}{\epsilon_0}.$$
My question is, how do I calculate this limit of $\Delta x \rightarrow 0$? I get that the integral basically computes all the scalar products of the electric field on the surface and adds them up, but how do I get the vector $ (\vec{E_a}-\vec{E_i})$? Why is the electric field given by this vector? It's probably really simple but I can't make sense of this right now, so some help would be appreciated.