Abstract
Some two-step almost P-stable methods with phase-lag of order infinite are developed for the numerical integration of second order periodic initial-value problem. One of them has algebraic order four and the other has algebraic order six. Extensive numerical testing indicates that these new methods are generally more accurate than other two-step methods.