Exact numerical implementation of the background (i.e., 1D) collisional Boltzmann equation for DM-baryon elastic scattering. Code used to generate results presented in S. Seher Gandhi & Y. Ali-Haïmoud (2022) [arXiv:2205.05536]. Please cite our paper if you use this code in whole or in part.