The present work includes an analytical investigation of a collapsing spherical star in f( R) gravity. The interior of the collapsing star admits a conformal flatness. Information regarding the fate of the collapse is extracted from the matching conditions of the extrinsic curvature and the Ricci curvature scalar across the boundary hypersurface of the star. The radial distribution of the physical quantities such as density, anisotropic pressure and radial heat flux are studied. The inhomogeneity of the collapsing interior leads to a non-zero acceleration. The divergence of this acceleration and the loss of energy through a heat conduction forces the rate of the collapse to die down and the formation of a zero proper volume singularity is realized only asymptotically.