Storing stationary points during solving

Just use a continuous callback on integrator(t,Val{1}) to get the derivatives and then save to some array in the affect!