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