d=3, M=8, ||Hc|| = 2.0000, cota pi/(4||Hc||) = 0.3927, amax = 40.0, parámetros 180, semillas 8
  T = 0.4500/J (0.286 x pi/2)  ->  1-F = 5.27e-03  | scipy: 1-F = 5.27e-03, fuga = 1.4e-04, max ||A_k|| = 35.6 J