Suppose $f:\mathbb{R}^2\rightarrow \mathbb{R}$ is differentiable at $(x_0,y_0)$. I want to show that
$$f(x,y) = f(x_0,y_0)+\frac{\partial f}{\partial x}\Delta x + \frac{\partial f}{\partial y}\Delta y +\epsilon (x,y)$$
where $e\rightarrow 0$ as $(x,y)\rightarrow (x_0,y_0)$.
The idea I have in mind is that of
$$f(x,y)\sim f(x_0,y) + \frac{\partial f}{\partial x}\Delta x\sim f(x_0,y_0)+\frac{\partial f}{\partial x}\Delta x + \frac{\partial f}{\partial y}\Delta y$$
The last approximation can be proven easily by regarding $f$ as a function of $y$ alone: let $f_{x_0}(y)=f(x_0,y)$ and note that $$f_{x_0}(y)=f_{x_0}(y_0)+\frac{df}{dy}\Delta y +\epsilon _y\Delta y$$ where $$\epsilon _y =\Big( \frac{f_{x_0}(y)-f_{x_0}(y_0)}{\Delta y}-\frac{df}{dy}\Big) \rightarrow 0 \ \ \text{as} \ \ \Delta y\rightarrow 0$$
Yet I can't prove the first approximation.
I would appreciate any help.