function [e]=doubled(ess,einf,eh,fr1,fr2,f)

e=einf+(ess-eh)./(1+(i*f./fr1))+(eh-einf)./((1+(i*f./fr2)))