function checkmxcondexp(m,n)
%checkmxcondexp(P,F,m,n)


F = fillsupsuma(m,n)


while 1
P=unifP(m,n)
  
for i=1:n
%for i=2

  di = 0;
  for yi = 1:m^(i-1)
    y = i2c(yi,m,i-1);
    for w1=1:m
      for w2=w1+1:m
        %fprintf('yi=%d w1=%d w2=%d  \n',yi-1,w1-1,w2-1);
        %Ef1 = condexphat(P,F,y,m,n,w1);
        %Ef2 = condexphat(P,F,y,m,n,w1);
        Ef1 = condexp(P,F,[y w1],m,n);
        Ef2 = condexp(P,F,[y w2],m,n);
        %fprintf('Ef1 = %1.5f;   Ef2 = %1.5f  \n',Ef1,Ef2);        
        di = max(di,abs(Ef1-Ef2));
      end
    end
  end

  mxdi = mxvi(P,F,m,n,i);
  err = abs(di-mxdi)
  if err>1e-9
    keyboard
  end
end
end
return



function i=c2i(c,m,n)
dims = repmat(m,[1 n]);
i = coord2ind(c,dims);
return

function c=i2c(i,m,n)
dims = repmat(m,[1 n]);
c = ind2coord(i,dims);
return