function t=tau(data);

[entropy]=mutual(data);
for i = 1:length(entropy)-1           
    if (entropy(i)<=entropy(i+1))
        t = i-1;            
        break;
    end
end
               
% plot(0:length(entropy)-1,entropy)
% xlabel('tau');
% title('互信息法求tau');
%