Asymptotic-numerical approximations for highly oscillatory second-order differential equations by the phase function method