In gravitational theory and astrophysical dynamics, singular initial value problems (IVPs) are frequently encountered. Finding the solutions to this class of IVPs can be challenging due to their complex nature. This study strives to circumvent the complexity by proposing a numerical method for solving such problems. The approach proposed in the current research seeks solutions to the IVP by partitioning the domain [0,L] of the problem into two intervals and solving the problem on each domain. The study seeks a closed-form solution to the IVP in the interval containing the singular point. A linearization technique and piecewise partitioning of the domain not containing the singularity are applied to the nonlinear IVP. The resulting linearized differential equation is solved using the Chebyshev spectral collocation method. Some examples are presented to illustrate the efficiency of the proposed method. Numerical analysis of the solution and residual errors are shown to ascertain convergence and accuracy. The results suggest that the technique gives accurate convergent solutions using a few collocation points.