De Vogelaere's method for the numerical integration of second-order differential equations without explicit first derivatives: Application to coupled equations arising from the schrödinger equation