solve_duffing_eqs.c