function F=control(u)
%u=unifrnd(0,1,1,n);
FN=callopt(u);
u=max(0,u-0.47);
CN=6*u.^2+u;
intg=0.53^2/2+0.53^3*2;
F=FN-CN+intg;
CVest=mean(F)
CVvar=var(FN-CN)
%CVest =
%    0.4614
%CVvar =
%  0.014297