Solving the XY Model Using Mixed Integer Optimization in Python

Author here. Glad you liked the post enough to submit it!