@objective(xstar=(-<NUM_LIT>, -<NUM_LIT>))<EOL>def mccormick(theta):
x, y = theta<EOL>obj = np.sin(x + y) + (x - y)**<NUM_LIT:2> - <NUM_LIT> * x + <NUM_LIT> * y + <NUM_LIT:1><EOL>grad = np.array([np.cos(x + y) + <NUM_LIT:2> * (x - y) - <NUM_LIT>,<EOL>np.cos(x + y) - <NUM_LIT:2> * (x - y) + <NUM_LIT>])<EOL>return obj, grad<EOL>
McCormick function