function w = dgaussian_weight_1d(N,par); 

w = zeros(1,N);

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