def METHOD_NAME(ct0, dt0, st0, et0, it0, logp_c, logp_d, beta, gamma, delta):
    # bt0 = trng.binomial(n=st0, p=beta)
    bt0 = st0 * beta
    bt0 = bt0.astype(st0.dtype)
    logp_c1 = binom_log_prob(et0, gamma, ct0).astype(logp_c.dtype)
    logp_d1 = binom_log_prob(it0, delta, dt0).astype(logp_d.dtype)
    st1 = st0 - bt0
    et1 = et0 + bt0 - ct0
    it1 = it0 + ct0 - dt0
    return st1, et1, it1, logp_c1, logp_d1