Ubern [SECS2D] [secs1D]
[bp,bn]=Ubern(x)

Bernoulli function
bp = B(x)=x/(exp(x)-1)
bn = B(-x)=x+B(x)