Numerical ODE integration in C