ising_mf_t.py¶
from __future__ import print_function
import numpy as np
from ising_mf import find_solution, calc_entropy, calc_specific_heat
def main():
h = 0
tmin = 0
tmax = 1
n = 201
for t in np.linspace(tmin, tmax, n):
try:
m = find_solution(1., t, h)
s = calc_entropy(m)
c = calc_specific_heat(m, 1., t)
except ValueError:
print(f"# {t:.4f}")
else:
print(f"{t:.4f} {m:.5e} {s:.5e} {c:.5e}")
if __name__ == '__main__':
main()