a=[0.075,0.225,0.075,0.075,1.125,12,0.5625];
b=[0.125,0.375,0.125,0.125,1.875,20,0.935];
x=ones(10,7);
for i=1:7
     x(:,i)=unifrnd(a(i),b(i),10,1);
end
D=[0.05,0.05,0.05,0.1,0.1,0.05,0.05]/3;
options=optimset('GradConstr','on');
[xx,zz]=fmincon(@Czfee,x(i,:),[],[],[],[],a,b);