function X = Poisson(lambda0)
% POISSON produces a Poisson distributed 
% random number X (e.g. an 'inter-spike interval') at rate lambda0

  X = -1./lambda0*log(1-rand);