function [d,dd] = sumsqmdif(P,F,m,n)

dd = zeros(1,n);
d=0;
for i=1:n
  di = mxvi(P,F,m,n,i);
  d = d + di^2;
  dd(i) = di;
end

return

% OLD school
function d = sumsqmdif0(P,F,m,n)

d = 0;
for i=1:n
  di = 0;
  for yi = 1:m^i
    y1 = i2c(yi,m,i);
    y0 = y1(1:end-1);
    Ef1 = condexp(P,F,y1,m,n);
    Ef0 = condexp(P,F,y0,m,n);
    di = max(di,abs(Ef1-Ef0));
  end
  d = d + di^2;
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