The linear approximation of a function z=f(x,y) at a point (px,py) is
z=f(x,y)≈f(px,py)+(x−px)∂z∂x(px,py)+(y−py)∂z∂y(px,py)
