Reproduce results from: Burkardt & Trenchea (2020) (Q7260366)

A **Snakemake** workflow that reproduces the numerical experiments of Burkardt and Trenchea (2020) comparing explicit and fixed-point implicit midpoint ODE solvers. Both solvers are implemented from scratch in **NumPy** following the pseudocode in the paper and applied to two test systems: the humps ODE (a stiff scalar problem) and a predator-prey Lotka-Volterra system.\N\N**Output**: Trajectory plots and phase-plane diagrams for each solver–problem combination, ensuring exact reproducibility of all six output figures.\N\N### Main steps\N\N**run_midpoint_explicit** — runs the explicit midpoint solver on both ODEs, producing a summary text file and two PNG plots.\N\N**run_midpoint_fixed** — runs the implicit fixed-point solver on both ODEs, producing an equivalent set of three outputs. Both rules are independent and can run in parallel.