function k=kernel(x)
%x must be a row vector.
k=exp(-x.^2/2)/sqrt(2*pi);