Abstract
In this paper, we introduce a new four-stages two-step fraught in phase method for the first time in this area. A complete theoretical study is presented, and the capability of the new method is sampled using systems of coupled differential equations which represent quantum chemistry problems.