Equivalent of numpy.gradient in Julia

I’m going to add a “Theory” page to the CoupledFields documentation, but for now there is a brief explanation here: Smooth delta functions to compute Jacobian and Hessian from samples - #10 by Mattriks