[A2,Plaq,Q_top,Wloop,theta2,stat]=HMC_KS([],0);
theta1=theta2;
[A2,Plaq,Q_top,Wloop,theta2,stat]=HMC_KS(theta1,1);
theta1=theta2;
[A2,Plaq,Q_top,Wloop,theta2,stat]=HMC_KS(theta1,1);
theta1=theta2;
[A2,Plaq,Q_top,Wloop,theta2,stat]=HMC_KS(theta1,1);
theta1=theta2;
[A2,Plaq,Q_top,Wloop,theta2,stat]=HMC_KS(theta1,1);
theta1=theta2;
[A2,Plaq,Q_top,Wloop,theta2,stat]=HMC_KS(theta1,1);
theta1=theta2;
[A2,Plaq,Q_top,Wloop,theta2,stat]=HMC_KS(theta1,1);
auto=Autocorel(Wloop(:,1),20);
semilogy(auto)
hold
auto=Autocorel(Wloop(:,2),20);
semilogy(auto)
auto=Autocorel(Wloop(:,3),20);
semilogy(auto)
auto=Autocorel(Wloop(:,4),20);
semilogy(auto)
auto=Autocorel(Wloop(:,5),20);
semilogy(auto)
xlabel('t')
ylabel('Autocorel')
replot
gset terminal postscript
gset out 'Auto.ps'
replot
gset terminal x11
hold
w=mean(Wloop);
sw=max(Binning(Wloop(:,1),10));
sw=[sw,max(Binning(Wloop(:,2),10))];
sw=[sw,max(Binning(Wloop(:,3),10))];
sw=[sw,max(Binning(Wloop(:,4),10))];
sw=[sw,max(Binning(Wloop(:,5),10))];
sw=[sw,max(Binning(Wloop(:,6),10))];
sw=[sw,max(Binning(Wloop(:,7),10))];
sw=[sw,max(Binning(Wloop(:,8),10))];
sw=[sw,max(Binning(Wloop(:,9),10))];
sw=[sw,max(Binning(Wloop(:,10),10))];
semilogyerr(w(1:5),sw(1:5))
hold
axis([0,6,1e-4,1])
semilogy(w(1:5))
axis([0,6,1e-3,1])
xlabel('Linear Size')
ylabel('Wloop')
replot
gset terminal postscript
gset out 'Wloop.ps'
replot
gset terminal x11
hold
hist(Q_top,30)
xlabel('Q_top')
ylabel('Frequency')
gset terminal postscript
gset out 'Q_top.ps'
replot
gset terminal x11
hold