function y=gauss(erg,w)
  y = exp(-(erg/w).^2);
end
