function [y]=bs_estimate_pdf2_condx(Pxy, X, Y, x)

pp=1;
if pp==2
yy=Y(:,1);
xx=ones(length(Y),1)*x;

Py_condx=griddata(X,Y,Pxy,xx,yy,'cubic');
in=find( Py_condx == max(Py_condx));
y=yy(in);
return
end

for i=1:2
   Pxy=utl_shapiro2(Pxy,2,2);
end
xtmp=X(1,:);
dx=abs(xtmp-x);
in =find ( dx(:) == min(dx(:)));
imax = find( Pxy(:,in) == max(Pxy(:,in)) );
imax=imax(1);
y=Y(imax,in);