function log_p = parzen_gaussian(A,B,h)

n=size(A,2);
h_n=h/sqrt(n);

log_p=[];
for i=1:size(h,2)
    log_p= [log_p log(sum(normpdf((A-B)/h_n(i)))/(n*h_n(i)))];
end
