max:
+5739.25  X1  +12.05  X2 +142.95  X3  +35.86  X4   +66.15  X5
   +9.49  X6 +105.38  X7 +185.09  X8  +64.59  X9  +327.09 X10
   +9.12 X11  +55.60 X12  +85.47 X13  +21.28 X14   +64.30 X15
+7416.42 X16  +38.54 X17  +53.04 X18 +121.11 X19 +3160.72 X20
 +336.60 X21  +57.59 X22  +25.67 X23 +504.66 X24   +63.94 X25
  +46.48 X26   +1.99 X27  +24.84 X28 +443.36 X29  +589.59 X30;
/* Número máximo de indivíduos por espécie */
 X1 - 100  Y1 <=0;
 X2 - 225  Y2 <=0;
 X3 - 225  Y3 <=0;
 X4 - 225  Y4 <=0;
 X5 - 225  Y5 <=0;
 X6 - 225  Y6 <=0;
 X7 - 225  Y7 <=0;
 X8 - 225  Y8 <=0;
 X9 - 225  Y9 <=0;
X10 - 225 Y10 <=0;
X11 - 225 Y11 <=0;
X12 - 225 Y12 <=0;
X13 - 350 Y13 <=0;
X14 - 225 Y14 <=0;
X15 - 225 Y15 <=0;
X16 - 225 Y16 <=0;
X17 - 225 Y17 <=0;
X18 - 225 Y18 <=0;
X19 - 225 Y19 <=0;
X20 - 185 Y20 <=0;
X21 - 350 Y21 <=0;
X22 - 225 Y22 <=0;
X23 - 225 Y23 <=0;
X24 - 350 Y24 <=0;
X25 - 225 Y25 <=0;
X26 - 225 Y26 <=0;
X27 - 225 Y27 <=0;
X28 - 225 Y28 <=0;
X29 - 225 Y29 <=0;
X30 - 225 Y30 <=0;
/* Contábeis para número de indivíduos por agrupamento */
MsNsPsZs-X14-X15-X23-X27=0;
MsNsPsZn-X25=0;
MsNsPnZs-X11-X17-X28-X29=0;
MsNsPnZn-X2-X3-X4-X5-X6-X7-X8-X12-X18-X19-X26=0;
MsNnPsZs=0;
MsNnPsZn=0;
MsNnPnZs=0;
MsNnPnZn=0;
MnNsPsZs-X13=0;
MnNsPsZn=0;
MnNsPnZs-X9-X10-X16-X21-X24-X30=0;
MnNsPnZn=0;
MnNnPsZs-X22=0;
MnNnPsZn=0;
MnNnPnZs-X1=0;
MnNnPnZn-X20=0;
/* Contábil para número total de indivíduos */
Xtotal-MsNsPsZs-MsNsPsZn
-MsNsPnZs-MsNsPnZn-MsNnPsZs-MsNnPsZn-MsNnPnZs-MsNnPnZn
-MnNsPsZs-MnNsPsZn-MnNsPnZs-MnNsPnZn-MnNnPsZs-MnNnPsZn
-MnNnPnZs-MnNnPnZn = 0;
/* Variáveis contábeis: número de espécies por categoria */
Ms-Y2-Y3-Y4-Y5-Y6-Y7-Y8-Y11-Y12-Y14-Y15-Y17-Y18-Y19-Y23-Y25-Y26-Y27-Y28-Y29 = 0;
Mn-Y1-Y9-Y10-Y13-Y16-Y20-Y21-Y22-Y24-Y30 = 0;
Ns-Y2-Y3-Y4-Y5-Y6-Y7-Y8-Y9-Y10-Y11-Y12-Y13-Y14-Y15-Y16-Y17-Y18-Y19-Y21-Y23-Y24-Y25-Y26-Y27-Y28-Y29-Y30 = 0;
Nn-Y1-Y20-Y22 = 0;
Ps-Y13-Y14-Y15-Y22-Y23-Y25-Y27 = 0;
Pn-Y1-Y2-Y3-Y4-Y5-Y6-Y7-Y8-Y9-Y10-Y11-Y12-Y16-Y17-Y18-Y19-Y20-Y21-Y23-Y24-Y26-Y28-Y29-Y30 = 0;
Zs-Y1-Y9-Y10-Y11-Y13-Y14-Y15-Y16-Y17-Y21-Y22-Y23-Y24-Y27-Y28-Y29-Y30 = 0;
Zn-Y2-Y3-Y4-Y5-Y6-Y7-Y8-Y12-Y18-Y19-Y20-Y25-Y26 = 0;
/* Contábil para número total de espécies */
Ytotal-Ms-Mn = 0;
/* Ocupação */
+9 X1+9 X2+ 9X3+ 9X4+ 9X5+ 9X6+ 9X7+ 9X8+ 9X9+9X10
+9X11+9X12+3X13+9X14+9X15+9X16+9X17+9X18+9X19+9X20
+9X21+9X22+9X23+3X24+9X25+9X26+9X27+9X28+9X29+9X30 <= 10000;
/* Metade das árvores deve ser de espécies nativas */
+MsNsPsZs+MsNsPsZn+MsNsPnZs+MsNsPnZn
+MnNsPsZs+MnNsPsZn+MnNsPnZs+MnNsPnZn - 0.50 Xtotal >= 0;
/* Pelo menos 1/3 das árvores deve ser de espécies madeireiras */
+MsNsPsZs+MsNsPsZn+MsNsPnZs+MsNsPnZn
+MsNnPsZs+MsNnPsZn+MsNnPnZs+MsNnPnZn - 0.33 Xtotal >= 0;
/* Pelo menos 1/3 das árvores deve ser de espécies zoocóricas */
+MsNsPsZs+MsNsPnZs+MsNnPsZs+MsNnPnZs
+MnNsPsZs+MnNsPnZs+MnNnPsZs+MnNnPnZs - 0.33 Xtotal >= 0;
/* Pelo menos 2 espécies não-pioneiras */
Pn >= 2;
/* Pelo menos 1 espécies pioneira */
Ps >= 1;
/* No mínimo 10 espécies */
/* Ytotal >= 10; */
/* Volume de madeira produzido no ano 30 */
+0.83 X2 +1.01 X3 +0.70 X4 +0.86 X5
+0.06 X6 +2.81 X7 +0.64 X8 +0.52X11
+0.62X12 +0.50X14 +1.16X15 +3.12X17
+0.85X18 +2.27X19 +0.62X23 +0.90X25
+0.55X26 +0.02X27 +0.32X28 +4.34X29 - VolMad30 = 0;
/* Declaração das variáveis binárias */
bin
 Y1, Y2, Y3, Y4, Y5, Y6, Y7, Y8, Y9,Y10,
Y11,Y12,Y13,Y14,Y15,Y16,Y17,Y18,Y19,Y20,
Y21,Y22,Y23,Y24,Y25,Y26,Y27,Y28,Y29,Y30;