Set j /0*1/;

Scalar PI / 3.14159265359 / ;

Variable x(j), z;

x.lo['0'] = -1.5 ;
x.up['0'] =  4 ;
x.lo['1'] = -3 ;
x.up['1'] =  3 ;

$macro defobj Sin((x['0']+x['1']))+power((x['0']-x['1']),2)-1.5*x['0']+2.5*x['1']+1.0
Equation
         Def_obj;
Def_obj..
    z =e= defobj ;

model m /all/;

solve m min z using nlp;

$onDotL
file fx / McCormic.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, put x.l(j)); put z.l /);