from gekko import GEKKO m = GEKKO() x = m.Var(0.5) y = m.if3(x-3,6*m.exp(-x)+1,x) m.Minimize(y) m.solve() print(x.value[0],y.value[0])