function DoG = differenceOfGaussians1D(a, x)

DoG = a(1) * exp(-((x - a(2)) / a(3)).^2) + ...
      a(4) * exp(-((x - a(2)) / a(5)).^2);
  
end