function checketakap(m,n)
P = pnormdim(rand(m^n,1),1);
F = 0*P;
for i=1:n
  xx = randxx(m,n); y=xx{1};
  xx = randxx(m,1); xhat=xx{1};
  [Yhat,K,TV] = condexphat(P,F,y(1:i),m,n,xhat);
  for j=i+1:n
    h1 = sum(pl(K));
    h2 = etaijX(P,m,n,i,j,y(1:i));
    er = abs(h1-h2)
    if er > 1e-6
      keyboard
    end
    K = getK1(K,m,n-j+1);
        
  end
end
  

return



function [hh,Ht,H] = gethhnX(P,m,n,X)
hh = zeros(n);
for i=1:n
  for j=i+1:n
    h = etaijX(P,m,n,i,j,X);
    hh(i,j) = h;
  end
end
Ht = 1+sum(hh,2)';
H = max(Ht);
return