$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.