dy = np.append(dy, (y[i] - y[i - 2]) / 2.0)