Set j /0*9/; Set i /i0*i29/; Scalar PI / 3.14159265359 / ; Variable x(j), z; x.lo['0'] = 0 ; x.up['0'] = 10 ; x.lo['1'] = 0 ; x.up['1'] = 10 ; x.lo['2'] = 0 ; x.up['2'] = 10 ; x.lo['3'] = 0 ; x.up['3'] = 10 ; x.lo['4'] = 0 ; x.up['4'] = 10 ; x.lo['5'] = 0 ; x.up['5'] = 10 ; x.lo['6'] = 0 ; x.up['6'] = 10 ; x.lo['7'] = 0 ; x.up['7'] = 10 ; x.lo['8'] = 0 ; x.up['8'] = 10 ; x.lo['9'] = 0 ; x.up['9'] = 10 ; Parameter c[i] / i0 0.806 , i1 0.517 , i2 0.100 , i3 0.908 , i4 0.965 , i5 0.669 , i6 0.524 , i7 0.902 , i8 0.531 , i9 0.876 , i10 0.462 , i11 0.491 , i12 0.463 , i13 0.714 , i14 0.352 , i15 0.869 , i16 0.813 , i17 0.811 , i18 0.828 , i19 0.964 , i20 0.789 , i21 0.360 , i22 0.369 , i23 0.992 , i24 0.332 , i25 0.817 , i26 0.632 , i27 0.883 , i28 0.608 , i29 0.326 / ; Parameter a[j,i] ; Table b[i,j] 0 1 2 3 4 5 6 7 8 9 i0 9.681 0.667 4.783 9.095 3.517 9.325 6.544 0.211 5.122 2.020 i1 9.400 2.041 3.788 7.931 2.882 2.672 3.568 1.284 7.033 7.374 i2 8.025 9.152 5.114 7.621 4.564 4.711 2.996 6.126 0.734 4.982 i3 2.196 0.415 5.649 6.979 9.510 9.166 6.304 6.054 9.377 1.426 i4 8.074 8.777 3.467 1.863 6.708 6.349 4.534 0.276 7.633 1.567 i5 7.650 5.658 0.720 2.764 3.278 5.283 7.474 6.274 1.409 8.208 i6 1.256 3.605 8.623 6.905 4.584 8.133 6.071 6.888 4.187 5.448 i7 8.314 2.261 4.224 1.781 4.124 0.932 8.129 8.658 1.208 5.762 i8 0.226 8.858 1.420 0.945 1.622 4.698 6.228 9.096 0.972 7.637 i9 7.305 2.228 1.242 5.928 9.133 1.826 4.060 5.204 8.713 8.247 i10 0.652 7.027 0.508 4.876 8.807 4.632 5.808 6.937 3.291 7.016 i11 2.699 3.516 5.874 4.119 4.461 7.496 8.817 0.690 6.593 9.789 i12 8.327 3.897 2.017 9.570 9.825 1.150 1.395 3.885 6.354 0.109 i13 2.132 7.006 7.136 2.641 1.882 5.943 7.273 7.691 2.880 0.564 i14 4.707 5.579 4.080 0.581 9.698 8.542 8.077 8.515 9.231 4.670 i15 8.304 7.559 8.567 0.322 7.128 8.392 1.472 8.524 2.277 7.826 i16 8.632 4.409 4.832 5.768 7.050 6.715 1.711 4.323 4.405 4.591 i17 4.887 9.112 0.170 8.967 9.693 9.867 7.508 7.770 8.382 6.740 i18 2.440 6.686 4.299 1.007 7.008 1.427 9.398 8.480 9.950 1.675 i19 6.306 8.583 6.084 1.138 4.350 3.134 7.853 6.061 7.457 2.258 i20 0.652 2.343 1.370 0.821 1.310 1.063 0.689 8.819 8.833 9.070 i21 5.558 1.272 5.756 9.857 2.279 2.764 1.284 1.677 1.244 1.234 i22 3.352 7.549 9.817 9.437 8.687 4.167 2.570 6.540 0.228 0.027 i23 8.798 0.880 2.370 0.168 1.701 3.680 1.231 2.390 2.499 0.064 i24 1.460 8.057 1.336 7.217 7.914 3.615 9.981 9.198 5.292 1.224 i25 0.432 8.645 8.774 0.249 8.081 7.461 4.416 0.652 4.002 4.644 i26 0.679 2.800 5.523 3.049 2.968 7.225 6.730 4.199 9.614 9.229 i27 4.263 1.074 7.286 5.599 8.291 5.200 9.214 8.272 4.398 4.506 i28 9.496 4.830 3.150 8.270 5.079 1.231 5.731 9.494 1.883 9.732 i29 4.138 2.562 2.532 9.661 5.611 5.500 6.886 2.341 9.699 6.500 ; a[j,i] = b[i,j] ; $macro defobj -Sum{i, 1.0/(Sum{j, Sqr(x[j]-a[j,i])} + c[i])} Equation Def_obj; Def_obj.. z =e= defobj ; Model m /all/; Solve m min z using nlp; $onDotL File fx / Shekelfox10.txt /; fx.nd=12; fx.nw=22; Put fx; Scalar cnt; for (cnt=1 to 100, x.l(j) = uniform(x.lo(j), x.up(j)); z.l = defobj ; Loop(j$(ord(j) le 10), Put x.l(j)); Put z.l /);