Add a derivative function to SolutionFunctor and use it in ODEs