Generalize oderkf to support orders != 5.
Add ode23_bs based on oderkf.
Add performance comparison for adaptive solvers.
Some cosmetics & new test added (B5).
Rebased. This gives the same tolerances in ode23 and oderkf.