So I know what the gradient of a (mathematical) function is, so I feel like I should know what
Return the gradient of an N-dimensional array.
The gradient is computed using central differences in the interior and first differences at the boundaries.
The default distance is 1
This means that in the interior it is computed as
where h = 1.0
and at the boundaries