function isis = isi(data,range)
% ISI列を返す関数
% ISI列 = isi(スパイク時刻列, 試行 [開始 終了] 配列)
% ISI列 = isi(スパイク時刻列, スカラー) 

isis=[];

if size(range,2) == 2
  trial=size(range,1);		%%% 試行数
  len=range(1,2)-range(1,1);	%%% 切り出す長さ

  for i=1:trial
    train=data(range(i,1)<=data & data<range(i,2));
			%%% 試行ごとに切り出す．
    tmp= [train; 0] - [0; train]; 	%%% 一個ずらして、引く
    isis=[isis; tmp(2:length(train))];	%%% 最初と最後以外をISI列に追加
  end
else
    tmp= [data; 0] - [0; data]; 	%%% 一個ずらして、引く
    isis=tmp(2:length(data));	%%% 最初と最後以外をISI列に追加
end
