function y=fun_eig(x)
y=cos(x)+1/cosh(x);

