Set j /0*4/; 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 ; 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,*] 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 / Shekelfox5.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 5), Put x.l(j)); Put z.l /);