$offlisting * * Equation counts * Total E G L N X C B * 2 1 0 1 0 0 0 0 * * Variable counts * x b i s1s s2s sc si * Total cont binary integer sos1 sos2 scont sint * 41 21 0 20 0 0 0 0 * FX 0 * * Nonzero counts * Total const NL DLL * 81 41 40 0 * * Solve m using MINLP minimizing objvar; Variables i1,i2,i3,i4,i5,i6,i7,i8,i9,i10,i11,i12,i13,i14,i15,i16,i17,i18,i19 ,i20,x21,x22,x23,x24,x25,x26,x27,x28,x29,x30,x31,x32,x33,x34,x35,x36 ,x37,x38,x39,x40,objvar; Integer Variables i1,i2,i3,i4,i5,i6,i7,i8,i9,i10,i11,i12,i13,i14,i15,i16,i17 ,i18,i19,i20; Equations e1,e2; e1.. 1.1*i1 + 1.86*i2 + 0.86*i3 + 1.23*i4 + 0.72*i5 + 0.47*i6 + 0.98*i7 + 0.39*i8 + 0.25*i9 + 0.41*i10 + 0.29*i11 + 0.38*i12 + 1.49*i13 + 1.27*i14 + 0.56*i15 + 1.08*i16 + 1.39*i17 + i18 + 1.37*i19 + 0.89*i20 + 0.25*x21 + 0.98*x22 + 1.71*x23 + 1.75*x24 + 0.54*x25 + 0.42*x26 + 1.13*x27 + 1.28*x28 + 0.83*x29 + 0.41*x30 + 1.9*x31 + 0.16*x32 + 0.21*x33 + 0.28*x34 + 0.33*x35 + 1.24*x36 + 1.15*x37 + 0.1*x38 + 1.86*x39 + 1.46*x40 - objvar =E= 0; e2.. -0.2*i1**0.035*i2**0.003*i3**0.04*i4**0.044*i5**0.046*i6**0.04*i7**0.037* i8**0.024*i9**0.008*i10**0.019*i11**0.006*i12**0.001*i13**0.044*i14**0.014 *i15**0.014*i16**0.016*i17**0.022*i18**0.03*i19**0.001*i20**0.039*x21** 0.026*x22**0.04*x23**0.016*x24**0.021*x25**0.003*x26**0.008*x27**0.031*x28 **0.015*x29**0.042*x30**0.006*x31**0.046*x32**0.025*x33**0.033*x34**0.047* x35**0.013*x36**0.019*x37**0.022*x38**0.036*x39**0.038*x40**0.005 =L= -1; * set non-default bounds i1.lo = 1; i1.up = 10; i2.lo = 1; i2.up = 10; i3.lo = 1; i3.up = 10; i4.lo = 1; i4.up = 10; i5.lo = 1; i5.up = 10; i6.lo = 1; i6.up = 10; i7.lo = 1; i7.up = 10; i8.lo = 1; i8.up = 10; i9.lo = 1; i9.up = 10; i10.lo = 1; i10.up = 10; i11.lo = 1; i11.up = 10; i12.lo = 1; i12.up = 10; i13.lo = 1; i13.up = 10; i14.lo = 1; i14.up = 10; i15.lo = 1; i15.up = 10; i16.lo = 1; i16.up = 10; i17.lo = 1; i17.up = 10; i18.lo = 1; i18.up = 10; i19.lo = 1; i19.up = 10; i20.lo = 1; i20.up = 10; x21.lo = 1E-5; x21.up = 10; x22.lo = 1E-5; x22.up = 10; x23.lo = 1E-5; x23.up = 10; x24.lo = 1E-5; x24.up = 10; x25.lo = 1E-5; x25.up = 10; x26.lo = 1E-5; x26.up = 10; x27.lo = 1E-5; x27.up = 10; x28.lo = 1E-5; x28.up = 10; x29.lo = 1E-5; x29.up = 10; x30.lo = 1E-5; x30.up = 10; x31.lo = 1E-5; x31.up = 10; x32.lo = 1E-5; x32.up = 10; x33.lo = 1E-5; x33.up = 10; x34.lo = 1E-5; x34.up = 10; x35.lo = 1E-5; x35.up = 10; x36.lo = 1E-5; x36.up = 10; x37.lo = 1E-5; x37.up = 10; x38.lo = 1E-5; x38.up = 10; x39.lo = 1E-5; x39.up = 10; x40.lo = 1E-5; x40.up = 10; Model m / all /; m.limrow=0; m.limcol=0; m.tolproj=0.0; $if NOT '%gams.u1%' == '' $include '%gams.u1%' $if not set MINLP $set MINLP MINLP Solve m using %MINLP% minimizing objvar;