code: def transRbend(theta=None, rho=None, gamma=None, incsym=-1):
      
      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
docs: 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