w = 0.1;
bins = 0.0:1.0:10.0;

counts{1} = spikecounts(poissonspikes, w);
counts{2} = spikecounts(pifouspikes, w);
counts{3} = spikecounts(lifadaptspikes, w);
titles = {'Poisson', 'PIF OU', 'LIF adapt'};
for k = 1:3
    subplot(1, 3, k);
    [h, b] = hist(counts{k}, bins);
    bar(b, h/sum(h));
    title(titles{k})
    xlabel('Spike count');
    xlim([-0.5, 8.5]);
    ylim([0.0 0.8]);
end
savefigpdf(gcf, 'spikecountshists.pdf', 20, 7);