	function y=rexaeck(a,E)
    b=8*a;
    k=sqrt(2*E);
    beta=sqrt(k^2-a);
    delta=sqrt(b-1/4);
	y = (cosh(2*pi*(k-beta))+cosh(2*pi*delta))/(cosh(2*pi*(k+beta))+cosh(2*pi*delta));
	end