function y = bigu(k)

   y = exp(-6*(k-2).^2)*sqrt(24*pi);