I am used to deriving consistency for approximations of derivatives such as
$\frac{\partial u(x,y)}{\partial x} = \frac{u(x+h,y) - u(x,y)}{h}$
and then going on with taylor series:
$u(x+h,y) = u(x,y) + hu'(x,y) + \frac{h^2}{2}u''(x,y) ... $
but how does this work for this for example:
$\frac{\partial^2 u(x,y)}{\partial x \partial y} = \frac{u(x+h,y+h) - u(x-h,y+h) - u(x+h, y-h) + u(x-h,y-h)}{4h^2}$
where I have $u(x+h,y+h)$? What would be the taylor series of this?