AbstractA numerical solution method based on the reduction of systems of ordinary differential equations to the Shannon form is considered. Shannon’s equations differ in that they contain only multiplication and summation operations. The absence of functional transformations makes it possible to simplify and unify the process of numerical integration of differential equations in the form of Shannon. To do this, it is sufficient in the initial equations given in the normal form of Cauchy to make a simple replacement of variables. In contrast to the classical fourth-order Runge-Kutta method, the numerical method under consideration may have a higher order of accuracy.