  In a recent paper we studied the behavior of the pairing gaps $\Delta_F$ as a
function of the Fermi momentum $k_F$ for neutron and nuclear matter in all
relevant angular momentum channels where superfluidity is believed to naturally
emerge. The calculations employed realistic chiral nucleon-nucleon potentials
with the inclusion of three-body forces and self-energy effects. In this
contribution, after a detailed description of the numerical method we employed
in the solution of the BCS equations, we will show a preliminary analysis of
the Cooper pair wavefunctions.
