b[0][0] = 2.0 * ((y[1] - y[0]) / (x[1] - x[0]) ** 3