load('ampullary.mat')
fig = figure();

%% create PSTH on the basis of the interspike intervals
fig.sub
% 1. get the interspike intervals for each trial
for i = 1:size(times,2)
   isi = diff(times{i});
    
end

%% create PSTH using the binning method



%% create PSTH using the kernel-convolution method