This repository contains two jupyter notebooks that reproduce the figures in Goulding+ 2019.
This makes Fig. 3 in the paper, which shows the timescales for J1010+1413 as a function of semi-major axis.
This makes Fig. 5 in the paper, which shows our number density assumptions and final GWB amplitude calculations.
Feel free to use these notebooks to calculate timescales or limits for your favorite object! If you have any questions or comments on this code, email me at: kpardo @ astro . princeton . edu.