function mV = getmv(P,F,m,n)

mV = 0;
for ni = 1:n
  for yi = 1:m^ni-1
    y = i2c(yi,m,ni-1);
    %for yi = 1:m^ni
    y = i2c(yi,m,ni);
    for w1 = 1:m
      for w2 = w1+1:m
        [Ef1,K1] = condexp(P,F,[y w1],m,n);
        [Ef2,K2] = condexp(P,F,[y w2],m,n);
        %[Ef1,K1] = condexp(P,F,[y         ],m,n);    
        %[Ef2,K2] = condexp(P,F,[y(1:end-1)],m,n);
        mV = max(mV,abs(Ef2-Ef1));
      end
    end
  end
end


return


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