Approximate equation of a surface/curves to be used as an objective function

Solution 1:

Calculate values of your function on a $3 \times 3$ Grid, and then use 2D quadratic Lagrange interpolation.