function value=esatta(e,x)
c1=(exp(-1/sqrt(e))-1)/(exp(1/sqrt(e))-exp(-1/sqrt(e)));
c0=-1-c1;
value=1+c0*exp(-1/sqrt(e).*x)+c1*exp(1/sqrt(e).*x);
