It is known, that many phenomena of neutrality are reduced to solving ordinary differential equation (ODE). There are several papers dedicated to solving ODE. In this paper, which compares many known algorithms applied to solving differential equations, also suggested an algorithm that uses hybrid methods and give a procedure for constructing higher order of accuracy hybrid methods. Concrete methods are constructed, and their advantages are indicated. In particularly suggested stable methods with the order of accuracy using two mesh points (for )