$Z_n$ elliptic Gaudin model with open boundaries
  The $Z_n$ elliptic Gaudin model with integrable boundaries specified by generic non-diagonal K-matrices with $n+1$ free boundary parameters is studied. The commuting families of Gaudin operators are diagonalized by the algebraic Bethe ansatz method. The eigenvalues and the corresponding Bethe ansatz equations are obtained. 
W. -L. Yang, R. Sasaki and Y. -Z. Zhang