Here is a minimal working example:
f = FunctionInterpolation[
Evaluate[Table[D[Sin[x - y], {{x, y}, k}], {k, 0, 1}]], {x, 0.,
20. Pi}, {y, 0., 20. Pi}]
Plot3D[{f[x, y] - Sin[x - y]}, {x, 0, 20 Pi}, {y, 0, 20 Pi}]
I don't know whether this is a maths problem or MMA problem, but MMA cannot capture the fact that the function is only a function of the difference of two variables. This thing can of course still be redundantly expressed as a function of two variables, so FunctionInterpolation
should be able to handle it anyways.
(The actual example has a function which is approximately dependent only on the difference and I do care about the deviance from this behaviour.)