def lagrange_interpolation(x, y, x_interp): """ Interpolate y at x_interp using Lagrange interpolation.
def df(x): return 2 * x