If I had to guess: 1. Get partial derivatives of each parameter's influence on the moved object by sampling the function twice. 2. Perform a bog standard gradient descent. 3. go back to 1 until convergence is achieved.

Or you could be lazy and just throw an off-the-shelf optimizer at it.