Changed from 2 programs to 5 by adding a run_hayford_model,
hayford_model_old, and hayford_model_modified_parameters (to change
around parameters and initial conditions), and updated hayford_model and
hayford_model_stochastic