function [bdvalue] =  epsilonbound(epsilon)

global g2

g =  g2;
epsilong = epsilon + g;

bdvalue = (1./epsilon.^2).*(   (1-(epsilong)).*log((1-(epsilong))./(1-g))  + (epsilong).*log((epsilong)/g)) ;