watersym2.gms

	[ MINLP World Home \| Board \| Solvers \| MINLPLib \| Gams World Forum \| Search \| Contact ]

watersym2.gms:

References:

GAMS Development, Reformulation of Models from GAMS Library.
Brooke, A, Drud, A S, and Meeraus, A, Modeling Systems and Nonlinear Programming in a Research Environment. In Ragavan, R, and Rohde, S M, Eds, Computers in Engineering, Vol. III. ACME, 1985.
Drud, A S, and Rosenborg, A, Dimensioning Water Distribution Networks. Masters thesis, Institute of Mathematical Statistics and Operations Research, Technical University of Denmark, 1973. (in Danish)
Original source: GAMS model waterx.gms from Gams Model Library

Points: p1 p2
Best known point (p2): Solution value 940.92 (global optimum?, SCIP certificate)

$offlisting * MINLP written by GAMS Convert at 04/17/01 16:42:12 * * Equation counts * Total E G L N X * 202 96 36 70 0 0 * * Variable counts * x b i s1s s2s sc si * Total cont binary integer sos1 sos2 scont sint * 322 98 28 0 196 0 0 0 * FX 0 0 0 0 0 0 0 0 * * Nonzero counts * Total const NL DLL * 1247 1187 60 0 * * Solve m using MINLP minimizing objvar; Sets s1 /127*133/,s2 /134*140/,s3 /141*147/,s4 /148*154/,s5 /155*161/ ,s6 /162*168/,s7 /169*175/,s8 /176*182/,s9 /183*189/,s10 /190*196/ ,s11 /197*203/,s12 /204*210/,s13 /211*217/,s14 /218*224/,s15 /225*231/ ,s16 /232*238/,s17 /239*245/,s18 /246*252/,s19 /253*259/,s20 /260*266/ ,s21 /267*273/,s22 /274*280/,s23 /281*287/,s24 /288*294/,s25 /295*301/ ,s26 /302*308/,s27 /309*315/,s28 /316*322/; Variables x1,x2,x3,x4,x5,x6,x7,x8,x9,x10,x11,x12,x13,x14,x15,x16,x17,x18,x19 ,x20,x21,x22,x23,x24,x25,x26,x27,x28,x29,x30,x31,x32,x33,x34,x35,x36 ,x37,x38,x39,x40,x41,x42,x43,x44,x45,x46,x47,x48,x49,x50,x51,x52,x53 ,x54,x55,x56,x57,x58,x59,x60,x61,x62,x63,x64,x65,x66,x67,x68,x69,x70 ,x71,x72,x73,x74,x75,x76,x77,x78,x79,x80,x81,x82,x83,x84,x85,x86,x87 ,x88,x89,x90,x91,x92,x93,x94,x95,x96,x97,objvar,b99,b100,b101,b102 ,b103,b104,b105,b106,b107,b108,b109,b110,b111,b112,b113,b114,b115 ,b116,b117,b118,b119,b120,b121,b122,b123,b124,b125,b126,s1s1(s1) ,s1s2(s2),s1s3(s3),s1s4(s4),s1s5(s5),s1s6(s6),s1s7(s7),s1s8(s8) ,s1s9(s9),s1s10(s10),s1s11(s11),s1s12(s12),s1s13(s13),s1s14(s14) ,s1s15(s15),s1s16(s16),s1s17(s17),s1s18(s18),s1s19(s19),s1s20(s20) ,s1s21(s21),s1s22(s22),s1s23(s23),s1s24(s24),s1s25(s25),s1s26(s26) ,s1s27(s27),s1s28(s28); Positive Variables x1,x2,x3,x4,x5,x6,x7,x8,x9,x10,x11,x12,x13,x14,x15,x16,x17 ,x18,x19,x20,x21,x22,x23,x24,x25,x26,x27,x28,x93,x94; Binary Variables b99,b100,b101,b102,b103,b104,b105,b106,b107,b108,b109,b110 ,b111,b112,b113,b114,b115,b116,b117,b118,b119,b120,b121,b122,b123 ,b124,b125,b126; SOS1 Variables s1s1(s1),s1s2(s2),s1s3(s3),s1s4(s4),s1s5(s5),s1s6(s6),s1s7(s7) ,s1s8(s8),s1s9(s9),s1s10(s10),s1s11(s11),s1s12(s12),s1s13(s13) ,s1s14(s14),s1s15(s15),s1s16(s16),s1s17(s17),s1s18(s18),s1s19(s19) ,s1s20(s20),s1s21(s21),s1s22(s22),s1s23(s23),s1s24(s24),s1s25(s25) ,s1s26(s26),s1s27(s27),s1s28(s28); Equations e1,e2,e3,e4,e5,e6,e7,e8,e9,e10,e11,e12,e13,e14,e15,e16,e17,e18,e19 ,e20,e21,e22,e23,e24,e25,e26,e27,e28,e29,e30,e31,e32,e33,e34,e35,e36 ,e37,e38,e39,e40,e41,e42,e43,e44,e45,e46,e47,e48,e49,e50,e51,e52,e53 ,e54,e55,e56,e57,e58,e59,e60,e61,e62,e63,e64,e65,e66,e67,e68,e69,e70 ,e71,e72,e73,e74,e75,e76,e77,e78,e79,e80,e81,e82,e83,e84,e85,e86,e87 ,e88,e89,e90,e91,e92,e93,e94,e95,e96,e97,e98,e99,e100,e101,e102,e103 ,e104,e105,e106,e107,e108,e109,e110,e111,e112,e113,e114,e115,e116 ,e117,e118,e119,e120,e121,e122,e123,e124,e125,e126,e127,e128,e129 ,e130,e131,e132,e133,e134,e135,e136,e137,e138,e139,e140,e141,e142 ,e143,e144,e145,e146,e147,e148,e149,e150,e151,e152,e153,e154,e155 ,e156,e157,e158,e159,e160,e161,e162,e163,e164,e165,e166,e167,e168 ,e169,e170,e171,e172,e173,e174,e175,e176,e177,e178,e179,e180,e181 ,e182,e183,e184,e185,e186,e187,e188,e189,e190,e191,e192,e193,e194 ,e195,e196,e197,e198,e199,e200,e201,e202; e1.. - x1 - x2 - x3 + x8 + x14 + x26 + x93 =E= 0; e2.. - x4 - x5 - x6 - x7 + x9 + x20 + x24 + x27 + x94 =E= 0; e3.. x1 + x4 - x8 - x9 - x10 - x11 - x12 - x13 + x15 + x17 + x21 + x28 =E= 1.212; e4.. x2 + x10 - x14 - x15 - x16 + x18 =E= 0.452; e5.. x11 + x16 - x17 - x18 - x19 + x22 =E= 0.245; e6.. x5 + x12 + x19 - x20 - x21 - x22 - x23 + x25 =E= 0.652; e7.. x6 + x23 - x24 - x25 =E= 0.252; e8.. x3 + x7 + x13 - x26 - x27 - x28 =E= 0.456; e9.. x29 - 38721.1970117411*s1s1('127') - 2543.8701482414*s1s1('128') - 207.747320703761*s1s1('129') - 23.9314504121258*s1s1('130') - 1.5722267648148*s1s1('131') - 0.181112645550961*s1s1('132') - 0.0390863672545667*s1s1('133') =E= 0; e10.. x30 - 32510.4890865135*s1s2('134') - 2135.84468132099*s1s2('135') - 174.425573683688*s1s2('136') - 20.0929521164322*s1s2('137') - 1.32004857865156*s1s2('138') - 0.152062982061963*s1s2('139') - 0.0328170876451919*s1s2('140') =E= 0; e11.. x31 - 63468.4628982673*s1s3('141') - 4169.69361956223*s1s3('142') - 340.521578201805*s1s3('143') - 39.2263796008983*s1s3('144') - 2.57705917665854*s1s3('145') - 0.296864304610023*s1s3('146') - 0.0640670186196026*s1s3('147') =E= 0; e12.. x32 - 50797.5773435889*s1s4('148') - 3337.25325093014*s1s4('149') - 272.539627020641*s1s4('150') - 31.3951994533022*s1s4('151') - 2.06257339263358*s1s4('152') - 0.237598120158509*s1s4('153') - 0.0512766370081929*s1s4('154') =E= 0; e13.. x33 - 59165.7349698592*s1s5('155') - 3887.01689524085*s1s5('156') - 317.436542928413*s1s5('157') - 36.5670992066393*s1s5('158') - 2.40235218067626*s1s5('159') - 0.27673893405488*s1s5('160') - 0.0597237127048799*s1s5('161') =E= 0; e14.. x34 - 32977.2294678044*s1s6('162') - 2166.50816836621*s1s6('163') - 176.929733450444*s1s6('164') - 20.3814187742893*s1s6('165') - 1.339*s1s6('166') - 0.154246090843839*s1s6('167') - 0.0332882297421199*s1s6('168') =E= 0; e15.. x35 - 33843.9321019273*s1s7('169') - 2223.4480134252*s1s7('170') - 181.579774357788*s1s7('171') - 20.9170801874496*s1s7('172') - 1.37419139860501*s1s7('173') - 0.158299963634093*s1s7('174') - 0.0341631060391402*s1s7('175') =E= 0; e16.. x36 - 38721.1970117411*s1s8('176') - 2543.8701482414*s1s8('177') - 207.747320703761*s1s8('178') - 23.9314504121258*s1s8('179') - 1.5722267648148*s1s8('180') - 0.181112645550961*s1s8('181') - 0.0390863672545667*s1s8('182') =E= 0; e17.. x37 - 50797.5773435889*s1s9('183') - 3337.25325093014*s1s9('184') - 272.539627020641*s1s9('185') - 31.3951994533022*s1s9('186') - 2.06257339263358*s1s9('187') - 0.237598120158509*s1s9('188') - 0.0512766370081929*s1s9('189') =E= 0; e18.. x38 - 31810.181054648*s1s10('190') - 2089.8364782095*s1s10('191') - 170.668274619734*s1s10('192') - 19.660130090483*s1s10('193') - 1.2916134290104*s1s10('194') - 0.148787395299671*s1s10('195') - 0.0321101751776739*s1s10('196') =E= 0; e19.. x39 - 39461.9459070343*s1s11('197') - 2592.53519858857*s1s11('198') - 211.721593458417*s1s11('199') - 24.3892667200816*s1s11('200') - 1.60230396616872*s1s11('201') - 0.184577388442944*s1s11('202') - 0.0398341019735132*s1s11('203') =E= 0; e20.. x40 - 32977.2294678044*s1s12('204') - 2166.50816836621*s1s12('205') - 176.929733450444*s1s12('206') - 20.3814187742893*s1s12('207') - 1.339*s1s12('208') - 0.154246090843839*s1s12('209') - 0.0332882297421199*s1s12('210') =E= 0; e21.. x41 - 52785.5148814787*s1s13('211') - 3467.85497167945*s1s13('212') - 283.205327698691*s1s13('213') - 32.6238347301504*s1s13('214') - 2.14329116080854*s1s13('215') - 0.246896402610059*s1s13('216') - 0.0532833223041444*s1s13('217') =E= 0; e22.. x42 - 32510.4890865135*s1s14('218') - 2135.84468132099*s1s14('219') - 174.425573683688*s1s14('220') - 20.0929521164322*s1s14('221') - 1.32004857865156*s1s14('222') - 0.152062982061963*s1s14('223') - 0.0328170876451919*s1s14('224') =E= 0; e23.. x43 - 31810.181054648*s1s15('225') - 2089.8364782095*s1s15('226') - 170.668274619734*s1s15('227') - 19.660130090483*s1s15('228') - 1.2916134290104*s1s15('229') - 0.148787395299671*s1s15('230') - 0.0321101751776739*s1s15('231') =E= 0; e24.. x44 - 30677.4142839491*s1s16('232') - 2015.41699236491*s1s16('233') - 164.590743970989*s1s16('234') - 18.9600290116536*s1s16('235') - 1.24561882211213*s1s16('236') - 0.143489047044288*s1s16('237') - 0.0309667255575633*s1s16('238') =E= 0; e25.. x45 - 39461.9459070343*s1s17('239') - 2592.53519858857*s1s17('240') - 211.721593458417*s1s17('241') - 24.3892667200816*s1s17('242') - 1.60230396616872*s1s17('243') - 0.184577388442944*s1s17('244') - 0.0398341019735132*s1s17('245') =E= 0; e26.. x46 - 30677.4142839491*s1s18('246') - 2015.41699236491*s1s18('247') - 164.590743970989*s1s18('248') - 18.9600290116536*s1s18('249') - 1.24561882211213*s1s18('250') - 0.143489047044288*s1s18('251') - 0.0309667255575633*s1s18('252') =E= 0; e27.. x47 - 28361.2795383154*s1s19('253') - 1863.25366856746*s1s19('254') - 152.164196629274*s1s19('255') - 17.5285530220005*s1s19('256') - 1.15157500841239*s1s19('257') - 0.132655670919396*s1s19('258') - 0.0286287479053886*s1s19('259') =E= 0; e28.. x48 - 59165.7349698592*s1s20('260') - 3887.01689524085*s1s20('261') - 317.436542928413*s1s20('262') - 36.5670992066393*s1s20('263') - 2.40235218067626*s1s20('264') - 0.27673893405488*s1s20('265') - 0.0597237127048799*s1s20('266') =E= 0; e29.. x49 - 32977.2294678044*s1s21('267') - 2166.50816836621*s1s21('268') - 176.929733450444*s1s21('269') - 20.3814187742893*s1s21('270') - 1.339*s1s21('271') - 0.154246090843839*s1s21('272') - 0.0332882297421199*s1s21('273') =E= 0; e30.. x50 - 28361.2795383154*s1s22('274') - 1863.25366856746*s1s22('275') - 152.164196629274*s1s22('276') - 17.5285530220005*s1s22('277') - 1.15157500841239*s1s22('278') - 0.132655670919396*s1s22('279') - 0.0286287479053886*s1s22('280') =E= 0; e31.. x51 - 50797.5773435889*s1s23('281') - 3337.25325093014*s1s23('282') - 272.539627020641*s1s23('283') - 31.3951994533022*s1s23('284') - 2.06257339263358*s1s23('285') - 0.237598120158509*s1s23('286') - 0.0512766370081929*s1s23('287') =E= 0; e32.. x52 - 32977.2294678044*s1s24('288') - 2166.50816836621*s1s24('289') - 176.929733450444*s1s24('290') - 20.3814187742893*s1s24('291') - 1.339*s1s24('292') - 0.154246090843839*s1s24('293') - 0.0332882297421199*s1s24('294') =E= 0; e33.. x53 - 50797.5773435889*s1s25('295') - 3337.25325093014*s1s25('296') - 272.539627020641*s1s25('297') - 31.3951994533022*s1s25('298') - 2.06257339263358*s1s25('299') - 0.237598120158509*s1s25('300') - 0.0512766370081929*s1s25('301') =E= 0; e34.. x54 - 63468.4628982673*s1s26('302') - 4169.69361956223*s1s26('303') - 340.521578201805*s1s26('304') - 39.2263796008983*s1s26('305') - 2.57705917665854*s1s26('306') - 0.296864304610023*s1s26('307') - 0.0640670186196026*s1s26('308') =E= 0; e35.. x55 - 33843.9321019273*s1s27('309') - 2223.4480134252*s1s27('310') - 181.579774357788*s1s27('311') - 20.9170801874496*s1s27('312') - 1.37419139860501*s1s27('313') - 0.158299963634093*s1s27('314') - 0.0341631060391402*s1s27('315') =E= 0; e36.. x56 - 52785.5148814787*s1s28('316') - 3467.85497167945*s1s28('317') - 283.205327698691*s1s28('318') - 32.6238347301504*s1s28('319') - 2.14329116080854*s1s28('320') - 0.246896402610059*s1s28('321') - 0.0532833223041444*s1s28('322') =E= 0; e37.. - sqr(x1)*x29 + x57 - x59 - x65 =E= 0; e38.. - sqr(x2)*x30 + x57 - x60 - x66 =E= 0; e39.. - sqr(x3)*x31 + x57 - x64 - x67 =E= 0; e40.. - sqr(x4)*x32 + x58 - x59 - x68 =E= 0; e41.. - sqr(x5)*x33 + x58 - x62 - x69 =E= 0; e42.. - sqr(x6)*x34 + x58 - x63 - x70 =E= 0; e43.. - sqr(x7)*x35 + x58 - x64 - x71 =E= 0; e44.. - sqr(x8)*x36 - x57 + x59 - x72 =E= 0; e45.. - sqr(x9)*x37 - x58 + x59 - x73 =E= 0; e46.. - sqr(x10)*x38 + x59 - x60 - x74 =E= 0; e47.. - sqr(x11)*x39 + x59 - x61 - x75 =E= 0; e48.. - sqr(x12)*x40 + x59 - x62 - x76 =E= 0; e49.. - sqr(x13)*x41 + x59 - x64 - x77 =E= 0; e50.. - sqr(x14)*x42 - x57 + x60 - x78 =E= 0; e51.. - sqr(x15)*x43 - x59 + x60 - x79 =E= 0; e52.. - sqr(x16)*x44 + x60 - x61 - x80 =E= 0; e53.. - sqr(x17)*x45 - x59 + x61 - x81 =E= 0; e54.. - sqr(x18)*x46 - x60 + x61 - x82 =E= 0; e55.. - sqr(x19)*x47 + x61 - x62 - x83 =E= 0; e56.. - sqr(x20)*x48 - x58 + x62 - x84 =E= 0; e57.. - sqr(x21)*x49 - x59 + x62 - x85 =E= 0; e58.. - sqr(x22)*x50 - x61 + x62 - x86 =E= 0; e59.. - sqr(x23)*x51 + x62 - x63 - x87 =E= 0; e60.. - sqr(x24)*x52 - x58 + x63 - x88 =E= 0; e61.. - sqr(x25)*x53 - x62 + x63 - x89 =E= 0; e62.. - sqr(x26)*x54 - x57 + x64 - x90 =E= 0; e63.. - sqr(x27)*x55 - x58 + x64 - x91 =E= 0; e64.. - sqr(x28)*x56 - x59 + x64 - x92 =E= 0; e65.. x65 + 12*b99 =L= 12; e66.. x66 + 12*b100 =L= 12; e67.. x67 + 12*b101 =L= 12; e68.. x68 + 12*b102 =L= 12; e69.. x69 + 12*b103 =L= 12; e70.. x70 + 12*b104 =L= 12; e71.. x71 + 12*b105 =L= 12; e72.. x72 + 12*b106 =L= 12; e73.. x73 + 12*b107 =L= 12; e74.. x74 + 12*b108 =L= 12; e75.. x75 + 12*b109 =L= 12; e76.. x76 + 12*b110 =L= 12; e77.. x77 + 12*b111 =L= 12; e78.. x78 + 12*b112 =L= 12; e79.. x79 + 12*b113 =L= 12; e80.. x80 + 12*b114 =L= 12; e81.. x81 + 12*b115 =L= 12; e82.. x82 + 12*b116 =L= 12; e83.. x83 + 12*b117 =L= 12; e84.. x84 + 12*b118 =L= 12; e85.. x85 + 12*b119 =L= 12; e86.. x86 + 12*b120 =L= 12; e87.. x87 + 12*b121 =L= 12; e88.. x88 + 12*b122 =L= 12; e89.. x89 + 12*b123 =L= 12; e90.. x90 + 12*b124 =L= 12; e91.. x91 + 12*b125 =L= 12; e92.. x92 + 12*b126 =L= 12; e93.. x65 - 12*b99 =G= -12; e94.. x66 - 12*b100 =G= -12; e95.. x67 - 12*b101 =G= -12; e96.. x68 - 12*b102 =G= -12; e97.. x69 - 12*b103 =G= -12; e98.. x70 - 12*b104 =G= -12; e99.. x71 - 12*b105 =G= -12; e100.. x72 - 12*b106 =G= -12; e101.. x73 - 12*b107 =G= -12; e102.. x74 - 12*b108 =G= -12; e103.. x75 - 12*b109 =G= -12; e104.. x76 - 12*b110 =G= -12; e105.. x77 - 12*b111 =G= -12; e106.. x78 - 12*b112 =G= -12; e107.. x79 - 12*b113 =G= -12; e108.. x80 - 12*b114 =G= -12; e109.. x81 - 12*b115 =G= -12; e110.. x82 - 12*b116 =G= -12; e111.. x83 - 12*b117 =G= -12; e112.. x84 - 12*b118 =G= -12; e113.. x85 - 12*b119 =G= -12; e114.. x86 - 12*b120 =G= -12; e115.. x87 - 12*b121 =G= -12; e116.. x88 - 12*b122 =G= -12; e117.. x89 - 12*b123 =G= -12; e118.. x90 - 12*b124 =G= -12; e119.. x91 - 12*b125 =G= -12; e120.. x92 - 12*b126 =G= -12; e121.. - (1.02*x93*(x57 - 6.5) + 1.02*x94*(x58 - 3.25)) + x95 =E= 0; e122.. x96 - 9.11349113439539*s1s1('127') - 17.6144733325531*s1s1('128') - 32.2986551864818*s1s1('129') - 54.4931814987685*s1s1('130') - 105.323928905069*s1s1('131') - 177.698914733437*s1s1('132') - 257.546555368226*s1s1('133') - 7.65172765642961*s1s2('134') - 14.7891900880288*s1s2('135') - 27.118094428506*s1s2('136') - 45.7527173518919*s1s2('137') - 88.4304387640365*s1s2('138') - 149.196798497086*s1s2('139') - 216.237232413786*s1s2('140') - 14.9380525029139*s1s3('141') - 28.8721329260735*s1s3('142') - 52.941183552398*s1s3('143') - 89.3205462402005*s1s3('144') - 172.637944844116*s1s3('145') - 291.268810037089*s1s3('146') - 422.148209648796*s1s3('147') - 11.9558099050809*s1s4('148') - 23.1080813747994*s1s4('149') - 42.3719709499612*s1s4('150') - 71.4885338137291*s1s4('151') - 138.172392322055*s1s4('152') - 233.119713791557*s1s4('153') - 337.870264236031*s1s4('154') - 13.9253546563734*s1s5('155') - 26.9147996770731*s1s5('156') - 49.3521332015331*s1s5('157') - 83.2652237802191*s1s5('158') - 160.93427229773*s1s5('159') - 271.522775764452*s1s5('160') - 393.529446744536*s1s5('161') - 7.76158051882097*s1s6('162') - 15.0015127080393*s1s6('163') - 27.5074183079396*s1s6('164') - 46.4095712271164*s1s6('165') - 89.7*s1s6('166') - 151.338758602103*s1s6('167') - 219.341665817957*s1s6('168') - 7.96556922221359*s1s7('169') - 15.3957802311063*s1s7('170') - 28.2303641796868*s1s7('171') - 47.6293006671023*s1s7('172') - 92.0574820424717*s1s7('173') - 155.316221319321*s1s7('174') - 225.10637081608*s1s7('175') - 9.11349113439539*s1s8('176') - 17.6144733325531*s1s8('177') - 32.2986551864818*s1s8('178') - 54.4931814987685*s1s8('179') - 105.323928905069*s1s8('180') - 177.698914733437*s1s8('181') - 257.546555368226*s1s8('182') - 11.9558099050809*s1s9('183') - 23.1080813747994*s1s9('184') - 42.3719709499612*s1s9('185') - 71.4885338137291*s1s9('186') - 138.172392322055*s1s9('187') - 233.119713791557*s1s9('188') - 337.870264236031*s1s9('189') - 7.48690188831565*s1s10('190') - 14.4706163324673*s1s10('191') - 26.5339439013751*s1s10('192') - 44.7671586494086*s1s10('193') - 86.5255598074927*s1s10('194') - 145.982952158506*s1s10('195') - 211.579268940989*s1s10('196') - 9.28783513744935*s1s11('197') - 17.9514438466182*s1s11('198') - 32.916538800503*s1s11('199') - 55.5356535066454*s1s11('200') - 107.338809384118*s1s11('201') - 181.098351861986*s1s11('202') - 262.473503425068*s1s11('203') - 7.76158051882097*s1s12('204') - 15.0015127080393*s1s12('205') - 27.5074183079396*s1s12('206') - 46.4095712271164*s1s12('207') - 89.7*s1s12('208') - 151.338758602103*s1s12('209') - 219.341665817957*s1s12('210') - 12.4236944883441*s1s13('211') - 24.0124044704238*s1s13('212') - 44.0301766363479*s1s13('213') - 74.2862014846846*s1s13('214') - 143.579699122125*s1s13('215') - 242.242736071415*s1s13('216') - 351.092646411238*s1s13('217') - 7.65172765642961*s1s14('218') - 14.7891900880288*s1s14('219') - 27.118094428506*s1s14('220') - 45.7527173518919*s1s14('221') - 88.4304387640365*s1s14('222') - 149.196798497086*s1s14('223') - 216.237232413786*s1s14('224') - 7.48690188831565*s1s15('225') - 14.4706163324673*s1s15('226') - 26.5339439013751*s1s15('227') - 44.7671586494086*s1s15('228') - 86.5255598074927*s1s15('229') - 145.982952158506*s1s15('230') - 211.579268940989*s1s15('231') - 7.22029184733547*s1s16('232') - 13.9553148538372*s1s16('233') - 25.5890649679471*s1s16('234') - 43.1729913716576*s1s16('235') - 83.44436769489*s1s16('236') - 140.784470672041*s1s16('237') - 204.044889780639*s1s16('238') - 9.28783513744935*s1s17('239') - 17.9514438466182*s1s17('240') - 32.916538800503*s1s17('241') - 55.5356535066454*s1s17('242') - 107.338809384118*s1s17('243') - 181.098351861986*s1s17('244') - 262.473503425068*s1s17('245') - 7.22029184733547*s1s18('246') - 13.9553148538372*s1s18('247') - 25.5890649679471*s1s18('248') - 43.1729913716576*s1s18('249') - 83.44436769489*s1s18('250') - 140.784470672041*s1s18('251') - 204.044889780639*s1s18('252') - 6.67516217420068*s1s19('253') - 12.9016931463472*s1s19('254') - 23.6570989315674*s1s19('255') - 39.913444642481*s1s19('256') - 77.1443452237428*s1s19('257') - 130.155289178744*s1s19('258') - 188.639567333459*s1s19('259') - 13.9253546563734*s1s20('260') - 26.9147996770731*s1s20('261') - 49.3521332015331*s1s20('262') - 83.2652237802191*s1s20('263') - 160.93427229773*s1s20('264') - 271.522775764452*s1s20('265') - 393.529446744536*s1s20('266') - 7.76158051882097*s1s21('267') - 15.0015127080393*s1s21('268') - 27.5074183079396*s1s21('269') - 46.4095712271164*s1s21('270') - 89.7*s1s21('271') - 151.338758602103*s1s21('272') - 219.341665817957*s1s21('273') - 6.67516217420068*s1s22('274') - 12.9016931463472*s1s22('275') - 23.6570989315674*s1s22('276') - 39.913444642481*s1s22('277') - 77.1443452237428*s1s22('278') - 130.155289178744*s1s22('279') - 188.639567333459*s1s22('280') - 11.9558099050809*s1s23('281') - 23.1080813747994*s1s23('282') - 42.3719709499612*s1s23('283') - 71.4885338137291*s1s23('284') - 138.172392322055*s1s23('285') - 233.119713791557*s1s23('286') - 337.870264236031*s1s23('287') - 7.76158051882097*s1s24('288') - 15.0015127080393*s1s24('289') - 27.5074183079396*s1s24('290') - 46.4095712271164*s1s24('291') - 89.7*s1s24('292') - 151.338758602103*s1s24('293') - 219.341665817957*s1s24('294') - 11.9558099050809*s1s25('295') - 23.1080813747994*s1s25('296') - 42.3719709499612*s1s25('297') - 71.4885338137291*s1s25('298') - 138.172392322055*s1s25('299') - 233.119713791557*s1s25('300') - 337.870264236031*s1s25('301') - 14.9380525029139*s1s26('302') - 28.8721329260735*s1s26('303') - 52.941183552398*s1s26('304') - 89.3205462402005*s1s26('305') - 172.637944844116*s1s26('306') - 291.268810037089*s1s26('307') - 422.148209648796*s1s26('308') - 7.96556922221359*s1s27('309') - 15.3957802311063*s1s27('310') - 28.2303641796868*s1s27('311') - 47.6293006671023*s1s27('312') - 92.0574820424717*s1s27('313') - 155.316221319321*s1s27('314') - 225.10637081608*s1s27('315') - 12.4236944883441*s1s28('316') - 24.0124044704238*s1s28('317') - 44.0301766363479*s1s28('318') - 74.2862014846846*s1s28('319') - 143.579699122125*s1s28('320') - 242.242736071415*s1s28('321') - 351.092646411238*s1s28('322') =E= 0; e123.. - 0.2*x93 - 0.17*x94 + x97 =E= 0; e124.. - 10*x95 - x96 - 10*x97 + objvar =E= 0; e125.. - b99 + s1s1('127') + s1s1('128') + s1s1('129') + s1s1('130') + s1s1('131') + s1s1('132') + s1s1('133') =E= 0; e126.. - b100 + s1s2('134') + s1s2('135') + s1s2('136') + s1s2('137') + s1s2('138') + s1s2('139') + s1s2('140') =E= 0; e127.. - b101 + s1s3('141') + s1s3('142') + s1s3('143') + s1s3('144') + s1s3('145') + s1s3('146') + s1s3('147') =E= 0; e128.. - b102 + s1s4('148') + s1s4('149') + s1s4('150') + s1s4('151') + s1s4('152') + s1s4('153') + s1s4('154') =E= 0; e129.. - b103 + s1s5('155') + s1s5('156') + s1s5('157') + s1s5('158') + s1s5('159') + s1s5('160') + s1s5('161') =E= 0; e130.. - b104 + s1s6('162') + s1s6('163') + s1s6('164') + s1s6('165') + s1s6('166') + s1s6('167') + s1s6('168') =E= 0; e131.. - b105 + s1s7('169') + s1s7('170') + s1s7('171') + s1s7('172') + s1s7('173') + s1s7('174') + s1s7('175') =E= 0; e132.. - b106 + s1s8('176') + s1s8('177') + s1s8('178') + s1s8('179') + s1s8('180') + s1s8('181') + s1s8('182') =E= 0; e133.. - b107 + s1s9('183') + s1s9('184') + s1s9('185') + s1s9('186') + s1s9('187') + s1s9('188') + s1s9('189') =E= 0; e134.. - b108 + s1s10('190') + s1s10('191') + s1s10('192') + s1s10('193') + s1s10('194') + s1s10('195') + s1s10('196') =E= 0; e135.. - b109 + s1s11('197') + s1s11('198') + s1s11('199') + s1s11('200') + s1s11('201') + s1s11('202') + s1s11('203') =E= 0; e136.. - b110 + s1s12('204') + s1s12('205') + s1s12('206') + s1s12('207') + s1s12('208') + s1s12('209') + s1s12('210') =E= 0; e137.. - b111 + s1s13('211') + s1s13('212') + s1s13('213') + s1s13('214') + s1s13('215') + s1s13('216') + s1s13('217') =E= 0; e138.. - b112 + s1s14('218') + s1s14('219') + s1s14('220') + s1s14('221') + s1s14('222') + s1s14('223') + s1s14('224') =E= 0; e139.. - b113 + s1s15('225') + s1s15('226') + s1s15('227') + s1s15('228') + s1s15('229') + s1s15('230') + s1s15('231') =E= 0; e140.. - b114 + s1s16('232') + s1s16('233') + s1s16('234') + s1s16('235') + s1s16('236') + s1s16('237') + s1s16('238') =E= 0; e141.. - b115 + s1s17('239') + s1s17('240') + s1s17('241') + s1s17('242') + s1s17('243') + s1s17('244') + s1s17('245') =E= 0; e142.. - b116 + s1s18('246') + s1s18('247') + s1s18('248') + s1s18('249') + s1s18('250') + s1s18('251') + s1s18('252') =E= 0; e143.. - b117 + s1s19('253') + s1s19('254') + s1s19('255') + s1s19('256') + s1s19('257') + s1s19('258') + s1s19('259') =E= 0; e144.. - b118 + s1s20('260') + s1s20('261') + s1s20('262') + s1s20('263') + s1s20('264') + s1s20('265') + s1s20('266') =E= 0; e145.. - b119 + s1s21('267') + s1s21('268') + s1s21('269') + s1s21('270') + s1s21('271') + s1s21('272') + s1s21('273') =E= 0; e146.. - b120 + s1s22('274') + s1s22('275') + s1s22('276') + s1s22('277') + s1s22('278') + s1s22('279') + s1s22('280') =E= 0; e147.. - b121 + s1s23('281') + s1s23('282') + s1s23('283') + s1s23('284') + s1s23('285') + s1s23('286') + s1s23('287') =E= 0; e148.. - b122 + s1s24('288') + s1s24('289') + s1s24('290') + s1s24('291') + s1s24('292') + s1s24('293') + s1s24('294') =E= 0; e149.. - b123 + s1s25('295') + s1s25('296') + s1s25('297') + s1s25('298') + s1s25('299') + s1s25('300') + s1s25('301') =E= 0; e150.. - b124 + s1s26('302') + s1s26('303') + s1s26('304') + s1s26('305') + s1s26('306') + s1s26('307') + s1s26('308') =E= 0; e151.. - b125 + s1s27('309') + s1s27('310') + s1s27('311') + s1s27('312') + s1s27('313') + s1s27('314') + s1s27('315') =E= 0; e152.. - b126 + s1s28('316') + s1s28('317') + s1s28('318') + s1s28('319') + s1s28('320') + s1s28('321') + s1s28('322') =E= 0; e153.. x1 - 0.0176041976445841*s1s1('127') - 0.0686820348432157*s1s1('128') - 0.240338257044582*s1s1('129') - 0.708118780382974*s1s1('130') - 2*s1s1('131') - 2*s1s1('132') - 2*s1s1('133') =L= 0; e154.. x2 - 0.0192122784105588*s1s2('134') - 0.0749558941482069*s1s2('135') - 0.262292300976835*s1s2('136') - 0.772802909347502*s1s2('137') - 2*s1s2('138') - 2*s1s2('139') - 2*s1s2('140') =L= 0; e155.. x3 - 0.0137502828767635*s1s3('141') - 0.0536461488738445*s1s3('142') - 0.187723353667753*s1s3('143') - 0.553097263345606*s1s3('144') - 2*s1s3('145') - 2*s1s3('146') - 2*s1s3('147') =L= 0; e156.. x4 - 0.0153698320860398*s1s4('148') - 0.0599647518268192*s1s4('149') - 0.209833968534382*s1s4('150') - 0.618242703881818*s1s4('151') - 2*s1s4('152') - 2*s1s4('153') - 2*s1s4('154') =L= 0; e157.. x5 - 0.0142414920290718*s1s5('155') - 0.0555625806701283*s1s5('156') - 0.194429501479406*s1s5('157') - 0.572855870518057*s1s5('158') - 2*s1s5('159') - 2*s1s5('160') - 2*s1s5('161') =L= 0; e158.. x6 - 0.0190758342372385*s1s6('162') - 0.0744235629590588*s1s6('163') - 0.260429520550158*s1s6('164') - 0.767314520523847*s1s6('165') - 2*s1s6('166') - 2*s1s6('167') - 2*s1s6('168') =L= 0; e159.. x7 - 0.0188299954674205*s1s7('169') - 0.0734644333642121*s1s7('170') - 0.257073249355929*s1s7('171') - 0.757425796631457*s1s7('172') - 2*s1s7('173') - 2*s1s7('174') - 2*s1s7('175') =L= 0; e160.. x8 - 0.0176041976445841*s1s8('176') - 0.0686820348432157*s1s8('177') - 0.240338257044582*s1s8('178') - 0.708118780382974*s1s8('179') - 2*s1s8('180') - 2*s1s8('181') - 2*s1s8('182') =L= 0; e161.. x9 - 0.0153698320860398*s1s9('183') - 0.0599647518268192*s1s9('184') - 0.209833968534382*s1s9('185') - 0.618242703881818*s1s9('186') - 2*s1s9('187') - 2*s1s9('188') - 2*s1s9('189') =L= 0; e162.. x10 - 0.0194226083350049*s1s10('190') - 0.0757764874800376*s1s10('191') - 0.265163793814297*s1s10('192') - 0.781263310246409*s1s10('193') - 2*s1s10('194') - 2*s1s10('195') - 2*s1s10('196') =L= 0; e163.. x11 - 0.0174381887671401*s1s11('197') - 0.0680343582075014*s1s11('198') - 0.238071849619242*s1s11('199') - 0.701441168247406*s1s11('200') - 2*s1s11('201') - 2*s1s11('202') - 2*s1s11('203') =L= 0; e164.. x12 - 0.0190758342372385*s1s12('204') - 0.0744235629590588*s1s12('205') - 0.260429520550158*s1s12('206') - 0.767314520523847*s1s12('207') - 2*s1s12('208') - 2*s1s12('209') - 2*s1s12('210') =L= 0; e165.. x13 - 0.0150776355652448*s1s13('211') - 0.0588247594211735*s1s13('212') - 0.205844806180028*s1s13('213') - 0.606489265973719*s1s13('214') - 2*s1s13('215') - 2*s1s13('216') - 2*s1s13('217') =L= 0; e166.. x14 - 0.0192122784105588*s1s14('218') - 0.0749558941482069*s1s14('219') - 0.262292300976835*s1s14('220') - 0.772802909347502*s1s14('221') - 2*s1s14('222') - 2*s1s14('223') - 2*s1s14('224') =L= 0; e167.. x15 - 0.0194226083350049*s1s15('225') - 0.0757764874800376*s1s15('226') - 0.265163793814297*s1s15('227') - 0.781263310246409*s1s15('228') - 2*s1s15('229') - 2*s1s15('230') - 2*s1s15('231') =L= 0; e168.. x16 - 0.0197779487583483*s1s16('232') - 0.0771628331590627*s1s16('233') - 0.270015017353593*s1s16('234') - 0.795556675515238*s1s16('235') - 2*s1s16('236') - 2*s1s16('237') - 2*s1s16('238') =L= 0; e169.. x17 - 0.0174381887671401*s1s17('239') - 0.0680343582075014*s1s17('240') - 0.238071849619242*s1s17('241') - 0.701441168247406*s1s17('242') - 2*s1s17('243') - 2*s1s17('244') - 2*s1s17('245') =L= 0; e170.. x18 - 0.0197779487583483*s1s18('246') - 0.0771628331590627*s1s18('247') - 0.270015017353593*s1s18('248') - 0.795556675515238*s1s18('249') - 2*s1s18('250') - 2*s1s18('251') - 2*s1s18('252') =L= 0; e171.. x19 - 0.02056968839856*s1s19('253') - 0.0802517719822704*s1s19('254') - 0.280824105561038*s1s19('255') - 0.827403949655566*s1s19('256') - 2*s1s19('257') - 2*s1s19('258') - 2*s1s19('259') =L= 0; e172.. x20 - 0.0142414920290718*s1s20('260') - 0.0555625806701283*s1s20('261') - 0.194429501479406*s1s20('262') - 0.572855870518057*s1s20('263') - 2*s1s20('264') - 2*s1s20('265') - 2*s1s20('266') =L= 0; e173.. x21 - 0.0190758342372385*s1s21('267') - 0.0744235629590588*s1s21('268') - 0.260429520550158*s1s21('269') - 0.767314520523847*s1s21('270') - 2*s1s21('271') - 2*s1s21('272') - 2*s1s21('273') =L= 0; e174.. x22 - 0.02056968839856*s1s22('274') - 0.0802517719822704*s1s22('275') - 0.280824105561038*s1s22('276') - 0.827403949655566*s1s22('277') - 2*s1s22('278') - 2*s1s22('279') - 2*s1s22('280') =L= 0; e175.. x23 - 0.0153698320860398*s1s23('281') - 0.0599647518268192*s1s23('282') - 0.209833968534382*s1s23('283') - 0.618242703881818*s1s23('284') - 2*s1s23('285') - 2*s1s23('286') - 2*s1s23('287') =L= 0; e176.. x24 - 0.0190758342372385*s1s24('288') - 0.0744235629590588*s1s24('289') - 0.260429520550158*s1s24('290') - 0.767314520523847*s1s24('291') - 2*s1s24('292') - 2*s1s24('293') - 2*s1s24('294') =L= 0; e177.. x25 - 0.0153698320860398*s1s25('295') - 0.0599647518268192*s1s25('296') - 0.209833968534382*s1s25('297') - 0.618242703881818*s1s25('298') - 2*s1s25('299') - 2*s1s25('300') - 2*s1s25('301') =L= 0; e178.. x26 - 0.0137502828767635*s1s26('302') - 0.0536461488738445*s1s26('303') - 0.187723353667753*s1s26('304') - 0.553097263345606*s1s26('305') - 2*s1s26('306') - 2*s1s26('307') - 2*s1s26('308') =L= 0; e179.. x27 - 0.0188299954674205*s1s27('309') - 0.0734644333642121*s1s27('310') - 0.257073249355929*s1s27('311') - 0.757425796631457*s1s27('312') - 2*s1s27('313') - 2*s1s27('314') - 2*s1s27('315') =L= 0; e180.. x28 - 0.0150776355652448*s1s28('316') - 0.0588247594211735*s1s28('317') - 0.205844806180028*s1s28('318') - 0.606489265973719*s1s28('319') - 2*s1s28('320') - 2*s1s28('321') - 2*s1s28('322') =L= 0; e181.. b99 + b106 =L= 1; e182.. b100 + b112 =L= 1; e183.. b101 + b124 =L= 1; e184.. b102 + b107 =L= 1; e185.. b103 + b118 =L= 1; e186.. b104 + b122 =L= 1; e187.. b105 + b125 =L= 1; e188.. b108 + b113 =L= 1; e189.. b109 + b115 =L= 1; e190.. b110 + b119 =L= 1; e191.. b111 + b126 =L= 1; e192.. b114 + b116 =L= 1; e193.. b117 + b120 =L= 1; e194.. b121 + b123 =L= 1; e195.. b99 + b100 + b101 =G= 2; e196.. b102 + b103 + b104 + b105 =G= 2; e197.. b99 + b102 + b108 + b109 + b110 + b111 =G= 2; e198.. b100 + b108 + b116 =G= 2; e199.. b109 + b116 + b120 =G= 2; e200.. b103 + b110 + b120 + b121 =G= 2; e201.. b104 + b121 =G= 2; e202.. b101 + b105 + b111 =G= 2; * set non default bounds x57.lo = 6.5; x58.lo = 3.25; x59.lo = 16.58; x60.lo = 14.92; x61.lo = 12.925; x62.lo = 12.26; x63.lo = 8.76; x64.lo = 16.08; x93.up = 2.5; x94.up = 6; $if set nostart $goto modeldef * set non default levels x57.l = 11.5; x58.l = 8.25; x59.l = 21.58; x60.l = 19.92; x61.l = 17.925; x62.l = 17.26; x63.l = 13.76; x64.l = 21.08; x93.l = 0.961470588235294; x94.l = 2.30752941176471; b99.l = 0.5; b100.l = 0.5; b101.l = 0.5; b102.l = 0.5; b103.l = 0.5; b104.l = 0.5; b105.l = 0.5; b106.l = 0.5; b107.l = 0.5; b108.l = 0.5; b109.l = 0.5; b110.l = 0.5; b111.l = 0.5; b112.l = 0.5; b113.l = 0.5; b114.l = 0.5; b115.l = 0.5; b116.l = 0.5; b117.l = 0.5; b118.l = 0.5; b119.l = 0.5; b120.l = 0.5; b121.l = 0.5; b122.l = 0.5; b123.l = 0.5; b124.l = 0.5; b125.l = 0.5; b126.l = 0.5; s1s1.l('127') = 0.0714285714285714; s1s1.l('128') = 0.0714285714285714; s1s1.l('129') = 0.0714285714285714; s1s1.l('130') = 0.0714285714285714; s1s1.l('131') = 0.0714285714285714; s1s1.l('132') = 0.0714285714285714; s1s1.l('133') = 0.0714285714285714; s1s2.l('134') = 0.0714285714285714; s1s2.l('135') = 0.0714285714285714; s1s2.l('136') = 0.0714285714285714; s1s2.l('137') = 0.0714285714285714; s1s2.l('138') = 0.0714285714285714; s1s2.l('139') = 0.0714285714285714; s1s2.l('140') = 0.0714285714285714; s1s3.l('141') = 0.0714285714285714; s1s3.l('142') = 0.0714285714285714; s1s3.l('143') = 0.0714285714285714; s1s3.l('144') = 0.0714285714285714; s1s3.l('145') = 0.0714285714285714; s1s3.l('146') = 0.0714285714285714; s1s3.l('147') = 0.0714285714285714; s1s4.l('148') = 0.0714285714285714; s1s4.l('149') = 0.0714285714285714; s1s4.l('150') = 0.0714285714285714; s1s4.l('151') = 0.0714285714285714; s1s4.l('152') = 0.0714285714285714; s1s4.l('153') = 0.0714285714285714; s1s4.l('154') = 0.0714285714285714; s1s5.l('155') = 0.0714285714285714; s1s5.l('156') = 0.0714285714285714; s1s5.l('157') = 0.0714285714285714; s1s5.l('158') = 0.0714285714285714; s1s5.l('159') = 0.0714285714285714; s1s5.l('160') = 0.0714285714285714; s1s5.l('161') = 0.0714285714285714; s1s6.l('162') = 0.0714285714285714; s1s6.l('163') = 0.0714285714285714; s1s6.l('164') = 0.0714285714285714; s1s6.l('165') = 0.0714285714285714; s1s6.l('166') = 0.0714285714285714; s1s6.l('167') = 0.0714285714285714; s1s6.l('168') = 0.0714285714285714; s1s7.l('169') = 0.0714285714285714; s1s7.l('170') = 0.0714285714285714; s1s7.l('171') = 0.0714285714285714; s1s7.l('172') = 0.0714285714285714; s1s7.l('173') = 0.0714285714285714; s1s7.l('174') = 0.0714285714285714; s1s7.l('175') = 0.0714285714285714; s1s10.l('190') = 0.0714285714285714; s1s10.l('191') = 0.0714285714285714; s1s10.l('192') = 0.0714285714285714; s1s10.l('193') = 0.0714285714285714; s1s10.l('194') = 0.0714285714285714; s1s10.l('195') = 0.0714285714285714; s1s10.l('196') = 0.0714285714285714; s1s11.l('197') = 0.0714285714285714; s1s11.l('198') = 0.0714285714285714; s1s11.l('199') = 0.0714285714285714; s1s11.l('200') = 0.0714285714285714; s1s11.l('201') = 0.0714285714285714; s1s11.l('202') = 0.0714285714285714; s1s11.l('203') = 0.0714285714285714; s1s12.l('204') = 0.0714285714285714; s1s12.l('205') = 0.0714285714285714; s1s12.l('206') = 0.0714285714285714; s1s12.l('207') = 0.0714285714285714; s1s12.l('208') = 0.0714285714285714; s1s12.l('209') = 0.0714285714285714; s1s12.l('210') = 0.0714285714285714; s1s13.l('211') = 0.0714285714285714; s1s13.l('212') = 0.0714285714285714; s1s13.l('213') = 0.0714285714285714; s1s13.l('214') = 0.0714285714285714; s1s13.l('215') = 0.0714285714285714; s1s13.l('216') = 0.0714285714285714; s1s13.l('217') = 0.0714285714285714; s1s18.l('246') = 0.0714285714285714; s1s18.l('247') = 0.0714285714285714; s1s18.l('248') = 0.0714285714285714; s1s18.l('249') = 0.0714285714285714; s1s18.l('250') = 0.0714285714285714; s1s18.l('251') = 0.0714285714285714; s1s18.l('252') = 0.0714285714285714; s1s22.l('274') = 0.0714285714285714; s1s22.l('275') = 0.0714285714285714; s1s22.l('276') = 0.0714285714285714; s1s22.l('277') = 0.0714285714285714; s1s22.l('278') = 0.0714285714285714; s1s22.l('279') = 0.0714285714285714; s1s22.l('280') = 0.0714285714285714; s1s23.l('281') = 0.0714285714285714; s1s23.l('282') = 0.0714285714285714; s1s23.l('283') = 0.0714285714285714; s1s23.l('284') = 0.0714285714285714; s1s23.l('285') = 0.0714285714285714; s1s23.l('286') = 0.0714285714285714; s1s23.l('287') = 0.0714285714285714; * set non default marginals * set non default priorities s1s1.prior('127') = 3; s1s1.prior('128') = 3; s1s1.prior('129') = 3; s1s1.prior('130') = 3; s1s1.prior('131') = 3; s1s1.prior('132') = 3; s1s1.prior('133') = 3; s1s2.prior('134') = 3; s1s2.prior('135') = 3; s1s2.prior('136') = 3; s1s2.prior('137') = 3; s1s2.prior('138') = 3; s1s2.prior('139') = 3; s1s2.prior('140') = 3; s1s3.prior('141') = 3; s1s3.prior('142') = 3; s1s3.prior('143') = 3; s1s3.prior('144') = 3; s1s3.prior('145') = 3; s1s3.prior('146') = 3; s1s3.prior('147') = 3; s1s4.prior('148') = 3; s1s4.prior('149') = 3; s1s4.prior('150') = 3; s1s4.prior('151') = 3; s1s4.prior('152') = 3; s1s4.prior('153') = 3; s1s4.prior('154') = 3; s1s5.prior('155') = 3; s1s5.prior('156') = 3; s1s5.prior('157') = 3; s1s5.prior('158') = 3; s1s5.prior('159') = 3; s1s5.prior('160') = 3; s1s5.prior('161') = 3; s1s6.prior('162') = 3; s1s6.prior('163') = 3; s1s6.prior('164') = 3; s1s6.prior('165') = 3; s1s6.prior('166') = 3; s1s6.prior('167') = 3; s1s6.prior('168') = 3; s1s7.prior('169') = 3; s1s7.prior('170') = 3; s1s7.prior('171') = 3; s1s7.prior('172') = 3; s1s7.prior('173') = 3; s1s7.prior('174') = 3; s1s7.prior('175') = 3; s1s8.prior('176') = 3; s1s8.prior('177') = 3; s1s8.prior('178') = 3; s1s8.prior('179') = 3; s1s8.prior('180') = 3; s1s8.prior('181') = 3; s1s8.prior('182') = 3; s1s9.prior('183') = 3; s1s9.prior('184') = 3; s1s9.prior('185') = 3; s1s9.prior('186') = 3; s1s9.prior('187') = 3; s1s9.prior('188') = 3; s1s9.prior('189') = 3; s1s10.prior('190') = 3; s1s10.prior('191') = 3; s1s10.prior('192') = 3; s1s10.prior('193') = 3; s1s10.prior('194') = 3; s1s10.prior('195') = 3; s1s10.prior('196') = 3; s1s11.prior('197') = 3; s1s11.prior('198') = 3; s1s11.prior('199') = 3; s1s11.prior('200') = 3; s1s11.prior('201') = 3; s1s11.prior('202') = 3; s1s11.prior('203') = 3; s1s12.prior('204') = 3; s1s12.prior('205') = 3; s1s12.prior('206') = 3; s1s12.prior('207') = 3; s1s12.prior('208') = 3; s1s12.prior('209') = 3; s1s12.prior('210') = 3; s1s13.prior('211') = 3; s1s13.prior('212') = 3; s1s13.prior('213') = 3; s1s13.prior('214') = 3; s1s13.prior('215') = 3; s1s13.prior('216') = 3; s1s13.prior('217') = 3; s1s14.prior('218') = 3; s1s14.prior('219') = 3; s1s14.prior('220') = 3; s1s14.prior('221') = 3; s1s14.prior('222') = 3; s1s14.prior('223') = 3; s1s14.prior('224') = 3; s1s15.prior('225') = 3; s1s15.prior('226') = 3; s1s15.prior('227') = 3; s1s15.prior('228') = 3; s1s15.prior('229') = 3; s1s15.prior('230') = 3; s1s15.prior('231') = 3; s1s16.prior('232') = 3; s1s16.prior('233') = 3; s1s16.prior('234') = 3; s1s16.prior('235') = 3; s1s16.prior('236') = 3; s1s16.prior('237') = 3; s1s16.prior('238') = 3; s1s17.prior('239') = 3; s1s17.prior('240') = 3; s1s17.prior('241') = 3; s1s17.prior('242') = 3; s1s17.prior('243') = 3; s1s17.prior('244') = 3; s1s17.prior('245') = 3; s1s18.prior('246') = 3; s1s18.prior('247') = 3; s1s18.prior('248') = 3; s1s18.prior('249') = 3; s1s18.prior('250') = 3; s1s18.prior('251') = 3; s1s18.prior('252') = 3; s1s19.prior('253') = 3; s1s19.prior('254') = 3; s1s19.prior('255') = 3; s1s19.prior('256') = 3; s1s19.prior('257') = 3; s1s19.prior('258') = 3; s1s19.prior('259') = 3; s1s20.prior('260') = 3; s1s20.prior('261') = 3; s1s20.prior('262') = 3; s1s20.prior('263') = 3; s1s20.prior('264') = 3; s1s20.prior('265') = 3; s1s20.prior('266') = 3; s1s21.prior('267') = 3; s1s21.prior('268') = 3; s1s21.prior('269') = 3; s1s21.prior('270') = 3; s1s21.prior('271') = 3; s1s21.prior('272') = 3; s1s21.prior('273') = 3; s1s22.prior('274') = 3; s1s22.prior('275') = 3; s1s22.prior('276') = 3; s1s22.prior('277') = 3; s1s22.prior('278') = 3; s1s22.prior('279') = 3; s1s22.prior('280') = 3; s1s23.prior('281') = 3; s1s23.prior('282') = 3; s1s23.prior('283') = 3; s1s23.prior('284') = 3; s1s23.prior('285') = 3; s1s23.prior('286') = 3; s1s23.prior('287') = 3; s1s24.prior('288') = 3; s1s24.prior('289') = 3; s1s24.prior('290') = 3; s1s24.prior('291') = 3; s1s24.prior('292') = 3; s1s24.prior('293') = 3; s1s24.prior('294') = 3; s1s25.prior('295') = 3; s1s25.prior('296') = 3; s1s25.prior('297') = 3; s1s25.prior('298') = 3; s1s25.prior('299') = 3; s1s25.prior('300') = 3; s1s25.prior('301') = 3; s1s26.prior('302') = 3; s1s26.prior('303') = 3; s1s26.prior('304') = 3; s1s26.prior('305') = 3; s1s26.prior('306') = 3; s1s26.prior('307') = 3; s1s26.prior('308') = 3; s1s27.prior('309') = 3; s1s27.prior('310') = 3; s1s27.prior('311') = 3; s1s27.prior('312') = 3; s1s27.prior('313') = 3; s1s27.prior('314') = 3; s1s27.prior('315') = 3; s1s28.prior('316') = 3; s1s28.prior('317') = 3; s1s28.prior('318') = 3; s1s28.prior('319') = 3; s1s28.prior('320') = 3; s1s28.prior('321') = 3; s1s28.prior('322') = 3; $label modeldef Model m / all /; m.limrow=0; m.limcol=0; m.prioropt = 1; $if NOT '%gams.u1%' == '' $include '%gams.u1%' $if not set MINLP $set MINLP MINLP Solve m using %MINLP% minimizing objvar;