function [h,hh] = neostrict(P,m,n)
h = 0;
hh = [];
for t=1:n-1
  hh(end+1) = strict(P,m,n,t);
end
h = max(hh);
return

function h = strict(P,m,n,t)
h = 0;
dimt = repmat([m],[1 t]);
for xi = 1:m^t
  Xt0 = ind2coord(xi,dimt);
  Xt1 = Xt0(1:t-1);
  
  P0 = getcondP(P,m,n,t+1,1:t  ,Xt0);
  P1 = getcondP(P,m,n,t+1,1:t-1,Xt1);  
  
  hxi = .5*sum(abs(P0-P1));
  h = max(h,hxi);
end
return