Hybrid numerical methods for periodic initial value problems involving second-order differential equations