Зададим на интервале 
 временную сетку 
, n = 0,1,...,N, 
с равномерным шагом 
.
Введем численную модель системы (1.1), обозначив 
приближение решения x(t) через 
.
Срезки функций x(t) и u(t) на полуинтервал 
  
будем обозначать 
 
и
, 
соответственно.
Будем также использовать обозначения 
, 
, 
, 
, 
, 
, 
, 
.
Для натурального m назовем m-этапным методом Рунге-Кутта 
( с кусочно-постоянной интерполяцией и экстраполяцией предыстории)
численную модель вида 
![]()

![]()
![]()
 
здесь 
, 
, 
 - параметры метода, 
функция-предыстория 
 при 
 кусочно-постоянна 
(интерполяция):
![]()
функция-предыстория 
 при 
 продолжена вправо 
(экстраполяция) постоянным вектором, равным 
: 

Заметим, что необходимая в методах (2.1) - (2.3) 
интерполяция (2.4) и экстраполяция (2.5) функции-предыстории 
модели 
 проводится простейшим методом, требующим минимального 
объема вычислений и памяти. 
Далее будут рассмотрены более точные способы интерполяции и экстраполяции  
функции 
.