I am trying to model a poisson process. When I histogram the inter-event times, the distribution is exponential. I want to prove that the simulation is Poissonian. How can I calculate the Fano-factor?I have a list of waiting times and spiketimes.

while duration <= Total_time

randNumbers = rand(1,1); IVI = -log(randNumbers)/rate;

SpikeTime = duration + IVI;

duration = duration + IVI; SpikeTimes = [SpikeTimes; SpikeTime];

counter = counter + 1;

IVI_hist(counter,1) = IVI;

end

figure(1) hist(IVI_hist, sqrt(length(IVI_hist)))