solve_henon_eqs.c