Testing the extremes of initial mass function variability using compact stellar systems: Code from MSc Thesis project at the University of Waterloo, 2020-2022