function res = gauss(x,mean,stddev)
% find the value of a Gauss function  at argument x  
% with mean value MEAN and standard deviation STDDEV
  res = exp(-1/2*(x - mean).^2/stddev^2)*1/(sqrt(2*pi)*stddev);
