Add solver ode23s for stiff systems (based on modified Rosenbrock triple) and
make fdjacobian available to all solvers.