rhs = [3 * ((f(x[i]) - f(x[i - 1])) / (x[i] - x[i - 1]) - (f(x[i - 1]) - f(x[i - 2])) / (x[i - 1] - x[i - 2]))