Numerical estimate of finite size corrections to the free energy of the SK model using Guerra--Toninelli interpolation

I use an interpolating formula introduced by Guerra and Toninelli to investigate numerically the finite size corrections to the free energy of the Sherrington--Kirkpatrick model. The results are compatible with a $(1/12 N) \ln(N/N_0)$ behavior at $T_c$, as predicted by Parisi, Ritort and Slanina, and a $1/N^{2/3}$ behavior below $T_c$.