I want to understand the correspondence between locally constant sheaves and vector bundles with flat connection on a manifold $X$.
Given a local system $\mathcal L$, it is clear how to define a vector bundle $E$ with locally constant transition functions. This allows us to define a connection $\nabla$ on the sheaf of smooth sections of $L$, i.e. on $\mathcal L \otimes \mathcal C^\infty_X$ it by demanding that $\nabla$ vanishes on $\mathcal L$, i.e. locally $\nabla = d$.
This connection is clearly flat, as $d^2=0$, and $\mathcal L = \ker \nabla$.
But what about the other direction. Assume we have a vector bundle $E$ which admits a flat connection $\nabla$. I read, that then $\ker \nabla$ defines a local system $\mathcal L$, that is we consider the sheaf of solutions to the differential equation $\nabla s=0$
How do I show that $\mathcal L$ is locally constant?
For any $x_0\in X$, we can find a neighbourhood $U$, where $\nabla = d + A$ for some matrix $A$ of $1$-forms.
Thus there, a section $s$ we have$$s=(f_1,\dots, f_n) \in \ker \nabla ~~~~~~\Leftrightarrow ~~~~~ ~ df_i +\sum_j a_{ij} f_j=0, ~~\forall i.$$
Assuming that for any initial value $s(x_0)=(f_1(x_0),\dots,f_n(x_0))\in \mathbb C^n$, there exists a unique solution to the above differential equation, we get an isomorphism $\mathcal L(U) \to \mathbb C^n$, which is consistent with restrictions.
But how can I show this? In the texts I read, it was always said that it follows from Frobenius theorem and then the proof was finished without expanding on it. So here is my question:
Why does the differential equation $$df_i +\sum_j a_{ij} f_j=0, ~~\forall i$$ have a unique solution for any initial value $s(x_0)=(f_1(x_0),\dots,f_n(x_0))\in \mathbb C^n$.
Why does it follow from Frobenius theorem and how do I use flatness of the bundle?