Mathematica code for joint inference of IBD via MCMC