Mean field self consistency single site.py

From Werner KRAUTH

(Difference between revisions)
Jump to: navigation, search
Revision as of 10:59, 8 November 2016
Werner (Talk | contribs)

← Previous diff
Revision as of 16:01, 8 November 2016
Werner (Talk | contribs)

Next diff →
Line 21: Line 21:
pylab.savefig('mean_field_self_consistency_Ising.png') pylab.savefig('mean_field_self_consistency_Ising.png')
pylab.show() pylab.show()
 +
 +[[Category:ICFP Lectures]]

Revision as of 16:01, 8 November 2016

import random, math, pylab

dim = 2
q = 2.0 * dim
Tc = 2.0 * dim
m_vec = [1.0]
T_vec = [0]
for iter in range(-1000, 1000):
    t = iter / float(1100.0)
    beta = (t * Tc + Tc) ** (-1)
    T_vec.append(1.0 / beta)
    m = 0.1
    for iter2 in range(1000):
        m = math.tanh(beta * q * m)
    m_vec.append(m)
pylab.plot(T_vec, m_vec)
pylab.title('Weiss self-consistency for the Ising model 
in $d$ = ' +  str(dim) + ' dimensions')
pylab.xlabel('$T$ (Temperature)', fontsize=18)
pylab.ylabel('$m$ (magnetization per site)', fontsize=18)
pylab.savefig('mean_field_self_consistency_Ising.png')
pylab.show()
Personal tools