%
% raster.m
%
% plot the firing diagram associated with iaf2net.c
%

load iaf2netout
n = length(iaf2netout);
dt = .04;

for j=1:n/7,

    t = j*dt;
    s = iaf2netout(1+(j-1)*7:j*7);
    sind = find(s>0.1);
    if ~isempty(sind),
       plot(t*ones(size(sind)),sind,'*');
       hold on
    end

end

axis([0 t+5*dt 0 8])
xlabel('t (ms)','fontsize',16)
ylabel('cell','fontsize',16)
hold off
