mk = (y[1:] - y[:-1]) / hk