Piecewise linear approximation of functions of two variables in MILP models

We consider three easy-to-implement methods for the piecewise linear approximation of functions of two variables. We experimentally evaluate their approximation quality, and give a detailed description of how the methods can be embedded in a MILP model. The advantages and drawbacks of the three methods are discussed on numerical examples.