function [fv,xhat,K] = boundYxhat(P,phi,Xt,m,n)
yy = [];
KK = {};
for xhat = 1:m
  [yy(xhat),K] = condexphat(P,phi,Xt,m,n,xhat);
  KK{xhat} = K;
end
[fv,xhat] = max(abs(yy));
K = KK{xhat};
return