function [x,y] = gaussian(mu,sigma)

% gaussian(mu,sigma) models the Gaussian probability density
% see, Bevington and Robinson, 2nd Edition, Section 2.3
% mu is the distribution mean
% gamma is the distribution standard deviation

x=[mu-5*sigma:.1:mu+5*sigma];
y=1/(sigma*sqrt(2*pi))*exp(-1/2*((x-mu)./sigma).^2);
plot(x,y);
