from gekko import GEKKO # define new GEKKO model m = GEKKO() # calculate x to minimize objective x = m.Var(1) # use abs to define y y = m.abs(x) # define objective to minimize m.Minimize((y+3)**2) # solve m.solve() # print solution print('x: ' + str(x.value)) print('y: ' + str(y.value))