integrate z1 + (z2 - z1) / (1 + exp (- k (z2 - z1) x) (z2 - z1) / (y0 - 1))