Adapted Falkner-type methods solving oscillatory second-order differential equations (Q1944742)

The authors extends the Falkner method with backward differences for the initial value problem of second order differential equation to solve the same problem for a system of second order differential equations of special form \[ u''(t) + M u(t) = g(t, u(t)), \] \[ u(t_{0})= u_{0},\qquad u'(t_{0})= u'_{0}, \] where \( M\) is a symmetric positive semi-definite matrix. The computing formulas of k-steps method have the usual form \[ u_{n+1}=u_{n}+ h u'_{n}+ h^{2} \sum^{k-1}_{j=0}\beta_{j}\nabla ^{j}f_{n}, \] \[ u'_{n+1}= u'_{n}+ h \sum^{k}_{j=0}\gamma_{j}\nabla ^{j}f_{n+1}. \] The convergence and stability of this method are proved, the order of error is illustrated for four model problems. The authors no indicate the concrete values of coefficients \(\beta_{j}\) and \(\gamma_{j}\). There are the other formulas to calculate these coefficients.