function MI = maxpairmutinfo(P,m,n)

MI = 0;
for i=1:m
  for j=i+1:m
    MI = max(MI,mutualinfo(P,m,n,i,j));
  end
end
return