If the initial configuration is strictly defined such that \(p(y,0) = \delta (y-y_0)\) , where \(\delta (\cdot )\)  is the Dirac delta, the solution of (REF ) is given by the following ODE, see p. 53–54 [1]}:
\(\dot{y} \equiv \frac{\mathrm {d}y}{\mathrm {d}t}=F(y)-y,\quad y(0)=y_0\;.\) 
