function [x]=debye_gamma(ess,eh,fr,gamma,f)

x = eh + (ess - eh)./(1 + i* f./fr) - i*f*1e-9*gamma;