Update ordinary_differential_eq.c