ampl/tp085.mod


If Cut/Copy and Paste fails, then click here for download.


# Model tp085
# Source version 1
# Simple APM canonicalizer version 1.3
# AMPL back end version 1.1
param coefy5a := .004324;
param coefy5b := .00423;
param coefy5 := coefy5a;
param a_2 := 17.505;
param a_3 := 11.275;
param a_4 := 214.228;
param a_5 := 7.458;
param a_6 := .961;
param a_7 := 1.612;
param a_8 := .146;
param a_9 := 107.99;
param a_10 := 922.693;
param a_11 := 926.832;
param a_12 := 18.766;
param a_13 := 1072.163;
param a_14 := 8961.448;
param a_15 := .063;
param a_16 := 71084.33;
param a_17 := 2802713;
param b_2 := 1053.6667;
param b_3 := 35.03;
param b_4 := 665.585;
param b_5 := 584.463;
param b_6 := 265.916;
param b_7 := 7.046;
param b_8 := .222;
param b_9 := 273.366;
param b_10 := 1286.105;
param b_11 := 1444.046;
param b_12 := 537.141;
param b_13 := 3247.039;
param b_14 := 26844.086;
param b_15 := .386;
param b_16 := 140000;
param b_17 := 12146108;
param c_10 := 123/7523;
var x_1 >= 704.4148, <= 906.3855;
let x_1 := 900;
var x_2 >= 68.6, <= 288.88;
let x_2 := 80;
var x_3 >= 0, <= 134.75;
let x_3 := 115;
var x_4 >= 193, <= 287.0966;
let x_4 := 267;
var x_5 >= 25, <= 84.1988;
let x_5 := 27;
var y_1 = x_2 + x_3 + 41.6;
var c_1 = .024*x_4 - 4.62;
var y_2 = 12.5/c_1 + 12;
var c_2 = .0003535*x_1^2 + .5311*x_1 + .08705*y_2*x_1;
var c_3 = .052*x_1 + 78 + .002377*y_2*x_1;
var y_3 = c_2/c_3;
var y_4 = 19*y_3;
var c_4 = .04782*(x_1 - y_3) + .1956*(x_1 - y_3)^2/x_2 + .6376*y_4 + 1.594*y_3;
var c_5 = 100*x_2;
var c_6 = x_1 - y_3 - y_4;
var c_7 = .95 - c_4/c_5;
var y_5 = c_6*c_7;
var y_6 = x_1 - y_5 - y_4 - y_3;
var c_8 = (y_5 + y_4)*.995;
var y_7 = c_8/y_1;
var y_8 = c_8/3798;
var c_9 = y_7 - .0663*y_7/y_8 - .3153;
var y_9 = 96.82/c_9 + .321*y_1;
var y_10 = 1.29*y_5 + 1.258*y_4 + 2.29*y_3 + 1.71*y_6;
var y_11 = 1.71*x_1 - .452*y_4 + .58*y_3;
var c_11 = (1.75*y_2)*(.995*x_1);
var c_12 = .995*y_10 + 1998;
var y_12 = c_10*x_1 + c_11/c_12;
var y_13 = c_12 - 1.75*y_2;
var y_14 = 3623 + 64.4*x_2 + 58.4*x_3 + 146312/(y_9 + x_5);
var c_13 = .995*y_10 + 60.8*x_2 + 48*x_4 - .1121*y_14 - 5095;
var y_15 = y_13/c_13;
var y_16 = 148000 - 331000*y_15 + 40*y_13 - 61*y_15*y_13;
var c_14 = 2324*y_10 - 28740000*y_2;
var y_17 = 14130000 - 1328*y_10 - 531*y_11 + c_14/c_12;
var c_15 = y_13/y_15 - y_13/.52;
var c_16 = 1.104 - .72*y_15;
var c_17 = y_9 + x_5;
var myminfun = (-5.843e-7)*y_17 + 1.17e-4*y_14 + 2.358e-5*y_13 + 1.502e-6*y_16 + .0321*y_12 + coefy5*y_5 +
1.0e-4*c_15/c_16 + 37.48*y_2/c_12 + .1365;
var myabsdevnod0_0 = abs(myminfun + 1.905155258534784);
var myreldevnod0_0 = abs(myminfun/1.905155258534784 + 1);
var myabsdevnod0_1 = abs(x_1 - 705.1745370700908);
var myreldevnod0_1 = abs(x_1/705.1745370700908 - 1);
var myabsdevnod0_2 = abs(x_2 - 68.6);
var myreldevnod0_2 = abs(x_2/68.6 - 1);
var myabsdevnod0_3 = abs(x_3 - 102.9);
var myreldevnod0_3 = abs(x_3/102.9 - 1);
var myabsdevnod0_4 = abs(x_4 - 282.3249315936603);
var myreldevnod0_4 = abs(x_4/282.3249315936603 - 1);
var myabsdevnod0_5 = abs(x_5 - 37.58411642580555);
var myreldevnod0_5 = abs(x_5/37.58411642580555 - 1);
var myabsdevnod1_0 = (1/2)*(abs(myabsdevnod0_0 - myabsdevnod0_1) + (myabsdevnod0_0 + myabsdevnod0_1));
var myreldevnod1_0 = (1/2)*(abs(myreldevnod0_0 - myreldevnod0_1) + (myreldevnod0_0 + myreldevnod0_1));
var myabsdevnod1_2 = (1/2)*(abs(myabsdevnod0_2 - myabsdevnod0_3) + (myabsdevnod0_2 + myabsdevnod0_3));
var myreldevnod1_2 = (1/2)*(abs(myreldevnod0_2 - myreldevnod0_3) + (myreldevnod0_2 + myreldevnod0_3));
var myabsdevnod1_4 = (1/2)*(abs(myabsdevnod0_4 - myabsdevnod0_5) + (myabsdevnod0_4 + myabsdevnod0_5));
var myreldevnod1_4 = (1/2)*(abs(myreldevnod0_4 - myreldevnod0_5) + (myreldevnod0_4 + myreldevnod0_5));
var myabsdevnod2_0 = (1/2)*(abs(myabsdevnod1_0 - myabsdevnod1_2) + (myabsdevnod1_0 + myabsdevnod1_2));
var myreldevnod2_0 = (1/2)*(abs(myreldevnod1_0 - myreldevnod1_2) + (myreldevnod1_0 + myreldevnod1_2));
var myabsdevnod2_4 = myabsdevnod1_4;
var myreldevnod2_4 = myreldevnod1_4;
var myabsdevnod3_0 = (1/2)*(abs(myabsdevnod2_0 - myabsdevnod2_4) + (myabsdevnod2_0 + myabsdevnod2_4));
var myreldevnod3_0 = (1/2)*(abs(myreldevnod2_0 - myreldevnod2_4) + (myreldevnod2_0 + myreldevnod2_4));
var zmyabsdevmax = myabsdevnod3_0;
var zmyreldevmax = myreldevnod3_0;
minimize obj: myminfun;
subject to ci_1: 1.5*x_2 - x_3 >= 0;
subject to ci_2: y_1 - 213.1 >= 0;
subject to ci_3: 405.23 - y_1 >= 0;
subject to ci_4: y_2 - a_2 >= 0;
subject to ci_5: y_3 - a_3 >= 0;
subject to ci_6: y_4 - a_4 >= 0;
subject to ci_7: y_5 - a_5 >= 0;
subject to ci_8: y_6 - a_6 >= 0;
subject to ci_9: y_7 - a_7 >= 0;
subject to ci_10: y_8 - a_8 >= 0;
subject to ci_11: y_9 - a_9 >= 0;
subject to ci_12: y_10 - a_10 >= 0;
subject to ci_13: y_11 - a_11 >= 0;
subject to ci_14: y_12 - a_12 >= 0;
subject to ci_15: y_13 - a_13 >= 0;
subject to ci_16: y_14 - a_14 >= 0;
subject to ci_17: y_15 - a_15 >= 0;
subject to ci_18: y_16 - a_16 >= 0;
subject to ci_19: y_17 - a_17 >= 0;
subject to ci_20: b_2 - y_2 >= 0;
subject to ci_21: b_3 - y_3 >= 0;
subject to ci_22: b_4 - y_4 >= 0;
subject to ci_23: b_5 - y_5 >= 0;
subject to ci_24: b_6 - y_6 >= 0;
subject to ci_25: b_7 - y_7 >= 0;
subject to ci_26: b_8 - y_8 >= 0;
subject to ci_27: b_9 - y_9 >= 0;
subject to ci_28: b_10 - y_10 >= 0;
subject to ci_29: b_11 - y_11 >= 0;
subject to ci_30: b_12 - y_12 >= 0;
subject to ci_31: b_13 - y_13 >= 0;
subject to ci_32: b_14 - y_14 >= 0;
subject to ci_33: b_15 - y_15 >= 0;
subject to ci_34: b_16 - y_16 >= 0;
subject to ci_35: b_17 - y_17 >= 0;
subject to ci_36: y_4 - (28/72)*y_5 >= 0;
subject to ci_37: 21 - 3496*y_2/c_12 >= 0;
subject to ci_38: 62212/c_17 - 110.6 - y_1 >= 0;
solve;
display x_1;
display x_2;
display x_3;
display x_4;
display x_5;
display myminfun;
display zmyabsdevmax;
display zmyreldevmax;
# End Model

Stephan K.H. Seidl