The bifurcations from degenerate homoclinics to linearly independent subharmonic solutions are considered for singular ordinary differential equations. The scalar ∈ R + and time-dependent function g are regarded as parameters. By means of the functional analytic method based on the Lyapunov-Schmidt reduction, we obtain various bifurcation functions where the zeros correspond to the various coexistences of subharmonic solutions at nonzero parameter values. These situations of various coexistences actually define different bifurcation manifolds with finite codimension in an appropriate parameter space.