A new multistep method with optimized characteristics for initial and/or boundary value problems

In this paper we introduce, for the first time in the literature, an optimized multistage symmetric two-step method. This method is considered as optimized due to the following reasons: (1) it is of tenth-algebraic order scheme, (2) it has obliterated the phase-lag and its first, second, third and fourth derivatives, (3) it has improved stability characteristics, (4) it is a P-stable method. For the new proposed multistage symmetric two-step method we present a full theoretical investigation consisted of: (1) local truncation error and comparative error analysis, (2) stability analysis and (3) interval of periodicity analysis. The effectiveness of the new builded multistage symmetric two-step method is evaluated on the solution of systems of coupled differential equations of the Schrödinger type.