function y=fun_prim(x)
y=(1./cosh(x)).^2; %fun(x);
%y=yy*(1-yy);
%y=(1-yy).*(1+yy);