function itr=itr(N,P,Dur);

if  P== 1
  itr = log2(N)*60/Dur;
else
  itr = (log2(N)+P*log2(P)+(1-P)*log2((1-P)/(N-1)))*60/Dur;
end