ode_solve function

solve ODEs