function x=gauss_fn(t), x=exp(-0.5*t.^2);
% gauss_fn - Gaussian kernel function, exp(-t^2/2)
%
% gauss_fn :: [D] ~'any size array' -> [D->noneg]~'same size result'.
