Update ordinary_differential_eq.h