code: def funTransEdgeY(theta, rho):
      
      return np.matrix([[1, 0], [-np.tan(theta) / rho, 1]], dtype=np.double)
docs: Fringe matrix in Y
  
      :param theta: fringe angle, in [rad]
      :param rho: bend radius, in [m]
      :return: 2x2 numpy array