颜色= optimvar (“颜色”, (“黑色”，“白色”，“红色”，“绿色”，“蓝色”]、[“青色”，“红色”，“黄色”，“黑色”]);

[idxadd, idxsub] = findindex(颜色,{“红色”，“黑”}, {“黑”}）

颜色= optimvar (“颜色”, (“黑色”，“白色”，“红色”，“绿色”，“蓝色”]、[“青色”，“红色”，“黄色”，“黑色”]);

idx = findindex(颜色,(“白色”，“红色”，“绿色”，“蓝色”]、[“黑色”，“青色”，“红色”，“黄色”]）

rng (0)%的再现性p = optimproblem (“ObjectiveSense”，“最大化”）;流= optimvar (“流”，．..｛“苹果”，“橘子”，“香蕉”，“浆果”}， {“纽约”，“bo”，“宽松”}，．..下界的，0，“类型”，“整数”）;p.Objective =总和(总和(兰德(4,3)。*流));p.Constraints.NYC =兰德(1、4)*流(:,“纽约”) < = 10;p.Constraints.BOS =兰德(1、4)*流(:,“bo”) < = 12;p.Constraints.LAX =兰德(1、4)*流(:,“宽松”) < = 35;索尔=解决(p);

[idxFruit,idxAirports] = finddex (flow， {“橘子”，“浆果”}， {“纽约”，“宽松”}）

橙莓= sol.flow(idxFruit, idxAirports)

Idx = finddex (flow， {“浆果”，“苹果”，“橘子”}， {“纽约”，“bo”，“宽松”}）

optimalFlow = sol.flow (idx)

土地= [“irr-good”，“irr-poor”，“dry-good”，“dry-poor”];作物= [“wheat-lentil”，“小麦加工”，“barley-chickpea”，“barley-lentil”，“wheat-onion”，“barley-onion”];犁= [“传统”，“机械化”];xcrop = optimvar (“xcrop”、土地、庄稼、犁,下界的, 0);

x0。xcrop = 0(大小(xcrop));

[idxLand, idxCrop, idxPlough] = finddex (xcrop， [“dry-good”，“dry-poor”]，．..［“wheat-onion”，“wheat-lentil”]，“传统”）;x0.xcrop (idxLand idxCrop idxPlough) = 3000;

idx = findindex (xcrop,．..［“dry-good”，“irr-poor”，“irr-good”]，．..［“小麦加工”，“barley-onion”，“barley-chickpea”]，．..［“机械化”，“传统”，“机械化”]);x0.xcrop (idx) =(2000、5000、3500);