function isivec = isis( spikes )
% returns a single list of isis computed from all trials in spikes
%
% isivec = isis( spikes )
%   spikes: a cell array of vectors of spike times in seconds
%   isivec: a column vector with all the interspike intervalls

    isivec = [];
    for k = 1:length(spikes)
        difftimes = diff( spikes{k} );
        % difftimes(:) ensures a column vector
        % regardless of the type of vector spikes{k}
        isivec = [ isivec; difftimes(:) ];
    end
end