function [spk_vect,v_vect] = inhom_pp(dt,stim_vect)
%
%[spk_vect,v_vect] = inhom_pp(dt,stim_vect)
%
%dt is the time step in  msec
%stim_vect is an instantaneous firing frequency in spk/sec 
%  
  
  v_vect = zeros(size(stim_vect));
  spk_vect = zeros(size(stim_vect));

  v_thres = exprnd(1);
  v(1) = 0;
  rho_vect = (stim_vect/1000); %in spk/msec

  %save some time
  dtrho_vect = dt*rho_vect;

  for i = 2:length(stim_vect)
    v(i) = v(i-1) + dtrho_vect(i);
    if (v(i) > v_thres)
      v(i) = 0;
      spk_vect(i) = 1;
      v_thres = exprnd(1);
    end;
  end;