Sector matrix in Y

    :param theta: bend angle, in [rad]
    :param rho: bend radius, in [m]
    :return: 2x2 numpy array

def funTransSectY(theta, rho):
    """ Sector matrix in Y

    :param theta: bend angle, in [rad]
    :param rho: bend radius, in [m]
    :return: 2x2 numpy array
    """
    return np.matrix([[1, rho * theta], [0, 1]], dtype=np.double)