Holman et al., 1999
This is the current method used to discern the stability of moons around binary planets. It's implemented in the Yakei-Odius System spreadsheet. It seems based off this, we can have stable moons past ~100,000km, avoiding n:1 resonances.
Lam et al., 2018
This paper provides a Python script written in Python 2 that requires a few adjustments to get working. Parker et al., 2021 mentions this but does not use it, instead relying on the methods described in Holman et al., 1999.
Parker et al., 2021
This program may be useful to try and simulate stability.