I'm having a bit of trouble with my thesis. I'm doing symbolic regression of implicit mathematical functions. The approach I went with uses estimation of derivatives provided by experimentally obtained data, which help in finding the underlying function of the data.
So far I have used fitting scattered data to conic sections for data table of any dimension size. For example, if I have 2D points, the least squares fit would look like: Ax^2 + By^2 + Cxy + Dx + Ey + 1 = 0. Out of the resulted cone, I would compare its derivatives to my wanted function which gives me a sanity check that the derivatives match and that I can undergo the symbolic regression with some other algorithms.
This method has served me well for spherical functions, for example Gauss function, circle, sphere, ellipse or any kind of other spherical object.
However for example, if I apply the same fitting algorithm to estimate the derivative of Z = X * Y, the derivatives are way off for each point.
I'm coming for a more technical and computer-science approach than the mathematical one so maybe I'm missing something. Can you please suggest some derivation estimation methods for multiple dimensional data, or give any advice to it?