function lags = lag(data1,data2,range)
% 全てのスパイクの組の時間差 l_ij = t_i - s_j 列を返す関数
% 時間差列 = lag(スパイク時刻列1, スパイク時刻列2, 試行 [開始 終了] 配列)

lags=[];

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

for i=1:trial
  train1=data1(range(i,1)<=data1 & data1<range(i,2));
  train2=data2(range(i,1)<=data2 & data2<range(i,2));
			%%% 試行ごとに切り出す．


  for j=1:length(train2)
%    lags=[lags; (train1-train2(j))];
    tmp=(train1-train2(j));
    lags=[lags; tmp(-200<tmp & tmp<200)];
	   %%% data2のスパイク時刻を基準とした、data1のスパイクの相対時刻列
  end
end