[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,'g-')
hold
auto=Autocorel(Wloop(:,2),20);
semilogy(auto,'r-')
auto=Autocorel(Wloop(:,3),20);
semilogy(auto,'b-')
auto=Autocorel(Wloop(:,4),20);
semilogy(auto,'k-')
auto=Autocorel(Wloop(:,5),20);
semilogy(auto,'m-')
xlabel('t')
ylabel('Autocorel')
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))];
plot(log(abs(w(1:5))));
hold
errorbar(log(abs(w(1:5))),sw(1:5)./abs(w(1:5)))
xlabel('Linear Size')
ylabel('LogWloop')
hold
hist(Q_top,30)
xlabel('Q_top')
ylabel('Frequency')
hold