The present study aims to design a second-order nonlinear Lane-Emden coupled functional differential model and numerically investigate by using the famous spectral collocation method. For validation of the newly designed model, three dissimilar variants have been considered and formulated numerically by applying a famous spectral collocation method. Moreover, a comparison of the obtained results with the exact/true results endorses the effectiveness and competency of the newly designed model, as well as the present technique.