Compute *func* and its derivatives along x and y at positions *xs*, *ys*, ...
returns (u, v) as a ``(2, shape)`` array.