<SYSTEM_TASK:>
Solve the following problem using Python, implementing the functions described below, one line at a time
<END_TASK>
<USER_TASK:>
Description:
 def transRbend(theta=None, rho=None, gamma=None, incsym=-1):
 """ Transport matrix of rectangle dipole :param theta: bending angle in [RAD] :param incsym: incident symmetry, -1 by default, available options: * -1: left half symmetry, * 0: full symmetry, * 1: right half symmetry :param rho: bending radius in [m] :param gamma: electron energy, gamma value :return: 6x6 numpy array """

    if None in (theta, rho, gamma):
        print("warning: 'theta', 'rho', 'gamma' should be positive float numbers.")
        m = np.eye(6, 6, dtype=np.float64)
        return m
    else:
        beta12d = {'-1': (0, theta), '0': (theta * 0.5, theta * 0.5), '1': (theta, 0)}
        (beta1, beta2) = beta12d[str(incsym)]
        mf1 = transFringe(beta=beta1, rho=rho)
        mf2 = transFringe(beta=beta2, rho=rho)
        ms = transSect(theta=theta, rho=rho, gamma=gamma)
        m = reduce(np.dot, [mf1, ms, mf2])
        return m