function x=epsc(ess,einf,fr,f,a)

x=einf+(ess-einf)./((1+(i*f./fr)).^(1-a))