Second Order Partial Derivatives Estimation Using a Modification of the Lagrange Quadratic Interpolant