
T = 0.001:0.01:5;
beta = 0:0.01:5;
EX =  ( 1 - (1./ (sinh(beta).^4) ) ).^(1/8);
EX = EX.*(beta>0.88);
plot(beta,EX)
axis([0,5,-0.5,1.5]);
xlabel('Inverse Temperature');
ylabel('Mean Field Value');

print -deps meanplot.eps

