runMuFit2D.m was used to calculate results for figure 4 of paper.