alpha = k1 * beta + k2 * poly.deriv(2)[0] + k3 * H0 + k4