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];
xx=ones(10,7);
for i=1:7
     xx(:,i)=unifrnd(a(i),b(i),10,1);
end
zzz=3000000;
for i=1:10
    options=optimset('GradConstr','on');
    [xx1,zz]=fmincon(@Czfee,xx(i,:),[],[],[],[],a,b);
    if zzz>zz
        zzz=zz;
        xxx=xx1;
    end
end