function y = pes(erg,eb,w)
  y = 0;
  for k=1:length(eb)
    if(eb(k)<0.0)
      y = y + gauss(erg-eb(k),w);
    end
  end
end

