function s = kernel(u)

%gaussian kernel

s=exp((-0.5)*(u.^2))/sqrt(2*pi);