function T=TriCorr(x)
%
% TriCorr Triple Correlation function 1-D
%
% Example T(d1,d2)=TriCorr(x);
%

M=length(x);
T=zeros(M,M);
for d1=0:M-1
   for d2=d1:M-1
      for i=0:M-1-d2
	 T(d1+1,d2+1)=T(d1+1,d2+1)+x(i+1)*x(i+1+d1)*x(i+1+d2);
      end;
      T(d1+1,d2+1)=T(d1+1,d2+1)/(M-d2);
   end;
end;
for d2=0:M-1
   for d1=d2+1:M-1
      for i=0:M-1-d1
	 T(d1+1,d2+1)=T(d1+1,d2+1)+x(i+1)*x(i+1+d1)*x(i+1+d2);
      end;
      T(d1+1,d2+1)=T(d1+1,d2+1)/(M-d1);
   end;
end;

