function y = f(x, a, b)

y = exp(-1/2 * ((x - a) / b) .^ 2);