function y = gauss(x,std)
% This function computes the probability density function 
% at x for a normal distribution with zero expectation and 
% standard deviation std

y = exp(-x^2/(2*std^2)) / (std*sqrt(2*pi));
