function kij=BIC(Tc)

% Tc is a vector

n=length(Tc);
kij=zeros(n);
kij(1:end-1,end)=1-((2*sqrt(Tc(end)*Tc(1:end-1)'))./(Tc(end)+Tc(1:end-1)')).^0.27;
kij(end,1:end-1)=1-((2*sqrt(Tc(end)*Tc(1:end-1)))./(Tc(end)+Tc(1:end-1))).^0.27;

end