The code for generating figures 12 and 13 is provided.
The ".py" file performs the experiment, and the jupyter notebook file processes the outputs and generates plots.
The Mathematica code to compute the roots of the cubic polynomial in eq. (83) is implemented in "q3q4.nb", and the result id used the .py file.