# Domains
var x1 >= -1.0e8, <= 1.0e8;
var x2 >= -1.0e8, <= 1.0e8;
var x3 >= -1.0e8, <= 1.0e8;
var x4 >= -1.0e8, <= 1.0e8;
var x5 >= -1.0e8, <= 1.0e8;
var x6 >= -1.0e8, <= 1.0e8;
var x7 >= -1.0e8, <= 1.0e8;
var x8 >= -1.0e8, <= 1.0e8;
var x9 >= -1.0e8, <= 1.0e8;
var x10 >= -1.0e8, <= 1.0e8;
var x11 >= -1.0e8, <= 1.0e8;
var x12 >= -1.0e8, <= 1.0e8;
var x13 >= -1.0e8, <= 1.0e8;
var x14 >= -1.0e8, <= 1.0e8;
var x15 >= -1.0e8, <= 1.0e8;
var x16 >= -1.0e8, <= 1.0e8;
var x17 >= -1.0e8, <= 1.0e8;
var x18 >= -1.0e8, <= 1.0e8;
var x19 >= -1.0e8, <= 1.0e8;
var x20 >= -1.0e8, <= 1.0e8;
var x21 >= -1.0e8, <= 1.0e8;
var x22 >= -1.0e8, <= 1.0e8;
var x23 >= -1.0e8, <= 1.0e8;
var x24 >= -1.0e8, <= 1.0e8;
var x25 >= -1.0e8, <= 1.0e8;
var x26 >= -1.0e8, <= 1.0e8;
var x27 >= -1.0e8, <= 1.0e8;
var x28 >= -1.0e8, <= 1.0e8;
var x29 >= -1.0e8, <= 1.0e8;
var x30 >= -1.0e8, <= 1.0e8;
var x31 >= -1.0e8, <= 1.0e8;
var x32 >= -1.0e8, <= 1.0e8;
var x33 >= -1.0e8, <= 1.0e8;
var x34 >= -1.0e8, <= 1.0e8;
var x35 >= -1.0e8, <= 1.0e8;
var x36 >= -1.0e8, <= 1.0e8;
var x37 >= -1.0e8, <= 1.0e8;
var x38 >= -1.0e8, <= 1.0e8;
var x39 >= -1.0e8, <= 1.0e8;
var x40 >= -1.0e8, <= 1.0e8;


subject to
cons1 : x1*(2+5*x1^2) + 1 - (x1*(1+x1) + x2*(1+x2)) = 0;
cons2 : x2*(2+5*x2^2) + 1 - (x1*(1+x1) + x3*(1+x3)) = 0;
cons3 : x3*(2+5*x3^2) + 1 - (x1*(1+x1) + x2*(1+x2) + x4*(1+x4)) = 0;
cons4 : x4*(2+5*x4^2) + 1 - (x1*(1+x1) + x2*(1+x2) + x3*(1+x3) + x5*(1+x5)) = 0;
cons5 : x5*(2+5*x5^2) + 1 - (x1*(1+x1) + x2*(1+x2) + x3*(1+x3) + x4*(1+x4) + x6*(1+x6)) = 0;
cons6 : x6*(2+5*x6^2) + 1 - (x1*(1+x1) + x2*(1+x2) + x3*(1+x3) + x4*(1+x4) + x5*(1+x5) + x7*(1+x7)) = 0;
cons7 : x7*(2+5*x7^2) + 1 - (x2*(1+x2) + x3*(1+x3) + x4*(1+x4) + x5*(1+x5) + x6*(1+x6) + x8*(1+x8)) = 0;
cons8 : x8*(2+5*x8^2) + 1 - (x3*(1+x3) + x4*(1+x4) + x5*(1+x5) + x6*(1+x6) + x7*(1+x7) + x9*(1+x9)) = 0;
cons9 : x9*(2+5*x9^2) + 1 - (x4*(1+x4) + x5*(1+x5) + x6*(1+x6) + x7*(1+x7) + x8*(1+x8) + x10*(1+x10)) = 0;
cons10 : x10*(2+5*x10^2) + 1 - (x5*(1+x5) + x6*(1+x6) + x7*(1+x7) + x8*(1+x8) + x9*(1+x9) + x11*(1+x11)) = 0;
cons11 : x11*(2+5*x11^2) + 1 - (x6*(1+x6) + x7*(1+x7) + x8*(1+x8) + x9*(1+x9) + x10*(1+x10) + x12*(1+x12)) = 0;
cons12 : x12*(2+5*x12^2) + 1 - (x7*(1+x7) + x8*(1+x8) + x9*(1+x9) + x10*(1+x10) + x11*(1+x11) + x13*(1+x13)) = 0;
cons13 : x13*(2+5*x13^2) + 1 - (x8*(1+x8) + x9*(1+x9) + x10*(1+x10) + x11*(1+x11) + x12*(1+x12) + x14*(1+x14)) = 0;
cons14 : x14*(2+5*x14^2) + 1 - (x9*(1+x9) + x10*(1+x10) + x11*(1+x11) + x12*(1+x12) + x13*(1+x13) + x15*(1+x15)) = 0;
cons15 : x15*(2+5*x15^2) + 1 - (x10*(1+x10) + x11*(1+x11) + x12*(1+x12) + x13*(1+x13) + x14*(1+x14) + x16*(1+x16)) = 0;
cons16 : x16*(2+5*x16^2) + 1 - (x11*(1+x11) + x12*(1+x12) + x13*(1+x13) + x14*(1+x14) + x15*(1+x15) + x17*(1+x17)) = 0;
cons17 : x17*(2+5*x17^2) + 1 - (x12*(1+x12) + x13*(1+x13) + x14*(1+x14) + x15*(1+x15) + x16*(1+x16) + x18*(1+x18)) = 0;
cons18 : x18*(2+5*x18^2) + 1 - (x13*(1+x13) + x14*(1+x14) + x15*(1+x15) + x16*(1+x16) + x17*(1+x17) + x19*(1+x19)) = 0;
cons19 : x19*(2+5*x19^2) + 1 - (x14*(1+x14) + x15*(1+x15) + x16*(1+x16) + x17*(1+x17) + x18*(1+x18) + x20*(1+x20)) = 0;
cons20 : x20*(2+5*x20^2) + 1 - (x15*(1+x15) + x16*(1+x16) + x17*(1+x17) + x18*(1+x18) + x19*(1+x19) + x21*(1+x21)) = 0;
cons21 : x21*(2+5*x21^2) + 1 - (x16*(1+x16) + x17*(1+x17) + x18*(1+x18) + x19*(1+x19) + x20*(1+x20) + x22*(1+x22)) = 0;
cons22 : x22*(2+5*x22^2) + 1 - (x17*(1+x17) + x18*(1+x18) + x19*(1+x19) + x20*(1+x20) + x21*(1+x21) + x23*(1+x23)) = 0;
cons23 : x23*(2+5*x23^2) + 1 - (x18*(1+x18) + x19*(1+x19) + x20*(1+x20) + x21*(1+x21) + x22*(1+x22) + x24*(1+x24)) = 0;
cons24 : x24*(2+5*x24^2) + 1 - (x19*(1+x19) + x20*(1+x20) + x21*(1+x21) + x22*(1+x22) + x23*(1+x23) + x25*(1+x25)) = 0;
cons25 : x25*(2+5*x25^2) + 1 - (x20*(1+x20) + x21*(1+x21) + x22*(1+x22) + x23*(1+x23) + x24*(1+x24) + x26*(1+x26)) = 0;
cons26 : x26*(2+5*x26^2) + 1 - (x21*(1+x21) + x22*(1+x22) + x23*(1+x23) + x24*(1+x24) + x25*(1+x25) + x27*(1+x27)) = 0;
cons27 : x27*(2+5*x27^2) + 1 - (x22*(1+x22) + x23*(1+x23) + x24*(1+x24) + x25*(1+x25) + x26*(1+x26) + x28*(1+x28)) = 0;
cons28 : x28*(2+5*x28^2) + 1 - (x23*(1+x23) + x24*(1+x24) + x25*(1+x25) + x26*(1+x26) + x27*(1+x27) + x29*(1+x29)) = 0;
cons29 : x29*(2+5*x29^2) + 1 - (x24*(1+x24) + x25*(1+x25) + x26*(1+x26) + x27*(1+x27) + x28*(1+x28) + x30*(1+x30)) = 0;
cons30 : x30*(2+5*x30^2) + 1 - (x25*(1+x25) + x26*(1+x26) + x27*(1+x27) + x28*(1+x28) + x29*(1+x29) + x31*(1+x31)) = 0;
cons31 : x31*(2+5*x31^2) + 1 - (x26*(1+x26) + x27*(1+x27) + x28*(1+x28) + x29*(1+x29) + x30*(1+x30) + x32*(1+x32)) = 0;
cons32 : x32*(2+5*x32^2) + 1 - (x27*(1+x27) + x28*(1+x28) + x29*(1+x29) + x30*(1+x30) + x31*(1+x31) + x33*(1+x33)) = 0;
cons33 : x33*(2+5*x33^2) + 1 - (x28*(1+x28) + x29*(1+x29) + x30*(1+x30) + x31*(1+x31) + x32*(1+x32) + x34*(1+x34)) = 0;
cons34 : x34*(2+5*x34^2) + 1 - (x29*(1+x29) + x30*(1+x30) + x31*(1+x31) + x32*(1+x32) + x33*(1+x33) + x35*(1+x35)) = 0;
cons35 : x35*(2+5*x35^2) + 1 - (x30*(1+x30) + x31*(1+x31) + x32*(1+x32) + x33*(1+x33) + x34*(1+x34) + x36*(1+x36)) = 0;
cons36 : x36*(2+5*x36^2) + 1 - (x31*(1+x31) + x32*(1+x32) + x33*(1+x33) + x34*(1+x34) + x35*(1+x35) + x37*(1+x37)) = 0;
cons37 : x37*(2+5*x37^2) + 1 - (x32*(1+x32) + x33*(1+x33) + x34*(1+x34) + x35*(1+x35) + x36*(1+x36) + x38*(1+x38)) = 0;
cons38 : x38*(2+5*x38^2) + 1 - (x33*(1+x33) + x34*(1+x34) + x35*(1+x35) + x36*(1+x36) + x37*(1+x37) + x39*(1+x39)) = 0;
cons39 : x39*(2+5*x39^2) + 1 - (x34*(1+x34) + x35*(1+x35) + x36*(1+x36) + x37*(1+x37) + x38*(1+x38) + x40*(1+x40)) = 0;
cons40 : x40*(2+5*x40^2) + 1 - (x35*(1+x35) + x36*(1+x36) + x37*(1+x37) + x38*(1+x38) + x39*(1+x39)) = 0;

solve;
display 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, x5, x6, x7, x8, x9;