function w = gaussian_1d(N,x,xsig); 

w = zeros(1,N);

for i = 1:N
   xv = i-x;
   if (xv>(N/2))  
      xv = i-N;
   end
   if (xv<(-N/2))
      xv = i+N;
   end
   di2 = xv^2;
   w(1,i) = exp( -(di2/(2.*(xsig^2))) );
 end