function hessf=hessf0(x)
hessf=1-((exp(x)-exp(-x))/(exp(x)+exp(-x)))^2;