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 and merged.