Задаване на променливи

[go6aih]

Как да направя sx или va променливи защото при този вид отговора на системата е -
NMaximize::bcons: The following constraints are not valid: {30. va[1]+39. va[2]+46. va[3]+20. va[4]<={256.}}. Constraints should be equalities, inequalities, or domain specifications involving the variables. >>

ClearAll;
n = 5;
m = 4;
c = N[RandomInteger[50, {n, m}]];
z = N[RandomInteger[500, {n, 1}]];
q = N[RandomInteger[100, {m, 1}]];
Table[c];
sx = Table[va[j], {j, 1, m}]
Do[ogr = Table[sx[[j]]*c[[All, j]], {j, 1, m}], {i, 1, n}];
s1p = Sum[ogr[[j]], {j, 1, m}];
qq = Table[q[[j]]*sx[[j]], {j, 1, m}];
zz = Table[s1p[[i, All]] <= z[[i]], {i, 1, n}];
sqq = Sum[qq[[j]], {j, 1, m}];
Maximize[sqq[[1]], zz[[1]], sx[[1]]]

Ако последният ред е
Maximize[sqq, zz, {sx[[j]], {j, 1, m}}]
връща - Maximize::vdom: Variable domain {{va[1],va[2],va[3],va[4]}[[j]],{j,1,4}} should be either Reals or Integers. >>

Някой да помогне...

[go6aih]

Това е решението:

"7 зад. "
"n - вида продукция"
"m - вида суровини "
ClearAll;
m = 5;
n = 4;
c = RandomInteger[10, {m, n}];
cc = c // MatrixForm
" ЗАПАСИ "
z = RandomInteger[{500, 1000}, {m, 1}]
" ПРОИЗВОДСТВЕНИ ЦЕНИ "
q = RandomInteger[{500, 1000}, {n, 1}]
Table[c];
tabx = Table[x[j], {j, 1, n}];
Do[sxc = Table[c[[All, j]]*tabx[[j]], {j, 1, n}], {i, 1, m}];
srow = Sum[sxc[[j]], {j, 1, n}];
qx = Table[q[[j]]*tabx[[j]], {j, 1, n}];
Table[qx];
conz = Table[srow[[i, All]] <= z[[i, 1]], {i, 1, m}];
suq = Sum[qx[[j]], {j, 1, n}];
w = Maximize[{suq[[1]], {conz,
tabx[[1]] >= 0 && tabx[[2]] >= 0 && tabx[[3]] >= 0 &&
tabx[[4]] >= 0}}, tabx]

Ако му трябва на някого....