## ams_version=1.0 Model Main_like { Variable x1; Variable x2; Variable x3; Variable x4 { Range: nonnegative; } Variable x5 { Range: nonnegative; } Variable x6 { Range: nonnegative; } Variable x7 { Range: nonnegative; } Variable x8 { Range: nonnegative; } Variable x9 { Range: nonnegative; } Variable objvar; Constraint e1 { Definition: { -(log(0.398942448887604*(x1/x7*exp(-0.5*sqr((95 - x4)/x7)) + x2/x8*exp(- 0.5*sqr((95 - x5)/x8)) + x3/x9*exp(-0.5*sqr((95 - x6)/x9)))) + log( 0.398942448887604*(x1/x7*exp(-0.5*sqr((105 - x4)/x7)) + x2/x8*exp(-0.5* sqr((105 - x5)/x8)) + x3/x9*exp(-0.5*sqr((105 - x6)/x9)))) + 4*log( 0.398942448887604*(x1/x7*exp(-0.5*sqr((110 - x4)/x7)) + x2/x8*exp(-0.5* sqr((110 - x5)/x8)) + x3/x9*exp(-0.5*sqr((110 - x6)/x9)))) + 4*log( 0.398942448887604*(x1/x7*exp(-0.5*sqr((115 - x4)/x7)) + x2/x8*exp(-0.5* sqr((115 - x5)/x8)) + x3/x9*exp(-0.5*sqr((115 - x6)/x9)))) + 15*log( 0.398942448887604*(x1/x7*exp(-0.5*sqr((120 - x4)/x7)) + x2/x8*exp(-0.5* sqr((120 - x5)/x8)) + x3/x9*exp(-0.5*sqr((120 - x6)/x9)))) + 15*log( 0.398942448887604*(x1/x7*exp(-0.5*sqr((125 - x4)/x7)) + x2/x8*exp(-0.5* sqr((125 - x5)/x8)) + x3/x9*exp(-0.5*sqr((125 - x6)/x9)))) + 15*log( 0.398942448887604*(x1/x7*exp(-0.5*sqr((130 - x4)/x7)) + x2/x8*exp(-0.5* sqr((130 - x5)/x8)) + x3/x9*exp(-0.5*sqr((130 - x6)/x9)))) + 13*log( 0.398942448887604*(x1/x7*exp(-0.5*sqr((135 - x4)/x7)) + x2/x8*exp(-0.5* sqr((135 - x5)/x8)) + x3/x9*exp(-0.5*sqr((135 - x6)/x9)))) + 21*log( 0.398942448887604*(x1/x7*exp(-0.5*sqr((140 - x4)/x7)) + x2/x8*exp(-0.5* sqr((140 - x5)/x8)) + x3/x9*exp(-0.5*sqr((140 - x6)/x9)))) + 12*log( 0.398942448887604*(x1/x7*exp(-0.5*sqr((145 - x4)/x7)) + x2/x8*exp(-0.5* sqr((145 - x5)/x8)) + x3/x9*exp(-0.5*sqr((145 - x6)/x9)))) + 17*log( 0.398942448887604*(x1/x7*exp(-0.5*sqr((150 - x4)/x7)) + x2/x8*exp(-0.5* sqr((150 - x5)/x8)) + x3/x9*exp(-0.5*sqr((150 - x6)/x9)))) + 4*log( 0.398942448887604*(x1/x7*exp(-0.5*sqr((155 - x4)/x7)) + x2/x8*exp(-0.5* sqr((155 - x5)/x8)) + x3/x9*exp(-0.5*sqr((155 - x6)/x9)))) + 20*log( 0.398942448887604*(x1/x7*exp(-0.5*sqr((160 - x4)/x7)) + x2/x8*exp(-0.5* sqr((160 - x5)/x8)) + x3/x9*exp(-0.5*sqr((160 - x6)/x9)))) + 8*log( 0.398942448887604*(x1/x7*exp(-0.5*sqr((165 - x4)/x7)) + x2/x8*exp(-0.5* sqr((165 - x5)/x8)) + x3/x9*exp(-0.5*sqr((165 - x6)/x9)))) + 17*log( 0.398942448887604*(x1/x7*exp(-0.5*sqr((170 - x4)/x7)) + x2/x8*exp(-0.5* sqr((170 - x5)/x8)) + x3/x9*exp(-0.5*sqr((170 - x6)/x9)))) + 8*log( 0.398942448887604*(x1/x7*exp(-0.5*sqr((175 - x4)/x7)) + x2/x8*exp(-0.5* sqr((175 - x5)/x8)) + x3/x9*exp(-0.5*sqr((175 - x6)/x9)))) + 6*log( 0.398942448887604*(x1/x7*exp(-0.5*sqr((180 - x4)/x7)) + x2/x8*exp(-0.5* sqr((180 - x5)/x8)) + x3/x9*exp(-0.5*sqr((180 - x6)/x9)))) + 6*log( 0.398942448887604*(x1/x7*exp(-0.5*sqr((185 - x4)/x7)) + x2/x8*exp(-0.5* sqr((185 - x5)/x8)) + x3/x9*exp(-0.5*sqr((185 - x6)/x9)))) + 7*log( 0.398942448887604*(x1/x7*exp(-0.5*sqr((190 - x4)/x7)) + x2/x8*exp(-0.5* sqr((190 - x5)/x8)) + x3/x9*exp(-0.5*sqr((190 - x6)/x9)))) + 4*log( 0.398942448887604*(x1/x7*exp(-0.5*sqr((195 - x4)/x7)) + x2/x8*exp(-0.5* sqr((195 - x5)/x8)) + x3/x9*exp(-0.5*sqr((195 - x6)/x9)))) + 3*log( 0.398942448887604*(x1/x7*exp(-0.5*sqr((200 - x4)/x7)) + x2/x8*exp(-0.5* sqr((200 - x5)/x8)) + x3/x9*exp(-0.5*sqr((200 - x6)/x9)))) + 3*log( 0.398942448887604*(x1/x7*exp(-0.5*sqr((205 - x4)/x7)) + x2/x8*exp(-0.5* sqr((205 - x5)/x8)) + x3/x9*exp(-0.5*sqr((205 - x6)/x9)))) + 8*log( 0.398942448887604*(x1/x7*exp(-0.5*sqr((210 - x4)/x7)) + x2/x8*exp(-0.5* sqr((210 - x5)/x8)) + x3/x9*exp(-0.5*sqr((210 - x6)/x9)))) + log( 0.398942448887604*(x1/x7*exp(-0.5*sqr((215 - x4)/x7)) + x2/x8*exp(-0.5* sqr((215 - x5)/x8)) + x3/x9*exp(-0.5*sqr((215 - x6)/x9)))) + 6*log( 0.398942448887604*(x1/x7*exp(-0.5*sqr((220 - x4)/x7)) + x2/x8*exp(-0.5* sqr((220 - x5)/x8)) + x3/x9*exp(-0.5*sqr((220 - x6)/x9)))) + 5*log( 0.398942448887604*(x1/x7*exp(-0.5*sqr((230 - x4)/x7)) + x2/x8*exp(-0.5* sqr((230 - x5)/x8)) + x3/x9*exp(-0.5*sqr((230 - x6)/x9)))) + log( 0.398942448887604*(x1/x7*exp(-0.5*sqr((235 - x4)/x7)) + x2/x8*exp(-0.5* sqr((235 - x5)/x8)) + x3/x9*exp(-0.5*sqr((235 - x6)/x9)))) + 7*log( 0.398942448887604*(x1/x7*exp(-0.5*sqr((240 - x4)/x7)) + x2/x8*exp(-0.5* sqr((240 - x5)/x8)) + x3/x9*exp(-0.5*sqr((240 - x6)/x9)))) + log( 0.398942448887604*(x1/x7*exp(-0.5*sqr((245 - x4)/x7)) + x2/x8*exp(-0.5* sqr((245 - x5)/x8)) + x3/x9*exp(-0.5*sqr((245 - x6)/x9)))) + 2*log( 0.398942448887604*(x1/x7*exp(-0.5*sqr((260 - x4)/x7)) + x2/x8*exp(-0.5* sqr((260 - x5)/x8)) + x3/x9*exp(-0.5*sqr((260 - x6)/x9))))) - objvar = 0 } } Constraint e2 { Definition: x1 + x2 + x3 = 1; } Constraint e3 { Definition: - x4 + x5 >= 0; } Constraint e4 { Definition: - x5 + x6 >= 0; } Procedure MainInitialization { Body: { x1.lower := 0.1; x2.lower := 0.1; x3.lower := 0.1; x1.l := 0.333333333333333; x2.l := 0.333333333333333; x3.l := 0.333333333333333; x4.l := 130; x5.l := 160; x6.l := 190; x7.l := 15; x8.l := 15; x9.l := 15; } } MathematicalProgram like { Objective: objvar; Direction: minimize; Constraints: AllConstraints; Variables: AllVariables; Type: NLP; } Procedure MainExecution { Body: { solve like; } } Procedure MainTermination { Body: { return 1; } } }