set exper /1*16/
    quantities /Z1, Z2, Z3, dZ1, dZ2, dZ3/;

variables dZ1_mod, dZ2_mod, dZ3_mod, b, P1, P2, P3, P4, P5, w;

table meas(exper,quantities)
         Z1           Z2           Z3           dZ1                 dZ2                dZ3
1        0.461        0.114        0.018        -5.963163056        1.704819365        0.432035893
2        0.426        0.135        0.035        -4.235380505        1.434204277        0.404761617
3        0.383        0.157        0.045        -2.957111609        1.196382402        0.380346428
4        0.305        0.194        0.047        -1.139800037        0.68750181         0.330465563
5        0.195        0.231        0.084        -0.514625236        0.077647376        0.289494016
6        0.17         0.234        0.095        -0.5233168          -0.004498955       0.29155085
7        0.139        0.228        0.111        -0.496843591        -0.056196706       0.291348409
8        0.112        0.228        0.134        -0.401976585        -0.055158335       0.276883213
9        0.112        0.226        0.168        -0.319485734        -0.062592085       0.261990297
10       0.09         0.22         0.148        -0.185254347        -0.097443168       0.228100469
11       0.082        0.214        0.157        -0.143860262        -0.118719804       0.230054365
12       0.066        0.198        0.206        -0.243698403        -0.132191801       0.172029741
13       0.053        0.188        0.206        -0.274762504        -0.095359082       0.103870025
14       0.043        0.183        0.214        0.193317728         0.017820456        -0.000417513
15       0.041        0.184        0.213        0.263109921         0.02082751         -0.002363041
16       0.029        0.166        0.23         -0.94069295         -0.329170417       0.259580405
;



equations E1, E2, E3, obj;

E1(exper)..  dZ1_mod(exper) =E=  -(2*P1-P1*meas(exper, 'Z2')/((P2+P5)*meas(exper, 'Z1')+meas(exper, 'Z2'))+P3+P4)*meas(exper, 'Z1');
E2(exper)..  dZ2_mod(exper) =E=   P1*meas(exper, 'Z1')*(P2*meas(exper, 'Z1')-meas(exper, 'Z2'))/((P2+P5)*meas(exper, 'Z1')+meas(exper, 'Z2'))+P3*meas(exper, 'Z1');
E3(exper)..  dZ3_mod(exper) =E=   P1*meas(exper, 'Z1')*(meas(exper, 'Z2')+P5*meas(exper, 'Z1'))/((P2+P5)*meas(exper, 'Z1')+meas(exper, 'Z2'))+P4*meas(exper, 'Z1');

obj..        w              =E= sum( exper, (dZ1_mod(exper) - meas(exper, 'dZ1'))*(dZ1_mod(exper) - meas(exper, 'dZ1')) + (dZ2_mod(exper) - meas(exper, 'dZ2'))*(dZ2_mod(exper) - meas(exper, 'dZ2')) + (dZ3_mod(exper) - meas(exper, 'dZ3'))*(dZ3_mod(exper) - meas(exper, 'dZ3')) );



P1.lo = 0;
P1.up = 100;
P2.lo = 0;
P2.up = 100;
P3.lo = 0;
P3.up = 100;
P4.lo = 0;
P4.up = 100;
P5.lo = 0;
P5.up = 100;

 MODEL TEST / E1, E2, E3, obj / ;

 OPTION LIMROW = 0.0;
 OPTION LIMCOL = 0.0;
*OPTION MINLP=SBB;
 Test.optfile = 1;

 SOLVE TEST USING NLP minimiZING w;