Метод численного решения линейных дифференциальных уравнений

Предлагается алгоритм численного решения задачи Коши и краевой задачи для линейных обыкновенных дифференциальных уравнений. При $q$ вычислениях правых частей получается приближение к решению с ошибкой $O(\tau^{2q+1})$ на шаге. Метод определен для любого натурального $q$.