function fx=f0(x)
fx=log(exp(x)+exp(-x));