*--------------------------------------------------------------*
* Univariate Polynomial - Test Problem 9                       *
*--------------------------------------------------------------*
File  
    res  / results / ;
Put res ;

Sets
    i    /1*2/;

Variables
    x(i)
    f;

x.lo('1') = 0;
x.lo('2') = 0;
x.up('1') = 3;
x.up('2') = 4;

Equations
    Obj    objective function
    Con1   constraint function 1
    Con2   constraint function 2;

Obj .. 
    f =e= - x('1') -x('2');

Con1 ..
    x('2') - 2 - 2*power(x('1'),4) + 8*power(x('1'),3) - 8*power(x('1'),2) =l= 0;

Con2 ..
    x('2') - 4*power(x('1'),4) + 32*power(x('1'),3) - 88*power(x('1'),2) + 96*x('1') - 36 =l= 0;

Model
    problem /Obj, Con1, Con2/;

     x.l('1') = 2.3295;
     x.l('2') = 3.17846;

    solve problem using nlp minimizing f;
    PUT "Min f",f.l:16:10//;
    Loop(i, PUT "x    ",x.l(i):16:10//);