I wrote my model code in MATLAB because it is convenient and fast. You are welcome to download and use the scripts as you like. They may also work with Octave. If you translate them to another language I'd be happy to share a link or host your files here — just get in touch.

This is the main script that repeatedly generates model data by integrating out the differential equations to see how well it “fits” with empirical data. It compares multiple fits and adjusts the parameters until it eventually finds the best fit.

• testfitode5.m

These fitted parameter values and their uncertainties are then used in the following scripts to interpret the model.

These scripts generate the graphs on these pages:

• plotr0firstderiv.m
• plotcfrh.m
• plotcases.m
• plotlogcases.m
• ultimateimpactlatest.m

These scripts grab the latest empirical data from online sources and compile them into timeseries data my model can use:

• loadMontanaflynnData.m
• loadWikipediaData.m
• loadCmriversData.m
• loadJonathanDuganData.m