Numerical and statistical analysis of NeuralODE with Runge-Kutta time integration