function [f] = gaussian(mu, sigma, x)
% Gaussian.m
% (c) 2007, University of Surrey, UK
%
%  Usage:
% function [f] = gaussian(mu, sigma, x)
% 
%  Input:  mu is the mean of the distribution and sigma is standard deviation,
%          x are the times at which the function is to be evaluated
%  Output: f is the function output, according to 
%
%  f(x) = (1 / [sqrt(2 * pi) * sigma]) * exp(-0.5 * ((x - mu) / sigma)^2).

x = x(:);
f = exp(-0.5 * ((x - mu) / sigma).^2) / (sqrt(2 * pi) * sigma);
