“…However, it turned out that in the Fermion case the solutions exhibit singularities which are hard to treat numerically [3,4]. Recently, exactly solvable pairing models have gained new attention [5], with applications to nanometalic grains [6] (for a review see [7]), Bose Einstein Condensates [8] and nuclear physics [9]. It was first shown that the pairing model was integrable by finding the complete set of commuting integrals of motions [10], and subsequently, three new families of fully integrable and exactly solvable models [11] giving rise to a large class of pairing Hamiltonians with non-uniform matrix elements [12] were presented.…”