'initial_integrator' may now be provided as an ordinary function handle that computes the exact solution.