Piecewise linear approximation of functions of two variables in MILP models