Abstract A new family of three-stage two-step methods are presented in this paper. These methods are of algebraic order 12 and have an important P-stability property. To make these methods, vanishing phase-lag and some of its derivatives have been used. The main structure of these methods are multiderivative, and the combined phases have been applied for expanding stability interval and for achieving P-stability. The advantage of the new methods in comparison with similar methods, in terms of efficiency, accuracy, and stability, has been showed by the implementation of them in some important problems, including the radial time-independent Schrodinger equation during the resonance problems with the use of the Woods-Saxon potential undamped Duffing equation, etc.