We present a method to calculate nuclear matter properties in the superfluid
phase. The method is based on the use of self-consistent off-shell nucleon
propagators in the T-matrix equation. Such a complete treatment of the spectral
function, is required below and around $T_c$ due to a pseudogap formation in
the spectral function. In the superfluid phase we introduce the anomalous
self-energy in the fermion propagators and in the T-matrix equation,
consistently with the strong coupling BCS equations. The equations for the
nucleon spectral function include both a contribution of condensed and
scattering pairs. The method is illustrated by numerical calculations. Above
$T_c$ pseudogap formation is visible in the spectral function and below $T_c$ a
superfluid gap also appears.Comment: correted version, appendix on numerical methods adde