Fit
for that purpose, which uses the least squares method. But, in each data point, in addition to the function value, I have the function derivative value. I need the fitting algorithm to try to match those as well (squares of value differences and square of derivative differences contribute to the total sum of squares being minimized with equal weights).