function F=fn(u)
F=(exp(u)-1)./(exp(1)-1);
%crude=mean(F)
%crudev=var(F)

