CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C This is an example call of MIDACO 6.0 C ------------------------------------- C C MIDACO solves Multi-Objective Mixed-Integer Non-Linear Problems: C C C Minimize F_1(X),... F_O(X) where X(1,...N-NI) is CONTINUOUS C and X(N-NI+1,...N) is DISCRETE C C subject to G_j(X) = 0 (j=1,...ME) equality constraints C G_j(X) >= 0 (j=ME+1,...M) inequality constraints C C and bounds XL <= X <= XU C C C The problem statement of this example is given below. You can use C this example as template to run your own problem. To do so: Replace C the objective functions 'F' (and in case the constraints 'G') given C here with your own problem and follow the below instruction steps. C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC CCCCCCCCCCCCCCCCCCCCCCCCC MAIN PROGRAM CCCCCCCCCCCCCCCCCCCCCCCCCCCCC CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC PROGRAM MAIN IMPLICIT NONE CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C Dimensions of the optimization problem INTEGER O, N, NI, M, ME C Lower and upper bounds ('XL' & 'XU') and optimization variable 'X' DOUBLE PRECISION XL(1000), XU(1000), X(1000) C Objectives 'F(X)' and constraints 'G(X)' DOUBLE PRECISION F(10), G(1000) C MIDACO information and stop flags INTEGER IFLAG, ISTOP C MIDACO parameter DOUBLE PRECISION PARAM(13) C MIDACO integer 'IW' and real'RW' workspace and pareto front 'PF' INTEGER LIW, LRW, LPF PARAMETER (LIW = 50000, LRW = 50000, LPF = 50000) INTEGER IW(LIW) DOUBLE PRECISION RW(LRW),PF(LPF) C Parameter for stopping criteria, printing and license INTEGER MAXTIME, MAXEVAL, PRINTEVAL, SAVE2FILE, I CHARACTER*60 KEY KEY='************************************************************' CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC CCC Step 1: Problem definition CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C Step 1.A : Problem dimensions C CCCCCCCCCCCCCCCCCCCCCCCCCCCCC O = 1 ! Number of objectives N = 205 ! Number of variables (in total) NI = 25 ! Number of integer variables (0 <= NI <= N) M = 115 ! Number of constraints (in total) ME = 81 ! Number of equality constraints (0 <= ME <= M) C C Step 1.B : Lower and upper bounds: 'XL' and 'XU' C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC DO I = 1,N XL(I) = 0.0D0 XU(I) = 1.0D0 END DO XL(26) = 1.0D0 XU(26) = 100.0D0 XL(180+1) = 1.0D0 XL(180+7) = 1.0D0 XL(180+13) = 1.0D0 XL(180+19) = 1.0D0 XL(180+25) = 1.0D0 C C Step 1.C : Starting point 'X' C CCCCCCCCCCCCCCCCCCCCCCCCCCCCC DO I = 1,N X(I) = XL(I) ! Here for example: starting point = lower bounds ENDDO CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC CCC Step 2: Choose stopping criteria and printing options CCCCCCCCCCC CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C Step 2.A : Stopping criteria C CCCCCCCCCCCCCCCCCCCCCCCCCCCC MAXEVAL = 999999999 ! Maximum evaluation budget (e.g. 1000000) MAXTIME = 60*60*24 ! Maximum time limit (e.g. 60*60*24 = 1 Day) C C Step 2.B : Printing options C CCCCCCCCCCCCCCCCCCCCCCCCCCC PRINTEVAL = 1000000 ! Print-Frequency for current best solution (e.g. 1000) SAVE2FILE = 1 ! Save SCREEN and SOLUTION to TXT-files [0=NO/1=YES] CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC CCC Step 3: Choose MIDACO parameters (FOR ADVANCED USERS) CCCCCCCCCCC CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC PARAM( 1) = 0.01D0 ! ACCURACY PARAM( 2) = 0.0D0 ! SEED PARAM( 3) = 21729.0d0 ! FSTOPP PARAM( 4) = 0.0D0 ! ALGOSTOP PARAM( 5) = 0.0D0 ! EVALSTOP PARAM( 6) = 1000.0D0 ! FOCUS PARAM( 7) = 0.0D0 ! ANTS PARAM( 8) = 0.0D0 ! KERNEL PARAM( 9) = 0.0D0 ! ORACLE PARAM(10) = 0.0D0 ! PARETOMAX PARAM(11) = 0.0D0 ! EPSILON PARAM(12) = 0.0D0 ! BALANCE PARAM(13) = 0.0D0 ! CHARACTER CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C Call MIDACO by Reverse Communication C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C Print MIDACO headline with basic information CALL MIDACO_PRINT(1,PRINTEVAL,SAVE2FILE,IFLAG,ISTOP,F,G,X,XL, & XU,O,N,NI,M,ME,RW,PF,MAXEVAL,MAXTIME,PARAM,1,0,KEY) DO WHILE(ISTOP.EQ.0) !~~~Start~of~reverse~communication~loop C Evaluate Objective F(X) and constraints G(X) CALL PROBLEM_FUNCTION( F, G , X) C Call MIDACO CALL MIDACO(1,O,N,NI,M,ME,X,F,G,XL,XU,IFLAG, & ISTOP,PARAM,RW,LRW,IW,LIW,PF,LPF,KEY) C Call MIDACO printing routine CALL MIDACO_PRINT(2,PRINTEVAL,SAVE2FILE,IFLAG,ISTOP,F,G,X, & XL,XU,O,N,NI,M,ME,RW,PF,MAXEVAL,MAXTIME,PARAM,1,0,KEY) ENDDO !~~~~~~~~~~~~~~~~~~~~End~of~reverse~communication~loop CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC ! PRINT*," Solution F(1) = ", F(1) ! PRINT*," Solution G(1) = ", G(1) ! PRINT*," Solution X(1) = ", X(1) CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC END CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC CCCCCCCCCCCCCCCCCCCCCCCCCCCC END OF MAIN CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC CCCCCCCCCCCCCCCCCCCCCC OPTIMIZATION PROBLEM CCCCCCCCCCCCCCCCCCCCCCCC CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC SUBROUTINE PROBLEM_FUNCTION(F,G,X) CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC IMPLICIT NONE CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC DOUBLE PRECISION F(*),G(*),X(*) CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC DOUBLE PRECISION Y(25), & B1, B2, B3, B4, B5, B6, B7, B8, B9, B10, & B11, B12, B13, B14, B15, B16, B17, B18, B19, B20, & B21, B22, B23, B24, B25, 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, X98, X99,X100, & X101,X102,X103,X104,X105,X106,X107,X108,X109,X110, & X111,X112,X113,X114,X115,X116,X117,X118,X119,X120, & X121,X122,X123,X124,X125,X126,X127,X128,X129,X130, & X131,X132,X133,X134,X135,X136,X137,X138,X139,X140, & X141,X142,X143,X144,X145,X146,X147,X148,X149,X150, & X151,X152,X153,X154,X155,X156,X157,X158,X159,X160, & X161,X162,X163,X164,X165,X166,X167,X168,X169,X170, & X171,X172,X173,X174,X175,X176,X177,X178,X179,X180 DOUBLE PRECISION & X181,X182,X183,X184,X185,X186,X187,X188,X189,X190, & X191,X192,X193,X194,X195,X196,X197,X198,X199,X200, & X201,X202,X203,X204,X205, & 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, & 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 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC INTEGER I CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC DO I=1,25 Y(I) = X(I) END DO X26 = X(1) X27 = X(2) X28 = X(3) X29 = X(4) X30 = X(5) X31 = X(6) X32 = X(7) X33 = X(8) X34 = X(9) X35 = X(10) X36 = X(11) X37 = X(12) X38 = X(13) X39 = X(14) X40 = X(15) X41 = X(16) X42 = X(17) X43 = X(18) X44 = X(19) X45 = X(20) X46 = X(21) X47 = X(22) X48 = X(23) X49 = X(24) X50 = X(25) X51 = X(26) X52 = X(27) X53 = X(28) X54 = X(29) X55 = X(30) X56 = X(31) X57 = X(32) X58 = X(33) X59 = X(34) X60 = X(35) X61 = X(36) X62 = X(37) X63 = X(38) X64 = X(39) X65 = X(40) X66 = X(41) X67 = X(42) X68 = X(43) X69 = X(44) X70 = X(45) X71 = X(46) X72 = X(47) X73 = X(48) X74 = X(49) X75 = X(50) X76 = X(51) X77 = X(52) X78 = X(53) X79 = X(54) X80 = X(55) X81 = X(56) X82 = X(57) X83 = X(58) X84 = X(59) X85 = X(60) X86 = X(61) X87 = X(62) X88 = X(63) X89 = X(64) X90 = X(65) X91 = X(66) X92 = X(67) X93 = X(68) X94 = X(69) X95 = X(70) X96 = X(71) X97 = X(72) X98 = X(73) X99 = X(74) X100 = X(75) X101 = X(76) X102 = X(77) X103 = X(78) X104 = X(79) X105 = X(80) X106 = X(81) X107 = X(82) X108 = X(83) X109 = X(84) X110 = X(85) X111 = X(86) X112 = X(87) X113 = X(88) X114 = X(89) X115 = X(90) X116 = X(91) X117 = X(92) X118 = X(93) X119 = X(94) X120 = X(95) X121 = X(96) X122 = X(97) X123 = X(98) X124 = X(99) X125 = X(100) X126 = X(101) X127 = X(102) X128 = X(103) X129 = X(104) X130 = X(105) X131 = X(106) X132 = X(107) X133 = X(108) X134 = X(109) X135 = X(110) X136 = X(111) X137 = X(112) X138 = X(113) X139 = X(114) X140 = X(115) X141 = X(116) X142 = X(117) X143 = X(118) X144 = X(119) X145 = X(120) X146 = X(121) X147 = X(122) X148 = X(123) X149 = X(124) X150 = X(125) X151 = X(126) X152 = X(127) X153 = X(128) X154 = X(129) X155 = X(130) X156 = X(131) X157 = X(132) X158 = X(133) X159 = X(134) X160 = X(135) X161 = X(136) X162 = X(137) X163 = X(138) X164 = X(139) X165 = X(140) X166 = X(141) X167 = X(142) X168 = X(143) X169 = X(144) X170 = X(145) X171 = X(146) X172 = X(147) X173 = X(148) X174 = X(149) X175 = X(150) X176 = X(151) X177 = X(152) X178 = X(153) X179 = X(154) X180 = X(155) X181 = X(156) X182 = X(157) X183 = X(158) X184 = X(159) X185 = X(160) X186 = X(161) X187 = X(162) X188 = X(163) X189 = X(164) X190 = X(165) X191 = X(166) X192 = X(167) X193 = X(168) X194 = X(169) X195 = X(170) X196 = X(171) X197 = X(172) X198 = X(173) X199 = X(174) X200 = X(175) X201 = X(176) X202 = X(177) X203 = X(178) X204 = X(179) X205 = X(180) B1 = Y(1) B2 = Y(2) B3 = Y(3) B4 = Y(4) B5 = Y(5) B6 = Y(6) B7 = Y(7) B8 = Y(8) B9 = Y(9) B10 = Y(10) B11 = Y(11) B12 = Y(12) B13 = Y(13) B14 = Y(14) B15 = Y(15) B16 = Y(16) B17 = Y(17) B18 = Y(18) B19 = Y(19) B20 = Y(20) B21 = Y(21) B22 = Y(22) B23 = Y(23) B24 = Y(24) B25 = Y(25) F(1) = 1013.33333333333D0*X26/X51 + 1013.33333333333D0*X27/X51 / + 1013.33333333333D0*X28/X51 + 1013.33333333333D0*X29/X51 / + 1013.33333333333D0*X30/X51 + 781.25D0*X31/X51 / + 781.25D0*X32/X51 + 781.25D0*X33/X51 + 781.25D0*X34/X51 / + 781.25D0*X35/X51 + 962.5D0*X36/X51 + 962.5D0*X37/X51 / + 962.5D0*X38/X51 + 962.5D0*X39/X51 + 962.5D0*X40/X51 / + 937.5D0*X41/X51 + 937.5D0*X42/X51 + 937.5D0*X43/X51 / + 937.5D0*X44/X51 + 937.5D0*X45/X51 + 923.203125D0*X46/X51 / + 923.203125D0*X47/X51 + 923.203125D0*X48/X51 / + 923.203125D0*X49/X51 + 923.203125D0*X50/X51 / - 1406.25D0*X82/X51 - 1406.25D0*X83/X51 - 1406.25D0*X84/X51 / - 1406.25D0*X85/X51 - 1406.25D0*X86/X51 - 4331.25D0*X87/X51 / - 4331.25D0*X88/X51 - 4331.25D0*X89/X51 - 4331.25D0*X90/X51 / - 4331.25D0*X91/X51 - 3375D0*X92/X51 - 3375D0*X93/X51 / - 3375D0*X94/X51 - 3375D0*X95/X51 - 3375D0*X96/X51 / - 4154.4140625D0*X97/X51 - 4154.4140625D0*X98/X51 / - 4154.4140625D0*X99/X51 - 4154.4140625D0*X100/X51 / - 4154.4140625D0*X101/X51 - 1824D0*X102/X51 - 1824D0*X103/X51 / - 1824D0*X104/X51 - 1824D0*X105/X51 - 1824D0*X106/X51 / - 2598.75D0*X112/X51 - 2598.75D0*X113/X51 F(1) = F(1) / - 2598.75D0*X114/X51 - 2598.75D0*X115/X51 / - 2598.75D0*X116/X51 - 1687.5D0*X117/X51 - 1687.5D0*X118/X51 / - 1687.5D0*X119/X51 - 1687.5D0*X120/X51 - 1687.5D0*X121/X51 / - 2492.6484375D0*X122/X51 - 2492.6484375D0*X123/X51 / - 2492.6484375D0*X124/X51 - 2492.6484375D0*X125/X51 / - 2492.6484375D0*X126/X51 - 4560.0D0*X127/X51 / - 4560.0D0*X128/X51 - 4560D0*X129/X51 - 4560D0*X130/X51 / - 4560D0*X131/X51 - 2109.375D0*X132/X51 - 2109.375D0*X133/X51 / - 2109.375D0*X134/X51 - 2109.375D0*X135/X51 / - 2109.375D0*X136/X51 - 1687.5D0*X142/X51 - 1687.5D0*X143/X51 / - 1687.5D0*X144/X51 - 1687.5D0*X145/X51 - 1687.5D0*X146/X51 / - 4985.296875D0*X147/X51 - 4985.296875D0*X148/X51 / - 4985.296875D0*X149/X51 - 4985.296875D0*X150/X51 / - 4985.296875D0*X151/X51 - 3648D0*X152/X51 - 3648D0*X153/X51 / - 3648D0*X154/X51 - 3648D0*X155/X51 - 3648D0*X156/X51 / - 1406.25D0*X157/X51 - 1406.25D0*X158/X51 / - 1406.25D0*X159/X51 - 1406.25D0*X160/X51 / - 1406.25D0*X161/X51 - 1732.5D0*X162/X51 - 1732.5D0*X163/X51 / - 1732.5D0*X164/X51 - 1732.5D0*X165/X51 - 1732.5D0*X166/X51 F(1) = F(1) / - 4154.4140625D0*X172/X51 - 4154.4140625D0*X173/X51 / - 4154.4140625D0*X174/X51 - 4154.4140625D0*X175/X51 / - 4154.4140625D0*X176/X51 - 4560D0*X177/X51 - 4560D0*X178/X51 / - 4560D0*X179/X51 - 4560D0*X180/X51 - 4560D0*X181/X51 / - 2109.375D0*X182/X51 - 2109.375D0*X183/X51 / - 2109.375D0*X184/X51 - 2109.375D0*X185/X51 / - 2109.375D0*X186/X51 - 5197.5D0*X187/X51 / - 5197.5D0*X188/X51 - 5197.5D0*X189/X51 - 5197.5D0*X190/X51 / - 5197.5D0*X191/X51 - 4218.75D0*X192/X51 - 4218.75D0*X193/X51 / - 4218.75D0*X194/X51 - 4218.75D0*X195/X51 / - 4218.75D0*X196/X51 / - 2500D0*X72/X51 - 2500D0*X73/X51 - 2500D0*X74/X51 / - 2500D0*X75/X51 - 2500D0*X76/X51 / + 3D0*X72 + 3D0*X73 + 3D0*X74 + 3D0*X75 + 3D0*X76 F(1) = F(1) / + 2.63768115942029D0*X52 + 2.63768115942029D0*X53 / + 2.63768115942029D0*X54 + 2.63768115942029D0*X55 / + 2.63768115942029D0*X56 + 2.81449275362319D0*X57 / + 2.81449275362319D0*X58 + 2.81449275362319D0*X59 / + 2.81449275362319D0*X60 + 2.81449275362319D0*X61 / + 2.88405797101449D0*X62 + 2.88405797101449D0*X63 / + 2.88405797101449D0*X64 + 2.88405797101449D0*X65 / + 2.88405797101449D0*X66 + 2.88405797101449D0*X67 / + 2.88405797101449D0*X68 + 2.88405797101449D0*X69 / + 2.88405797101449D0*X70 + 2.88405797101449D0*X71 / + 100000D0*X202 + 100000D0*X203 + 100000D0*X204 / + 100000D0*X205 E1 = -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 E2 = -(X52/X51 + X53/X51 + X54/X51 + X55/X51 + X56/X51) - X202 / + 0.483091787439614D0 E3 = -(X57/X51 + X58/X51 + X59/X51 + X60/X51 + X61/X51) - X203 / + 0.193236714975845D0 E4 = -(X62/X51 + X63/X51 + X64/X51 + X65/X51 + X66/X51) - X204 / + 0.144927536231884D0 E5 = -(X67/X51 + X68/X51 + X69/X51 + X70/X51 + X71/X51) - X205 / + 0.144927536231884D0 E6 = -1.33333333333333D0*X26 + X52 + 2.66666666666667D0*X102 / + 6.66666666666667D0*X127 + 5.33333333333333D0*X152 / + 6.66666666666667D0*X177 E7 = -1.33333333333333D0*X27 + X53 + 2.66666666666667D0*X103 / + 6.66666666666667D0*X128 + 5.33333333333333D0*X153 / + 6.66666666666667D0*X178 E8 = -1.33333333333333D0*X28 + X54 + 2.66666666666667D0*X104 / + 6.66666666666667D0*X129 + 5.33333333333333D0*X154 / + 6.66666666666667D0*X179 E9 = -1.33333333333333D0*X29 + X55 + 2.66666666666667D0*X105 / + 6.66666666666667D0*X130 + 5.33333333333333D0*X155 / + 6.66666666666667D0*X180 E10 = -1.33333333333333D0*X30 + X56 + 2.66666666666667D0*X106 / + 6.66666666666667D0*X131 + 5.33333333333333D0*X156 / + 6.66666666666667D0*X181 E11 = -1.04166666666667D0*X31 + X57 + 2.08333333333333D0*X82 / + 3.125D0*X132 + 2.08333333333333D0*X157 + 3.125D0*X182 E12 = -1.04166666666667D0*X32 + X58 + 2.08333333333333D0*X83 / + 3.125D0*X133 + 2.08333333333333D0*X158 + 3.125D0*X183 E13 = -1.04166666666667D0*X33 + X59 + 2.08333333333333D0*X84 / + 3.125D0*X134 + 2.08333333333333D0*X159 + 3.125D0*X184 E14 = -1.04166666666667D0*X34 + X60 + 2.08333333333333D0*X85 / + 3.125D0*X135 + 2.08333333333333D0*X160 + 3.125D0*X185 E15 = -1.04166666666667D0*X35 + X61 + 2.08333333333333D0*X86 / + 3.125D0*X136 + 2.08333333333333D0*X161 + 3.125D0*X186 E16 = -1.25D0*X36 + X62 + 6.25D0*X87 + 3.75D0*X112 + 2.5D0*X162 / + 7.5D0*X187 E17 = -1.25D0*X37 + X63 + 6.25D0*X88 + 3.75D0*X113 + 2.5D0*X163 / + 7.5D0*X188 E18 = -1.25D0*X38 + X64 + 6.25D0*X89 + 3.75D0*X114 + 2.5D0*X164 / + 7.5D0*X189 E19 = -1.25D0*X39 + X65 + 6.25D0*X90 + 3.75D0*X115 + 2.5D0*X165 / + 7.5D0*X190 E20 = -1.25D0*X40 + X66 + 6.25D0*X91 + 3.75D0*X116 + 2.5D0*X166 / + 7.5D0*X191 E21 = -1.25D0*X41 + X67 + 5D0*X92 + 2.5D0*X117 + 2.5D0*X142 / + 6.25D0*X192 E22 = -1.25D0*X42 + X68 + 5D0*X93 + 2.5D0*X118 + 2.5D0*X143 / + 6.25D0*X193 E23 = -1.25D0*X43 + X69 + 5D0*X94 + 2.5D0*X119 + 2.5D0*X144 / + 6.25D0*X194 E24 = -1.25D0*X44 + X70 + 5D0*X95 + 2.5D0*X120 + 2.5D0*X145 / + 6.25D0*X195 E25 = -1.25D0*X45 + X71 + 5D0*X96 + 2.5D0*X121 + 2.5D0*X146 / + 6.25D0*X196 E26 = -1.21875D0*X46 + X72 + 6.09375D0*X97 + 3.65625D0*X122 / + 7.3125D0*X147 + 6.09375D0*X172 E27 = -1.21875D0*X47 + X73 + 6.09375D0*X98 + 3.65625D0*X123 / + 7.3125D0*X148 + 6.09375D0*X173 E28 = -1.21875D0*X48 + X74 + 6.09375D0*X99 + 3.65625D0*X124 / + 7.3125D0*X149 + 6.09375D0*X174 E29 = -1.21875D0*X49 + X75 + 6.09375D0*X100 + 3.65625D0*X125 / + 7.3125D0*X150 + 6.09375D0*X175 E30 = -1.21875D0*X50 + X76 + 6.09375D0*X101 + 3.65625D0*X126 / + 7.3125D0*X151 + 6.09375D0*X176 E31 = B1 + B6 + B11 + B16 - 1.0D0 E32 = B2 + B7 + B12 + B17 - 1.0D0 E33 = B3 + B8 + B13 + B18 - 1.0D0 E34 = B4 + B9 + B14 + B19 - 1.0D0 E35 = B1 + B2 + B3 + B4 + B5 - 1.0D0 E36 = B6 + B7 + B8 + B9 + B10 - 1.0D0 E37 = B11 + B12 + B13 + B14 + B15 - 1.0D0 E38 = B16 + B17 + B18 + B19 + B20 - 1.0D0 E39 = B21 + B22 + B23 + B24 + B25 - 1.0D0 E41 = -256.0D0*B1 + X26 E42 = -256.0D0*B6 + X27 E43 = -256.0D0*B11 + X28 E44 = -256.0D0*B16 + X29 E45 = -256.0D0*B21 + X30 E46 = -131.0D0*B2 + X31 E47 = -131.0D0*B7 + X32 E48 = -131.0D0*B12 + X33 E49 = -131.0D0*B17 + X34 E50 = -131.0D0*B22 + X35 E51 = -86.0D0*B3 + X36 E52 = -86.0D0*B8 + X37 E53 = -86.0D0*B13 + X38 E54 = -86.0D0*B18 + X39 E55 = -86.0D0*B23 + X40 E56 = -85.0D0*B4 + X41 E57 = -85.0D0*B9 + X42 E58 = -85.0D0*B14 + X43 E59 = -85.0D0*B19 + X44 E60 = -85.0D0*B24 + X45 E61 = -690.0D0*B5 + X46 E62 = -690.0D0*B10 + X47 E63 = -690.0D0*B15 + X48 E64 = -690.0D0*B20 + X49 E65 = -690.0D0*B25 + X50 E66 = -B1 + X77 + X102 + X127 + X152 + X177 E67 = -B6 + X78 + X103 + X128 + X153 + X178 E68 = -B11 + X79 + X104 + X129 + X154 + X179 E69 = -B16 + X80 + X105 + X130 + X155 + X180 E70 = -B21 + X81 + X106 + X131 + X156 + X181 E71 = -B2 + X82 + X107 + X132 + X157 + X182 E72 = -B7 + X83 + X108 + X133 + X158 + X183 E73 = -B12 + X84 + X109 + X134 + X159 + X184 E74 = -B17 + X85 + X110 + X135 + X160 + X185 E75 = -B22 + X86 + X111 + X136 + X161 + X186 E76 = -B3 + X87 + X112 + X137 + X162 + X187 E77 = -B8 + X88 + X113 + X138 + X163 + X188 E78 = -B13 + X89 + X114 + X139 + X164 + X189 E79 = -B18 + X90 + X115 + X140 + X165 + X190 E80 = -B23 + X91 + X116 + X141 + X166 + X191 E81 = -B4 + X92 + X117 + X142 + X167 + X192 E82 = -B9 + X93 + X118 + X143 + X168 + X193 E83 = -B14 + X94 + X119 + X144 + X169 + X194 E84 = -B19 + X95 + X120 + X145 + X170 + X195 E85 = -B24 + X96 + X121 + X146 + X171 + X196 E86 = -B5 + X97 + X122 + X147 + X172 + X197 E87 = -B10 + X98 + X123 + X148 + X173 + X198 E88 = -B15 + X99 + X124 + X149 + X174 + X199 E89 = -B20 + X100 + X125 + X150 + X175 + X200 E90 = -B25 + X101 + X126 + X151 + X176 + X201 E91 = -B21 + X77 + X82 + X87 + X92 + X97 E92 = -B1 + X78 + X83 + X88 + X93 + X98 E93 = -B6 + X79 + X84 + X89 + X94 + X99 E94 = -B11 + X80 + X85 + X90 + X95 + X100 E95 = -B16 + X81 + X86 + X91 + X96 + X101 E96 = -B22 + X102 + X107 + X112 + X117 + X122 E97 = -B2 + X103 + X108 + X113 + X118 + X123 E98 = -B7 + X104 + X109 + X114 + X119 + X124 E99 = -B12 + X105 + X110 + X115 + X120 + X125 E100 = -B17 + X106 + X111 + X116 + X121 + X126 E101 = -B23 + X127 + X132 + X137 + X142 + X147 E102 = -B3 + X128 + X133 + X138 + X143 + X148 E103 = -B8 + X129 + X134 + X139 + X144 + X149 E104 = -B13 + X130 + X135 + X140 + X145 + X150 E105 = -B18 + X131 + X136 + X141 + X146 + X151 E106 = -B24 + X152 + X157 + X162 + X167 + X172 E107 = -B4 + X153 + X158 + X163 + X168 + X173 E108 = -B9 + X154 + X159 + X164 + X169 + X174 E109 = -B14 + X155 + X160 + X165 + X170 + X175 E110 = -B19 + X156 + X161 + X166 + X171 + X176 E111 = -B25 + X177 + X182 + X187 + X192 + X197 E112 = -B5 + X178 + X183 + X188 + X193 + X198 E113 = -B10 + X179 + X184 + X189 + X194 + X199 E114 = -B15 + X180 + X185 + X190 + X195 + X200 E115 = -B20 + X181 + X186 + X191 + X196 + X201 G(1) = E1 G(2) = E6 G(3) = E7 G(4) = E8 G(5) = E9 G(6) = E10 G(7) = E11 G(8) = E12 G(9) = E13 G(10) = E14 G(11) = E15 G(12) = E16 G(13) = E17 G(14) = E18 G(15) = E19 G(16) = E20 G(17) = E21 G(18) = E22 G(19) = E23 G(20) = E24 G(21) = E25 G(22) = E26 G(23) = E27 G(24) = E28 G(25) = E29 G(26) = E30 G(27) = E35 G(28) = E36 G(29) = E37 G(30) = E38 G(31) = E39 G(32) = E66 G(33) = E67 G(34) = E68 G(35) = E69 G(36) = E70 G(37) = E71 G(38) = E72 G(39) = E73 G(40) = E74 G(41) = E75 G(42) = E76 G(43) = E77 G(44) = E78 G(45) = E79 G(46) = E80 G(47) = E81 G(48) = E82 G(49) = E83 G(50) = E84 G(51) = E85 G(52) = E86 G(53) = E87 G(54) = E88 G(55) = E89 G(56) = E90 G(57) = E91 G(58) = E92 G(59) = E93 G(60) = E94 G(61) = E95 G(62) = E96 G(63) = E97 G(64) = E98 G(65) = E99 G(66) = E100 G(67) = E101 G(68) = E102 G(69) = E103 G(70) = E104 G(71) = E105 G(72) = E106 G(73) = E107 G(74) = E108 G(75) = E109 G(76) = E110 G(77) = E111 G(78) = E112 G(79) = E113 G(80) = E114 G(81) = E115 G(82) = -E2 G(83) = -E3 G(84) = -E4 G(85) = -E5 G(86) = E31 G(87) = E32 G(88) = E33 G(89) = E34 G(91) = -E41 G(92) = -E42 G(93) = -E43 G(94) = -E44 G(95) = -E45 G(96) = -E46 G(97) = -E47 G(98) = -E48 G(99) = -E49 G(100) = -E50 G(101) = -E51 G(102) = -E52 G(103) = -E53 G(104) = -E54 G(105) = -E55 G(106) = -E56 G(107) = -E57 G(108) = -E58 G(109) = -E59 G(110) = -E60 G(111) = -E61 G(112) = -E62 G(113) = -E63 G(114) = -E64 G(115) = -E65 END CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC CCCCCCCCCCCCCCCCCCCCCCCCCCCC END OF FILE CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC