function res=F(x,theta)
res=x.*theta;
res = res / sum(res);