% make plot
clf

b = 4; % block to depict

% plot block PMF, error bars, and PMF from using all data
subplot(2,1,1);
% PMF using all data
plot(xi_i, F_i, 'k-');
hold on
errorbar(xi_i, F_ib(:,b), dF_ib(:,b), dF_ib(:,b), 'k.');
legend('20 ns sim', '2 ns block');
legend boxoff;
oldaxis = axis;
axis([-180 +180 oldaxis(3:4)]);
xlabel('\psi (degrees)');
ylabel('F(\psi) (kcal/mol)');
title('explicit solvent');

% plot bar chart comparing estimated error from block std
subplot(2,1,2);


dF_i_blockstd = std(F_ib')';

bar(xi_i, [dF_ib(:,b) dF_i_blockstd]);
legh = legend('estimated', 'std over blocks', 'Location', 'NorthWest');
legend boxoff;
oldaxis = axis;
axis([-180 +180 oldaxis(3:4)]);
xlabel('\psi (degrees)');
ylabel('uncertainty (kcal/mol)');