Implementation of the Taylor series method for solving ordinary differential equations

New variable-step size variable-order algorithm and software implementations of the explicit Taylor series method for solving nonstiff ordinary differential equations with polynomial right-hand sides are proposed. The version of the method being used is based on new simple formulas for the recursive computation of the Taylor coefficients of the solutions and on new rigorous a priori local error estimates combined with conventional nonstrict a posteriori considerations. The author's Fortran 95 program is compared with three other existing Fortran programs of good quality that implement Dormand–Prince, Gragg–Bulirsch–Stoer and Taylor series explicit methods, respectively. Numerical experiments proves the competitive abilities of the program, its applicability and reliability to solve real problems of dynamics.