source: ETALON/CLIO/Analysis/Data19July/KK/KK_phase.m @ 723

 function theta=KK_phase(Pt)
 len=length(Pt);
 n=1: len;
 
for m=1: len
q=2*m/pi*log(Pt(n)/Pt(m));
g=(m.*m-n.*n)';
lm=logical(ones([1  len]));
lm(m)=0;
Th(m)=q(lm)*(1./g(lm));
%     if (m~=1)&&(m~=len)
%         Th(m)=Th(m)+(Th(m-1)+Th(m+1))/2;
%     end
end
% lx=2:(len-1);
% Th(lx)=Th(lx)+(Th(lx-1)+Th(lx+1))/2;

 theta=[-Th,0,fliplr(Th(2:length(Th)))];
 
% m=1: len;
% q=2/pi*repmat(m,len,1)*log((Pt(n)')*(1./Pt(m)));
% g=(repmat(m.*m,len,1)-repmat((n.*n)',1,len));
% lm=~logical(diag(ones([1  len])));
% Th(m)=q(lm)*(1./g(lm));%sum for some dimension