function w = gaussian_weight_1d(N,par); 

w = zeros(1,N);
w(1,:) = par(1);

for i = 1:N
   xv = i-1;
   if (xv>(N/2))  
      xv = N+1-i;
   end
   di2 = xv^2;
   eval = par(2) * exp( -(di2/(2.*(0.88^2))) );
   ival = 0.2 * par(2) * exp( -(di2/(2.*(2.03^2))) );
   w(1,i) = w(1,i)+eval-ival;
end