% Plot 1D PMF in psi.
clear

% Read data
[bin, centers, f, df] = textread('pmf-psi.out');

% Plot PMF
clf;
hold on

fill([centers fliplr(centers)], [(f+df) fliplr(f-df)], 0.7 * [1 1 1]);
plot(centers, f, 'k-', 'LineWidth', 1.5);
axis([min(centers) max(centers) min(f-df) max(f+df)]);
xlabel('\psi (degrees)');
ylabel('PMF (kT)');

%% Write plot.
outfilename = 'pmf-psi';
exportfig(gcf, sprintf('%s.eps', outfilename), 'width', 3.25, 'height', 3.75);
unix(sprintf('epstopdf %s.eps', outfilename));