function convh
m=2; n=3;

while 1

P1 = pnormdim(rand(m^n,1),1);
P2 = pnormdim(rand(m^n,1),1);

[hh1,Ht1,H1] = gethhn(P1,m,n);
[hh2,Ht2,H2] = gethhn(P2,m,n);

hold on
good = 1;
for lam=0:.01:1
    P3 = lam*P1 + (1-lam)*P2;
    HL = lam*H1 + (1-lam)*H2;
    [hh3,Ht3,H3] = gethhn(P3,m,n);
    plot(lam,HL,'.b');
    plot(lam,H3,'.r');
    drawnow
    
    good = good & (H3 <= HL);
end
    if ~good
        keyboard
    end


close
end
return