ParameterDomain = [0,];
Constant Q=[1600,1500,1600,1700,2000,1600],
         C=[100,130,120,210,240,640];
ConstStr b1=(x1+x2+x3)/q1, b2=(x1+x3)/q2,
         b3=(x1+x6+x3)/q3, b4=(x4+x5+x6)/q4,
         b5=(x2+x4+x5+x6)/q5,b6=(x4+x5+x2)/q6,
         D=x1*c1+x2*c2+x3*c3+x4*c4+x5*c5+x6*c6;
ConstStr aa=Sum(i=1:6)((x[i]*c[i]/D)*ln(x[i]*c[i]/D)),
         bb=Sum(i=1:6)(b[i]*ln(b[i]));
MinFunction 0.7*bb+0.3*aa;
            x1+x2+x3+x4+x5+x6-2400=0;