b=Ubernoulli(x,sg) Bernoulli function b = B(x)=x/(exp(x)-1) if sg==1 b = B(-x)=x+B(x) if sg==0 also works if x is a vector