Approximate solution of the Cauchy problem for ordinary differential equations by the method of Chebyshev series

An approximate analytical method of solving the systems of ordinary differential equations resolved with respect to the derivatives of unknown functions is considered. This method is based on the approximation of the solution to the Cauchy problem and its derivatives by partial sums of shifted Chebyshev series. The coefficients of the series are determined by an iterative process with the use of Markov's quadrature formulas. This approach can be used to solve ordinary differential equations with a higher accuracy and with a larger discretization step compared to the known Runge–Kutta and Adams methods.