function R = gaussian_1d(sigma, N); a = linspace (-N/2, N/2, N); R = zeros(size(a)); for i = 1:size(a,2) R(i) = gaussian(0, sigma, a(i)); end