fix(dsolvers/ode_solver.cpp): fix split function

Make sure that the splitted intervals are the subsets of input domain.

fix #61