I posted the same question to Math Stackexchange earlier without much luck, so I am posting here. I have a space $\mathcal{F}$ of vector-valued functions $f(\in\mathcal{F}):\mathbb{R}^d \rightarrow \mathbb{R}^n$. Consider a function $\Pi: \mathcal{F} \rightarrow \mathcal{F}$: $$ \Pi (f) (x) = \sum_{i=1}^{P} K \left(x,x_i\right) f(x_i) $$ where $x\in\mathbb{R}^d$, and $K:\mathbb{R}^d \times \mathbb{R}^d \rightarrow \mathbb{R}^{n\times n}$. To clarify, $\Pi (f) (x)$ means "evaluate $\Pi (f) \in \mathcal{F}$ at $x$". Our dynamical equation for a function $f$ is: $$ \frac{\partial f_t}{\partial t}=\Pi(f_t) $$ $t$ is a time in $\mathbb{R}$. Apparently, assuming $f_0$ as the initial condition, the solution is $$ f_t = \exp(t\Pi) f_0$$ but I am not sure how to prove this. The solution for a matrix differential equation has the same form and I can prove it for that, but I am not sure how to do it for this functional case. To provide the context, this is from the neural tangent kernel paper by Jacot et al.