m=2,n=3
P=unifP(m,n)
U = prodmeas(P,m,n);

Py = getky(P,1,m,n);
Py = Py / sum(Py);
ytil = prodmeas(Py,m,n-1)

Uy = getky(U,1,m,n);
Uy = Uy / sum(Uy);
tily = Uy

totdif(ytil,tily)



if 0

%function tmp
aa = 0:.01:.99;
nn=1:100000;
xx=[];yy=[];
for a=aa
    xx(end+1)=a;
    yy(end+1) = sum(a.^sqrt(nn));
    plot(xx,yy)
    drawnow
end
%return
end