Aqui les dejo una corrida de 32×32 sitios, con temperaturas de 1 a 5 con paso 0.1 y npre=2000 npasos=8000:ising_dat.tar

con este grafican: print_ising

IMPORTANTE… a veces surge un error al generar los vectores "uposta" y "cvposta" (solucion de Onsager). Si ese es el caso, simplemente usen los archivos homonimos en el zip para "loadear". (esta accion esta comentada en el codigo ising_print, solo hace falta descomentar, con ctrl+5)

La solucion de Onsager esta en terminos de integrales elipticas… no logré una manera simple de implementar estas, pero si un paper (https://arxiv.org/abs/1411.2495) donde explican como ponerlas en funcion de hipergeométricas, las cuales si pude implementar.

Finalmente aqui dejo un programita para visualizar la convergencia …convergencia_cumsum

Aqui la energia:

aqui la diferencia entre onsager y metropolis:

vemos que es maxima cerca de Tc. Esto se debe es donde mas le cuesta converger al algoritmo.

Aqui el cv:

y la magnetizacion:

Finalemente vemos como difieren para un sistema bien chico (4×4) donde claramente metropolis difiere muchisimo del cv analitico (correspondiente al limite termodinámico)