# Domains
var x1 >= -4, <= 5;
var x2 >= -4, <= 5;
var x3 >= -4, <= 5;
var x4 >= -4, <= 5;
var x5 >= -4, <= 5;
var x6 >= -4, <= 5;
var x7 >= -4, <= 5;
var x8 >= -4, <= 5;
var x9 >= -4, <= 5;
var x10 >= -4, <= 5;
var x11 >= -4, <= 5;
var x12 >= -4, <= 5;
var x13 >= -4, <= 5;
var x14 >= -4, <= 5;
var x15 >= -4, <= 5;
var x16 >= -4, <= 5;
var x17 >= -4, <= 5;
var x18 >= -4, <= 5;
var x19 >= -4, <= 5;
var x20 >= -4, <= 5;
var x21 >= -4, <= 5;
var x22 >= -4, <= 5;
var x23 >= -4, <= 5;
var x24 >= -4, <= 5;
var x25 >= -4, <= 5;
var x26 >= -4, <= 5;
var x27 >= -4, <= 5;
var x28 >= -4, <= 5;
var x29 >= -4, <= 5;
var x30 >= -4, <= 5;
var x31 >= -4, <= 5;
var x32 >= -4, <= 5;
var x33 >= -4, <= 5;
var x34 >= -4, <= 5;
var x35 >= -4, <= 5;
var x36 >= -4, <= 5;
var x37 >= -4, <= 5;
var x38 >= -4, <= 5;
var x39 >= -4, <= 5;
var x40 >= -4, <= 5;
var x41 >= -4, <= 5;
var x42 >= -4, <= 5;
var x43 >= -4, <= 5;
var x44 >= -4, <= 5;
var x45 >= -4, <= 5;
var x46 >= -4, <= 5;
var x47 >= -4, <= 5;
var x48 >= -4, <= 5;
var x49 >= -4, <= 5;
var x50 >= -4, <= 5;
var x51 >= -4, <= 5;
var x52 >= -4, <= 5;
var x53 >= -4, <= 5;
var x54 >= -4, <= 5;
var x55 >= -4, <= 5;
var x56 >= -4, <= 5;
var x57 >= -4, <= 5;
var x58 >= -4, <= 5;
var x59 >= -4, <= 5;
var x60 >= -4, <= 5;
var x61 >= -4, <= 5;
var x62 >= -4, <= 5;
var x63 >= -4, <= 5;
var x64 >= -4, <= 5;
var x65 >= -4, <= 5;
var x66 >= -4, <= 5;
var x67 >= -4, <= 5;
var x68 >= -4, <= 5;
var x69 >= -4, <= 5;
var x70 >= -4, <= 5;
var x71 >= -4, <= 5;
var x72 >= -4, <= 5;
var x73 >= -4, <= 5;
var x74 >= -4, <= 5;
var x75 >= -4, <= 5;
var x76 >= -4, <= 5;
var x77 >= -4, <= 5;
var x78 >= -4, <= 5;
var x79 >= -4, <= 5;
var x80 >= -4, <= 5;
var l0 >= -1.0e16, <= 1.0e16;
var l1 >= -1.0e16, <= 1.0e16;
var l2 >= -1.0e16, <= 1.0e16;
var l3 >= -1.0e16, <= 1.0e16;
var l4 >= -1.0e16, <= 1.0e16;
var l5 >= -1.0e16, <= 1.0e16;
var l6 >= -1.0e16, <= 1.0e16;
var l7 >= -1.0e16, <= 1.0e16;
var l8 >= -1.0e16, <= 1.0e16;
var l9 >= -1.0e16, <= 1.0e16;
var l10 >= -1.0e16, <= 1.0e16;
var l11 >= -1.0e16, <= 1.0e16;
var l12 >= -1.0e16, <= 1.0e16;
var l13 >= -1.0e16, <= 1.0e16;
var l14 >= -1.0e16, <= 1.0e16;
var l15 >= -1.0e16, <= 1.0e16;
var l16 >= -1.0e16, <= 1.0e16;
var l17 >= -1.0e16, <= 1.0e16;
var l18 >= -1.0e16, <= 1.0e16;
var l19 >= -1.0e16, <= 1.0e16;
var l20 >= -1.0e16, <= 1.0e16;
var l21 >= -1.0e16, <= 1.0e16;
var l22 >= -1.0e16, <= 1.0e16;
var l23 >= -1.0e16, <= 1.0e16;
var l24 >= -1.0e16, <= 1.0e16;
var l25 >= -1.0e16, <= 1.0e16;
var l26 >= -1.0e16, <= 1.0e16;
var l27 >= -1.0e16, <= 1.0e16;
var l28 >= -1.0e16, <= 1.0e16;
var l29 >= -1.0e16, <= 1.0e16;
var l30 >= -1.0e16, <= 1.0e16;
var l31 >= -1.0e16, <= 1.0e16;
var l32 >= -1.0e16, <= 1.0e16;
var l33 >= -1.0e16, <= 1.0e16;
var l34 >= -1.0e16, <= 1.0e16;
var l35 >= -1.0e16, <= 1.0e16;
var l36 >= -1.0e16, <= 1.0e16;
var l37 >= -1.0e16, <= 1.0e16;
var l38 >= -1.0e16, <= 1.0e16;
var l39 >= -1.0e16, <= 1.0e16;
var l40 >= -1.0e16, <= 1.0e16;
var l41 >= -1.0e16, <= 1.0e16;
var l42 >= -1.0e16, <= 1.0e16;
var l43 >= -1.0e16, <= 1.0e16;
var l44 >= -1.0e16, <= 1.0e16;
var l45 >= -1.0e16, <= 1.0e16;
var l46 >= -1.0e16, <= 1.0e16;
var l47 >= -1.0e16, <= 1.0e16;
var l48 >= -1.0e16, <= 1.0e16;
var l49 >= -1.0e16, <= 1.0e16;
var l50 >= -1.0e16, <= 1.0e16;
var l51 >= -1.0e16, <= 1.0e16;
var l52 >= -1.0e16, <= 1.0e16;
var l53 >= -1.0e16, <= 1.0e16;
var l54 >= -1.0e16, <= 1.0e16;
var l55 >= -1.0e16, <= 1.0e16;
var l56 >= -1.0e16, <= 1.0e16;
var l57 >= -1.0e16, <= 1.0e16;
var l58 >= -1.0e16, <= 1.0e16;
var l59 >= -1.0e16, <= 1.0e16;
var l60 >= -1.0e16, <= 1.0e16;
var l61 >= -1.0e16, <= 1.0e16;
var l62 >= -1.0e16, <= 1.0e16;
var l63 >= -1.0e16, <= 1.0e16;
var l64 >= -1.0e16, <= 1.0e16;
var l65 >= -1.0e16, <= 1.0e16;
var l66 >= -1.0e16, <= 1.0e16;
var l67 >= -1.0e16, <= 1.0e16;
var l68 >= -1.0e16, <= 1.0e16;
var l69 >= -1.0e16, <= 1.0e16;
var l70 >= -1.0e16, <= 1.0e16;
var l71 >= -1.0e16, <= 1.0e16;
var l72 >= -1.0e16, <= 1.0e16;
var l73 >= -1.0e16, <= 1.0e16;
var l74 >= -1.0e16, <= 1.0e16;
var l75 >= -1.0e16, <= 1.0e16;
var l76 >= -1.0e16, <= 1.0e16;
var l77 >= -1.0e16, <= 1.0e16;
var l78 >= -1.0e16, <= 1.0e16;
var l79 >= -1.0e16, <= 1.0e16;
var l80 >= -1.0e16, <= 1.0e16;
var u1 >= -1.0e16, <= 1.0e16;
var u2 >= -1.0e16, <= 1.0e16;
var u3 >= -1.0e16, <= 1.0e16;
var u4 >= -1.0e16, <= 1.0e16;
var u5 >= -1.0e16, <= 1.0e16;
var u6 >= -1.0e16, <= 1.0e16;
var u7 >= -1.0e16, <= 1.0e16;
var u8 >= -1.0e16, <= 1.0e16;
var u9 >= -1.0e16, <= 1.0e16;
var u10 >= -1.0e16, <= 1.0e16;
var u11 >= -1.0e16, <= 1.0e16;
var u12 >= -1.0e16, <= 1.0e16;
var u13 >= -1.0e16, <= 1.0e16;
var u14 >= -1.0e16, <= 1.0e16;
var u15 >= -1.0e16, <= 1.0e16;
var u16 >= -1.0e16, <= 1.0e16;
var u17 >= -1.0e16, <= 1.0e16;
var u18 >= -1.0e16, <= 1.0e16;
var u19 >= -1.0e16, <= 1.0e16;
var u20 >= -1.0e16, <= 1.0e16;
var u21 >= -1.0e16, <= 1.0e16;
var u22 >= -1.0e16, <= 1.0e16;
var u23 >= -1.0e16, <= 1.0e16;
var u24 >= -1.0e16, <= 1.0e16;
var u25 >= -1.0e16, <= 1.0e16;
var u26 >= -1.0e16, <= 1.0e16;
var u27 >= -1.0e16, <= 1.0e16;
var u28 >= -1.0e16, <= 1.0e16;
var u29 >= -1.0e16, <= 1.0e16;
var u30 >= -1.0e16, <= 1.0e16;
var u31 >= -1.0e16, <= 1.0e16;
var u32 >= -1.0e16, <= 1.0e16;
var u33 >= -1.0e16, <= 1.0e16;
var u34 >= -1.0e16, <= 1.0e16;
var u35 >= -1.0e16, <= 1.0e16;
var u36 >= -1.0e16, <= 1.0e16;
var u37 >= -1.0e16, <= 1.0e16;
var u38 >= -1.0e16, <= 1.0e16;
var u39 >= -1.0e16, <= 1.0e16;
var u40 >= -1.0e16, <= 1.0e16;
var u41 >= -1.0e16, <= 1.0e16;
var u42 >= -1.0e16, <= 1.0e16;
var u43 >= -1.0e16, <= 1.0e16;
var u44 >= -1.0e16, <= 1.0e16;
var u45 >= -1.0e16, <= 1.0e16;
var u46 >= -1.0e16, <= 1.0e16;
var u47 >= -1.0e16, <= 1.0e16;
var u48 >= -1.0e16, <= 1.0e16;
var u49 >= -1.0e16, <= 1.0e16;
var u50 >= -1.0e16, <= 1.0e16;
var u51 >= -1.0e16, <= 1.0e16;
var u52 >= -1.0e16, <= 1.0e16;
var u53 >= -1.0e16, <= 1.0e16;
var u54 >= -1.0e16, <= 1.0e16;
var u55 >= -1.0e16, <= 1.0e16;
var u56 >= -1.0e16, <= 1.0e16;
var u57 >= -1.0e16, <= 1.0e16;
var u58 >= -1.0e16, <= 1.0e16;
var u59 >= -1.0e16, <= 1.0e16;
var u60 >= -1.0e16, <= 1.0e16;
var u61 >= -1.0e16, <= 1.0e16;
var u62 >= -1.0e16, <= 1.0e16;
var u63 >= -1.0e16, <= 1.0e16;
var u64 >= -1.0e16, <= 1.0e16;
var u65 >= -1.0e16, <= 1.0e16;
var u66 >= -1.0e16, <= 1.0e16;
var u67 >= -1.0e16, <= 1.0e16;
var u68 >= -1.0e16, <= 1.0e16;
var u69 >= -1.0e16, <= 1.0e16;
var u70 >= -1.0e16, <= 1.0e16;
var u71 >= -1.0e16, <= 1.0e16;
var u72 >= -1.0e16, <= 1.0e16;
var u73 >= -1.0e16, <= 1.0e16;
var u74 >= -1.0e16, <= 1.0e16;
var u75 >= -1.0e16, <= 1.0e16;
var u76 >= -1.0e16, <= 1.0e16;
var u77 >= -1.0e16, <= 1.0e16;
var u78 >= -1.0e16, <= 1.0e16;
var u79 >= -1.0e16, <= 1.0e16;
var u80 >= -1.0e16, <= 1.0e16;
var u81 >= -1.0e16, <= 1.0e16;

# Constants
param h := 1/(80+1);
param t1 := 1 * h;
param t2 := 2 * h;
param t3 := 3 * h;
param t4 := 4 * h;
param t5 := 5 * h;
param t6 := 6 * h;
param t7 := 7 * h;
param t8 := 8 * h;
param t9 := 9 * h;
param t10 := 10 * h;
param t11 := 11 * h;
param t12 := 12 * h;
param t13 := 13 * h;
param t14 := 14 * h;
param t15 := 15 * h;
param t16 := 16 * h;
param t17 := 17 * h;
param t18 := 18 * h;
param t19 := 19 * h;
param t20 := 20 * h;
param t21 := 21 * h;
param t22 := 22 * h;
param t23 := 23 * h;
param t24 := 24 * h;
param t25 := 25 * h;
param t26 := 26 * h;
param t27 := 27 * h;
param t28 := 28 * h;
param t29 := 29 * h;
param t30 := 30 * h;
param t31 := 31 * h;
param t32 := 32 * h;
param t33 := 33 * h;
param t34 := 34 * h;
param t35 := 35 * h;
param t36 := 36 * h;
param t37 := 37 * h;
param t38 := 38 * h;
param t39 := 39 * h;
param t40 := 40 * h;
param t41 := 41 * h;
param t42 := 42 * h;
param t43 := 43 * h;
param t44 := 44 * h;
param t45 := 45 * h;
param t46 := 46 * h;
param t47 := 47 * h;
param t48 := 48 * h;
param t49 := 49 * h;
param t50 := 50 * h;
param t51 := 51 * h;
param t52 := 52 * h;
param t53 := 53 * h;
param t54 := 54 * h;
param t55 := 55 * h;
param t56 := 56 * h;
param t57 := 57 * h;
param t58 := 58 * h;
param t59 := 59 * h;
param t60 := 60 * h;
param t61 := 61 * h;
param t62 := 62 * h;
param t63 := 63 * h;
param t64 := 64 * h;
param t65 := 65 * h;
param t66 := 66 * h;
param t67 := 67 * h;
param t68 := 68 * h;
param t69 := 69 * h;
param t70 := 70 * h;
param t71 := 71 * h;
param t72 := 72 * h;
param t73 := 73 * h;
param t74 := 74 * h;
param t75 := 75 * h;
param t76 := 76 * h;
param t77 := 77 * h;
param t78 := 78 * h;
param t79 := 79 * h;
param t80 := 80 * h;

subject to
cons1 : l0 = 0;
cons2 : l1 = l0 + t1*(x1+t1+1)^3;
cons3 : l2 = l1 + t2*(x2+t2+1)^3;
cons4 : l3 = l2 + t3*(x3+t3+1)^3;
cons5 : l4 = l3 + t4*(x4+t4+1)^3;
cons6 : l5 = l4 + t5*(x5+t5+1)^3;
cons7 : l6 = l5 + t6*(x6+t6+1)^3;
cons8 : l7 = l6 + t7*(x7+t7+1)^3;
cons9 : l8 = l7 + t8*(x8+t8+1)^3;
cons10 : l9 = l8 + t9*(x9+t9+1)^3;
cons11 : l10 = l9 + t10*(x10+t10+1)^3;
cons12 : l11 = l10 + t11*(x11+t11+1)^3;
cons13 : l12 = l11 + t12*(x12+t12+1)^3;
cons14 : l13 = l12 + t13*(x13+t13+1)^3;
cons15 : l14 = l13 + t14*(x14+t14+1)^3;
cons16 : l15 = l14 + t15*(x15+t15+1)^3;
cons17 : l16 = l15 + t16*(x16+t16+1)^3;
cons18 : l17 = l16 + t17*(x17+t17+1)^3;
cons19 : l18 = l17 + t18*(x18+t18+1)^3;
cons20 : l19 = l18 + t19*(x19+t19+1)^3;
cons21 : l20 = l19 + t20*(x20+t20+1)^3;
cons22 : l21 = l20 + t21*(x21+t21+1)^3;
cons23 : l22 = l21 + t22*(x22+t22+1)^3;
cons24 : l23 = l22 + t23*(x23+t23+1)^3;
cons25 : l24 = l23 + t24*(x24+t24+1)^3;
cons26 : l25 = l24 + t25*(x25+t25+1)^3;
cons27 : l26 = l25 + t26*(x26+t26+1)^3;
cons28 : l27 = l26 + t27*(x27+t27+1)^3;
cons29 : l28 = l27 + t28*(x28+t28+1)^3;
cons30 : l29 = l28 + t29*(x29+t29+1)^3;
cons31 : l30 = l29 + t30*(x30+t30+1)^3;
cons32 : l31 = l30 + t31*(x31+t31+1)^3;
cons33 : l32 = l31 + t32*(x32+t32+1)^3;
cons34 : l33 = l32 + t33*(x33+t33+1)^3;
cons35 : l34 = l33 + t34*(x34+t34+1)^3;
cons36 : l35 = l34 + t35*(x35+t35+1)^3;
cons37 : l36 = l35 + t36*(x36+t36+1)^3;
cons38 : l37 = l36 + t37*(x37+t37+1)^3;
cons39 : l38 = l37 + t38*(x38+t38+1)^3;
cons40 : l39 = l38 + t39*(x39+t39+1)^3;
cons41 : l40 = l39 + t40*(x40+t40+1)^3;
cons42 : l41 = l40 + t41*(x41+t41+1)^3;
cons43 : l42 = l41 + t42*(x42+t42+1)^3;
cons44 : l43 = l42 + t43*(x43+t43+1)^3;
cons45 : l44 = l43 + t44*(x44+t44+1)^3;
cons46 : l45 = l44 + t45*(x45+t45+1)^3;
cons47 : l46 = l45 + t46*(x46+t46+1)^3;
cons48 : l47 = l46 + t47*(x47+t47+1)^3;
cons49 : l48 = l47 + t48*(x48+t48+1)^3;
cons50 : l49 = l48 + t49*(x49+t49+1)^3;
cons51 : l50 = l49 + t50*(x50+t50+1)^3;
cons52 : l51 = l50 + t51*(x51+t51+1)^3;
cons53 : l52 = l51 + t52*(x52+t52+1)^3;
cons54 : l53 = l52 + t53*(x53+t53+1)^3;
cons55 : l54 = l53 + t54*(x54+t54+1)^3;
cons56 : l55 = l54 + t55*(x55+t55+1)^3;
cons57 : l56 = l55 + t56*(x56+t56+1)^3;
cons58 : l57 = l56 + t57*(x57+t57+1)^3;
cons59 : l58 = l57 + t58*(x58+t58+1)^3;
cons60 : l59 = l58 + t59*(x59+t59+1)^3;
cons61 : l60 = l59 + t60*(x60+t60+1)^3;
cons62 : l61 = l60 + t61*(x61+t61+1)^3;
cons63 : l62 = l61 + t62*(x62+t62+1)^3;
cons64 : l63 = l62 + t63*(x63+t63+1)^3;
cons65 : l64 = l63 + t64*(x64+t64+1)^3;
cons66 : l65 = l64 + t65*(x65+t65+1)^3;
cons67 : l66 = l65 + t66*(x66+t66+1)^3;
cons68 : l67 = l66 + t67*(x67+t67+1)^3;
cons69 : l68 = l67 + t68*(x68+t68+1)^3;
cons70 : l69 = l68 + t69*(x69+t69+1)^3;
cons71 : l70 = l69 + t70*(x70+t70+1)^3;
cons72 : l71 = l70 + t71*(x71+t71+1)^3;
cons73 : l72 = l71 + t72*(x72+t72+1)^3;
cons74 : l73 = l72 + t73*(x73+t73+1)^3;
cons75 : l74 = l73 + t74*(x74+t74+1)^3;
cons76 : l75 = l74 + t75*(x75+t75+1)^3;
cons77 : l76 = l75 + t76*(x76+t76+1)^3;
cons78 : l77 = l76 + t77*(x77+t77+1)^3;
cons79 : l78 = l77 + t78*(x78+t78+1)^3;
cons80 : l79 = l78 + t79*(x79+t79+1)^3;
cons81 : l80 = l79 + t80*(x80+t80+1)^3;
cons82 : u81 = 0;
cons83 : u1 = u2 + (1-t1)*(x1+t1+1)^3;
cons84 : u2 = u3 + (1-t2)*(x2+t2+1)^3;
cons85 : u3 = u4 + (1-t3)*(x3+t3+1)^3;
cons86 : u4 = u5 + (1-t4)*(x4+t4+1)^3;
cons87 : u5 = u6 + (1-t5)*(x5+t5+1)^3;
cons88 : u6 = u7 + (1-t6)*(x6+t6+1)^3;
cons89 : u7 = u8 + (1-t7)*(x7+t7+1)^3;
cons90 : u8 = u9 + (1-t8)*(x8+t8+1)^3;
cons91 : u9 = u10 + (1-t9)*(x9+t9+1)^3;
cons92 : u10 = u11 + (1-t10)*(x10+t10+1)^3;
cons93 : u11 = u12 + (1-t11)*(x11+t11+1)^3;
cons94 : u12 = u13 + (1-t12)*(x12+t12+1)^3;
cons95 : u13 = u14 + (1-t13)*(x13+t13+1)^3;
cons96 : u14 = u15 + (1-t14)*(x14+t14+1)^3;
cons97 : u15 = u16 + (1-t15)*(x15+t15+1)^3;
cons98 : u16 = u17 + (1-t16)*(x16+t16+1)^3;
cons99 : u17 = u18 + (1-t17)*(x17+t17+1)^3;
cons100 : u18 = u19 + (1-t18)*(x18+t18+1)^3;
cons101 : u19 = u20 + (1-t19)*(x19+t19+1)^3;
cons102 : u20 = u21 + (1-t20)*(x20+t20+1)^3;
cons103 : u21 = u22 + (1-t21)*(x21+t21+1)^3;
cons104 : u22 = u23 + (1-t22)*(x22+t22+1)^3;
cons105 : u23 = u24 + (1-t23)*(x23+t23+1)^3;
cons106 : u24 = u25 + (1-t24)*(x24+t24+1)^3;
cons107 : u25 = u26 + (1-t25)*(x25+t25+1)^3;
cons108 : u26 = u27 + (1-t26)*(x26+t26+1)^3;
cons109 : u27 = u28 + (1-t27)*(x27+t27+1)^3;
cons110 : u28 = u29 + (1-t28)*(x28+t28+1)^3;
cons111 : u29 = u30 + (1-t29)*(x29+t29+1)^3;
cons112 : u30 = u31 + (1-t30)*(x30+t30+1)^3;
cons113 : u31 = u32 + (1-t31)*(x31+t31+1)^3;
cons114 : u32 = u33 + (1-t32)*(x32+t32+1)^3;
cons115 : u33 = u34 + (1-t33)*(x33+t33+1)^3;
cons116 : u34 = u35 + (1-t34)*(x34+t34+1)^3;
cons117 : u35 = u36 + (1-t35)*(x35+t35+1)^3;
cons118 : u36 = u37 + (1-t36)*(x36+t36+1)^3;
cons119 : u37 = u38 + (1-t37)*(x37+t37+1)^3;
cons120 : u38 = u39 + (1-t38)*(x38+t38+1)^3;
cons121 : u39 = u40 + (1-t39)*(x39+t39+1)^3;
cons122 : u40 = u41 + (1-t40)*(x40+t40+1)^3;
cons123 : u41 = u42 + (1-t41)*(x41+t41+1)^3;
cons124 : u42 = u43 + (1-t42)*(x42+t42+1)^3;
cons125 : u43 = u44 + (1-t43)*(x43+t43+1)^3;
cons126 : u44 = u45 + (1-t44)*(x44+t44+1)^3;
cons127 : u45 = u46 + (1-t45)*(x45+t45+1)^3;
cons128 : u46 = u47 + (1-t46)*(x46+t46+1)^3;
cons129 : u47 = u48 + (1-t47)*(x47+t47+1)^3;
cons130 : u48 = u49 + (1-t48)*(x48+t48+1)^3;
cons131 : u49 = u50 + (1-t49)*(x49+t49+1)^3;
cons132 : u50 = u51 + (1-t50)*(x50+t50+1)^3;
cons133 : u51 = u52 + (1-t51)*(x51+t51+1)^3;
cons134 : u52 = u53 + (1-t52)*(x52+t52+1)^3;
cons135 : u53 = u54 + (1-t53)*(x53+t53+1)^3;
cons136 : u54 = u55 + (1-t54)*(x54+t54+1)^3;
cons137 : u55 = u56 + (1-t55)*(x55+t55+1)^3;
cons138 : u56 = u57 + (1-t56)*(x56+t56+1)^3;
cons139 : u57 = u58 + (1-t57)*(x57+t57+1)^3;
cons140 : u58 = u59 + (1-t58)*(x58+t58+1)^3;
cons141 : u59 = u60 + (1-t59)*(x59+t59+1)^3;
cons142 : u60 = u61 + (1-t60)*(x60+t60+1)^3;
cons143 : u61 = u62 + (1-t61)*(x61+t61+1)^3;
cons144 : u62 = u63 + (1-t62)*(x62+t62+1)^3;
cons145 : u63 = u64 + (1-t63)*(x63+t63+1)^3;
cons146 : u64 = u65 + (1-t64)*(x64+t64+1)^3;
cons147 : u65 = u66 + (1-t65)*(x65+t65+1)^3;
cons148 : u66 = u67 + (1-t66)*(x66+t66+1)^3;
cons149 : u67 = u68 + (1-t67)*(x67+t67+1)^3;
cons150 : u68 = u69 + (1-t68)*(x68+t68+1)^3;
cons151 : u69 = u70 + (1-t69)*(x69+t69+1)^3;
cons152 : u70 = u71 + (1-t70)*(x70+t70+1)^3;
cons153 : u71 = u72 + (1-t71)*(x71+t71+1)^3;
cons154 : u72 = u73 + (1-t72)*(x72+t72+1)^3;
cons155 : u73 = u74 + (1-t73)*(x73+t73+1)^3;
cons156 : u74 = u75 + (1-t74)*(x74+t74+1)^3;
cons157 : u75 = u76 + (1-t75)*(x75+t75+1)^3;
cons158 : u76 = u77 + (1-t76)*(x76+t76+1)^3;
cons159 : u77 = u78 + (1-t77)*(x77+t77+1)^3;
cons160 : u78 = u79 + (1-t78)*(x78+t78+1)^3;
cons161 : u79 = u80 + (1-t79)*(x79+t79+1)^3;
cons162 : u80 = u81 + (1-t80)*(x80+t80+1)^3;
cons163 : 0 = x1+0.5*h*((1-t1)*l1 + t1*u2);
cons164 : 0 = x2+0.5*h*((1-t2)*l2 + t2*u3);
cons165 : 0 = x3+0.5*h*((1-t3)*l3 + t3*u4);
cons166 : 0 = x4+0.5*h*((1-t4)*l4 + t4*u5);
cons167 : 0 = x5+0.5*h*((1-t5)*l5 + t5*u6);
cons168 : 0 = x6+0.5*h*((1-t6)*l6 + t6*u7);
cons169 : 0 = x7+0.5*h*((1-t7)*l7 + t7*u8);
cons170 : 0 = x8+0.5*h*((1-t8)*l8 + t8*u9);
cons171 : 0 = x9+0.5*h*((1-t9)*l9 + t9*u10);
cons172 : 0 = x10+0.5*h*((1-t10)*l10 + t10*u11);
cons173 : 0 = x11+0.5*h*((1-t11)*l11 + t11*u12);
cons174 : 0 = x12+0.5*h*((1-t12)*l12 + t12*u13);
cons175 : 0 = x13+0.5*h*((1-t13)*l13 + t13*u14);
cons176 : 0 = x14+0.5*h*((1-t14)*l14 + t14*u15);
cons177 : 0 = x15+0.5*h*((1-t15)*l15 + t15*u16);
cons178 : 0 = x16+0.5*h*((1-t16)*l16 + t16*u17);
cons179 : 0 = x17+0.5*h*((1-t17)*l17 + t17*u18);
cons180 : 0 = x18+0.5*h*((1-t18)*l18 + t18*u19);
cons181 : 0 = x19+0.5*h*((1-t19)*l19 + t19*u20);
cons182 : 0 = x20+0.5*h*((1-t20)*l20 + t20*u21);
cons183 : 0 = x21+0.5*h*((1-t21)*l21 + t21*u22);
cons184 : 0 = x22+0.5*h*((1-t22)*l22 + t22*u23);
cons185 : 0 = x23+0.5*h*((1-t23)*l23 + t23*u24);
cons186 : 0 = x24+0.5*h*((1-t24)*l24 + t24*u25);
cons187 : 0 = x25+0.5*h*((1-t25)*l25 + t25*u26);
cons188 : 0 = x26+0.5*h*((1-t26)*l26 + t26*u27);
cons189 : 0 = x27+0.5*h*((1-t27)*l27 + t27*u28);
cons190 : 0 = x28+0.5*h*((1-t28)*l28 + t28*u29);
cons191 : 0 = x29+0.5*h*((1-t29)*l29 + t29*u30);
cons192 : 0 = x30+0.5*h*((1-t30)*l30 + t30*u31);
cons193 : 0 = x31+0.5*h*((1-t31)*l31 + t31*u32);
cons194 : 0 = x32+0.5*h*((1-t32)*l32 + t32*u33);
cons195 : 0 = x33+0.5*h*((1-t33)*l33 + t33*u34);
cons196 : 0 = x34+0.5*h*((1-t34)*l34 + t34*u35);
cons197 : 0 = x35+0.5*h*((1-t35)*l35 + t35*u36);
cons198 : 0 = x36+0.5*h*((1-t36)*l36 + t36*u37);
cons199 : 0 = x37+0.5*h*((1-t37)*l37 + t37*u38);
cons200 : 0 = x38+0.5*h*((1-t38)*l38 + t38*u39);
cons201 : 0 = x39+0.5*h*((1-t39)*l39 + t39*u40);
cons202 : 0 = x40+0.5*h*((1-t40)*l40 + t40*u41);
cons203 : 0 = x41+0.5*h*((1-t41)*l41 + t41*u42);
cons204 : 0 = x42+0.5*h*((1-t42)*l42 + t42*u43);
cons205 : 0 = x43+0.5*h*((1-t43)*l43 + t43*u44);
cons206 : 0 = x44+0.5*h*((1-t44)*l44 + t44*u45);
cons207 : 0 = x45+0.5*h*((1-t45)*l45 + t45*u46);
cons208 : 0 = x46+0.5*h*((1-t46)*l46 + t46*u47);
cons209 : 0 = x47+0.5*h*((1-t47)*l47 + t47*u48);
cons210 : 0 = x48+0.5*h*((1-t48)*l48 + t48*u49);
cons211 : 0 = x49+0.5*h*((1-t49)*l49 + t49*u50);
cons212 : 0 = x50+0.5*h*((1-t50)*l50 + t50*u51);
cons213 : 0 = x51+0.5*h*((1-t51)*l51 + t51*u52);
cons214 : 0 = x52+0.5*h*((1-t52)*l52 + t52*u53);
cons215 : 0 = x53+0.5*h*((1-t53)*l53 + t53*u54);
cons216 : 0 = x54+0.5*h*((1-t54)*l54 + t54*u55);
cons217 : 0 = x55+0.5*h*((1-t55)*l55 + t55*u56);
cons218 : 0 = x56+0.5*h*((1-t56)*l56 + t56*u57);
cons219 : 0 = x57+0.5*h*((1-t57)*l57 + t57*u58);
cons220 : 0 = x58+0.5*h*((1-t58)*l58 + t58*u59);
cons221 : 0 = x59+0.5*h*((1-t59)*l59 + t59*u60);
cons222 : 0 = x60+0.5*h*((1-t60)*l60 + t60*u61);
cons223 : 0 = x61+0.5*h*((1-t61)*l61 + t61*u62);
cons224 : 0 = x62+0.5*h*((1-t62)*l62 + t62*u63);
cons225 : 0 = x63+0.5*h*((1-t63)*l63 + t63*u64);
cons226 : 0 = x64+0.5*h*((1-t64)*l64 + t64*u65);
cons227 : 0 = x65+0.5*h*((1-t65)*l65 + t65*u66);
cons228 : 0 = x66+0.5*h*((1-t66)*l66 + t66*u67);
cons229 : 0 = x67+0.5*h*((1-t67)*l67 + t67*u68);
cons230 : 0 = x68+0.5*h*((1-t68)*l68 + t68*u69);
cons231 : 0 = x69+0.5*h*((1-t69)*l69 + t69*u70);
cons232 : 0 = x70+0.5*h*((1-t70)*l70 + t70*u71);
cons233 : 0 = x71+0.5*h*((1-t71)*l71 + t71*u72);
cons234 : 0 = x72+0.5*h*((1-t72)*l72 + t72*u73);
cons235 : 0 = x73+0.5*h*((1-t73)*l73 + t73*u74);
cons236 : 0 = x74+0.5*h*((1-t74)*l74 + t74*u75);
cons237 : 0 = x75+0.5*h*((1-t75)*l75 + t75*u76);
cons238 : 0 = x76+0.5*h*((1-t76)*l76 + t76*u77);
cons239 : 0 = x77+0.5*h*((1-t77)*l77 + t77*u78);
cons240 : 0 = x78+0.5*h*((1-t78)*l78 + t78*u79);
cons241 : 0 = x79+0.5*h*((1-t79)*l79 + t79*u80);
cons242 : 0 = x80+0.5*h*((1-t80)*l80 + t80*u81);

solve;
display l0, l1, l10, l11, l12, l13, l14, l15, l16, l17, l18, l19, l2, l20, l21, l22, l23, l24, l25, l26, l27, l28, l29, l3, l30, l31, l32, l33, l34, l35, l36, l37, l38, l39, l4, l40, l41, l42, l43, l44, l45, l46, l47, l48, l49, l5, l50, l51, l52, l53, l54, l55, l56, l57, l58, l59, l6, l60, l61, l62, l63, l64, l65, l66, l67, l68, l69, l7, l70, l71, l72, l73, l74, l75, l76, l77, l78, l79, l8, l80, l9, u1, u10, u11, u12, u13, u14, u15, u16, u17, u18, u19, u2, u20, u21, u22, u23, u24, u25, u26, u27, u28, u29, u3, u30, u31, u32, u33, u34, u35, u36, u37, u38, u39, u4, u40, u41, u42, u43, u44, u45, u46, u47, u48, u49, u5, u50, u51, u52, u53, u54, u55, u56, u57, u58, u59, u6, u60, u61, u62, u63, u64, u65, u66, u67, u68, u69, u7, u70, u71, u72, u73, u74, u75, u76, u77, u78, u79, u8, u80, u81, u9, x1, x10, x11, x12, x13, x14, x15, x16, x17, x18, x19, x2, x20, x21, x22, x23, x24, x25, x26, x27, x28, x29, x3, x30, x31, x32, x33, x34, x35, x36, x37, x38, x39, x4, x40, x41, x42, x43, x44, x45, x46, x47, x48, x49, x5, x50, x51, x52, x53, x54, x55, x56, x57, x58, x59, x6, x60, x61, x62, x63, x64, x65, x66, x67, x68, x69, x7, x70, x71, x72, x73, x74, x75, x76, x77, x78, x79, x8, x80, x9;