Initial public release of a code to calculate the luminosity function of Milky Way (MW) satellite galaxies using Approximate Bayesian Computation. The analysis uses as input only the observed satellite galaxies in the SDSS and DES surveys, and five Dark Matter-only simulations of MW-like haloes from the Aquarius simulation suite (Aq-A - Aq-F). Input data of the correct format is supplied for all known confirmed and candidate satellite galaxies as of February 2018.