g = 1 / (s*sqrt(2*pi))*exp(-((x-x_0).^2) / (2*s^2));
h = 1 / (pi*s) * sech(-(x-x_0) ./ (s));