function f=rexas(e)
k1=sqrt(2*e); k2=sqrt(2*(e-1.0)); 
alfa=5;
f=((sinh(pi*(k2-k1)/alfa))/(sinh(pi*(k2+k1)/alfa)))^2;