trajectories_boundaries_combined.m produces figure 1

resilience_estimate_combined.m produces figure 2 using the data found in resilience_estimate_matrixDD_high_res.mat (which is produced by resilience_estimate_matrixDD.m using the function in resilienceEstimateDD.m) and resilience_estimate_matrixP_high_res.mat (which is produced by resilience_estimate_matrix_P.m using the function in resilienceEstimateP.m)

simulation_results_combined.m produces figure 3 using the data found in spkcountoutput_high_res.mat (which is produced by spk_crashcounter.m) and sddkcountoutput_high_res.mat (which is produced by sddk_crashcounter.m)

Several of the scripts use the function in dxdtAllee.m