function [o,v,b,t1,V1]=optimal
u=rand(1,10000);
t1=cputime;
Y1=(0.53/2)*(callopt(0.47+0.53*u)+callopt(1-0.53*u));
t1=[t1 cputime];
Y2=0.37/2*(callopt(0.47+0.37*u)+callopt(0.84-0.37*u))+0.16/2*(callopt(0.84+0.16*u)+callopt(1-0.16*u));
t1=[t1 cputime];
Y3=0.37*callopt(0.47+0.37*u)+0.16*callopt(1-0.16*u);
t1=[t1 cputime];
Y4=control(u);
t1=[t1 cputime];
Y5=importance(u);
t1=[t1 cputime];
X=[Y1',Y2',Y3',Y4',Y5'];
%mean(X);
V=cov(X)
Z=ones(5,1);
V1=inv(V);
b=V1*Z/(Z'*V1*Z)
o=mean(X*b)
t1=[t1 cputime];
v=1/(Z'*V1*Z)
t1=diff(t1)
Crude=callopt(u);
Efficiency=var(u)/v
