  Supersymmetric t-J Gaudin models with both periodic and open boundary
conditions are constructed and diagonalized by means of the algebraic Bethe
ansatz method. Off-shell Bethe ansatz equations of the Gaudin systems are
derived, and used to construct and solve the KZ equations associated with
$sl(2|1)^{(1)}$ superalgebra.
