clear
dt=0.0001;
t=0:dt:5;
ntrials=10;
spikes=zeros(ntrials,length(t));
for j = 1:ntrials
    j
    for i=1:length(t)
        rate = 5+10*t(i);
        if (rand < dt*rate) 
            spikes(j,i) = 1;
        end
    end
    subplot(ntrials,1,j)
    j
    plot(t,spikes(j,:))
end