I'm attempting to perform interpolation in 3 dimensions and have a question that hopefully someone can answer.
The derivative approximation is simple in a single direction:
df/dx(i,j,k)= [f(i+1,j,k) - f(i-1,j,k)] / 2
And I know that in the second order:
d2f/dxdy(i,j,k)= [f(i+1,j+1,k)...