% function np=interbin(x,bins,nt)%% calculate the number of interpoint distances of the embedded time% series in the bins specified by bins.% I.e.%  np(k) is the number of pairs of points x(:,i) and x(:,j)%  (with |j-i|>nt) such that%                bins(k-1)<||x(:,i)-x(:,j)||<bins(k)           %  [assuming bins(0)=0]%% nt is the Theiler window (default nt=0)% bins are the k bin sizes (must be ascending)% x is a d-by-n matrix: the d-dimensional embedding of n-points%%% For more info, read README%% Michael Small% ensmall@polyu.edu.hk% 26/2/02