clear
dt=0.0001;
t=0:dt:60;
rate = 10;
spikes=zeros(size(t));
for i=1:length(t)
    if (rand < dt*rate) 
        spikes(i) = 1;
    end
end