function y=HeLe(x)

y = log( 1+abs(x).^2 )./(log( 1+abs(1).^2 ) );
