Das ist eine etwas ungewöhnliche Variante des Mittelwertsatzs:
$$f(y_1,y_2)-f(x_1,x_2)=(f(y_1,y_2)-f(x_1,y_2))+(f(x_1,y_2)-f(x_1,x_2)$$
Andwendung des 1-dim-Mittelwertsatzes auf die erste Differenz liefert die Existenz eines \(z_1 \in \R\) mit
$$f(y_1,y_2)-f(x_1,y_2)=(x_1-x_1)\partial_1f(z_1,y_2)$$
Analog für die 2. Differenz.